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The  generation  of  underwater  sound  by  a  high-power  laser  pulse  is  analyzed  both 
theoretically  and  experimentally.  The  mechanism  of  sound  generation  is  the 
thermal  expansion  of  the  water  caused  by  heat  imparted  by  the  laser.  The  basic 
physical  properties  of  such  a  laser-induced  thermoacoustic  source  are  investi¬ 
gated  with  a  time  domain  approach  which  is  valid  in  the  nearfield  of  the  source. 
Results  are  then  extended  to  the  case  of  a  thermoacoustic  source  moving  at 
subsonic,  transonic,  or  supersonic  velocities  on  the  surface  of  the  water. 
Special  attention  is  paid  to  the  case  of  a  source  roving  at  a  velocity  close  to 


dd 


FORM 
JAN  73 


EOITION  Or  I  NOV  S5  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CL  ASIlFIC  ATION  OF  THIS  PAOEfWtM*  DM*  Bnt#r#«; 


UNCLASSIFIED 


security  classification  or  this  rAOcrariaa  o««  biwMO 


20.  (cont'd) 

the  speed  of  sound  in  the  water,  because,  at  such  a  velocity,  high  pressure 
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generation  more  attractive  for  practical  applications.  Experimental  results 
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Most  of  the  experimental  results  presented  in  this  study  were  obtained  with  the 
Neodymium :G1 ass  laser.  Four  kinds  of  results  were  obtained:  pressure  waveforms, 
directivity  patterns,  sound  level  dependence  on  source  velocity,  and  spreading 
curves.  In  general,  the  experimental  results  are  in  good  agreement  with  the 
theoretical  predictions. 
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CHAPTER  I 


INTRODUCTION 

The  generation  of  sound  by  a  source  of  light  was  investigated 
experimentally  more  than  one  hundred  years  ago  by  Bell,1  who  come  to  the 
conclusion  that  "sonorousness,  under  the  influence  of  intermittent  light,  is 
a  property  common  to  all  matter".  Bells  experiment  resulted  in  the 
construction  of  a  "photophone"  or  apparatus  for  the  production  of  sound  by 
light.  Bell's  invention  of  the  photophone  was  neglected  for  many  years,  but 
recently  it  has  received  more  attention  because  new  developments  in  high 
power  laser  technology  moke  this  type  of  sound  generation  more  suitable 
for  practical  applications.  There  is  an  obvious  advantage  to  opto-ocoustic 
generation  of  sound:  it  does  not  require  any  physical  transducer  in  the 
medium  in  which  one  wants  to  generate  sound.  This  unique  property  of 
laser-induced  sound  has  prompted  intensive  research  in  the  past  twenty 
years  in  various  fields  such  as  nondestructive  testing,  molecular 
spectroscopy,  and  sonar  applications  to  name  only  a  few.  Tam2  lists  more 
than  400  references  on  laser-induced  sound  and  Pierce's  recent 
bibliography3  on  the  same  subject  is  more  than  twelve  pages  long. 

There  are  several  ways  to  produce  sound  with  a  laser.  The  most 
common  mechanisms  of  optical  to  acoustical  transduction  are,  by  order  of 
increasing  efficiency,  electrostriction,  thermal  expansion,  surface 
evaporation,  explosive  boiling  and  optical  breakdown.  Electrostriction  is 


the  result  of  polarization  of  molecules  due  to  the  electromagnetic  input  in 
the  medium.  The  polarization  of  molecules  induces  changes  in  density  end 
therefore  produces  sound  which  propagates  in  the  medium.  The  thermal 
mechanism  relies  on  heating  of  the  medium  by  the  laser.  Changes  in 
temperature  produce  changes  in  density  and  subsequently  an  acoustic  wove. 
Surface  evaporation  may  occur  if  the  energy  density  of  the  laser  et  the 
point  of  impact  on  the  medium  is  sufficiently  high.  In  this  case  momentum 
transferred  into  the  medium  generates  a  sound  wave.  At  even  higher  energy 
densities,  it  is  possible  to  hove  bubble  formation  in  the  medium,  and  the 
collapse  of  the  bubbles  radiates  noise.  Finally,  for  extremely  high  loser 
energy  densities,  such  as  in  the  case  of  a  focused  high  power  laser,  a  very 
hot  plasma  is  formed  locally  around  impurities  in  the  liquid  and  ‘rapidly 
expanding  cavities  appear  in  the  focal  region,  followed  by  shock  wave 
propagation,  cavity  deceleration,  localized  cavitation,  end  eventual  bubble 
collapse.-4  This  is  known  as  sound  generation  by  optical  breakdown.  The 
earliest  experimental  detections  of  laser-induced  sound  were  reported  in 
1963. 5>6 

In  this  study,  we  restrict  our  attention  to  the  specific  mechanism  of 
thermal  expansion.  A  sound  source  which  relies  on  heating  of  the  medium  by 
a  light  source  is  referred  to  as  a  thermoacoustic  source  (TS)  or 
optoecoustic  antenna.  Most  thermoacoustic  sources  are  created  by  shining  a 
laser  into  a  medium  (usually  water).  The  intensity  of  the  laser  is  amplitude 
modulated  at  a  fixed  frequency  so  that  it  induces  a  periodic  heating  of  the 
medium  and  therefore  a  fluctuation  of  density.  This  in  turn  generates  an 
acoustic  wove  whose  frequency  is  equal  to  the  frequency  at  which  the 


laser  intensity  is  modulated.  The  wave  equation  describing  the  sound  field 
of  a  lossless  medium  containing  heat  sources  was  first  derived  by  Ingard.7 
A  few  years  later  Westervelt  end  Larson8  showed  theoretically  that  very 
directive  sound  beams  of  low  frequency  can  be  achieved  by  exploiting  the 
thermoecoustic  conversion  of  energy  in  water.  The  first  experimental 
verifications  of  the  directional  properties  of  thermoacoustic  sources  were 
made  by  Muir,  Culbertson  and  Clynch,9  who  confirmed  the  validity  of  the 
theoretical  model.  As  expected,  the  efficiency  of  the  conversion  of 
electromagnetic  energy  into  acoustic  energy  was  found  to  be  very  small  (of 
the  order  of  10'8)  and  they  concluded  that,  “for  potential  applications  of 
practical  interest,  (...)  megawatts  of  optical  power  would  probably  be 
required  in  the  megahertz  frequency  region,  with  gigawatts  of  power 
required  in  the  kilohertz  region,  and  terawetts  of  power  necessary  in  the 
low  audio  band."  This  low  efficiency  barrier  was  investigated  in  detail  by 
Soviet  physicists  and  they  showed10'11'12  that  motion  of  the  TS  is 
expected  to  significantly  increase  the  peak  amplitude  of  the  thermoacoustic 
signal,  especially  for  a  source  moving  at  velocities  close  to  the  speed  of 
sound  for  the  medium.  Such  high  source  velocities  are  easily  achieved 
experimentally.  Since  there  is  no  physical  transducer  in  the  medium,  drag 
force  end  flow  noise  are  non-existent. 

The  main  objective  of  the  present  study  is  to  investigate  the  basic 
physical  properties  of  a  laser-induced  TS  and  to  extend  the  results  to  the 
cose  of  a  moving  thermoecoustic  source  (MTS).  Previous  studies  of  MTS 
have  been  limited  by  experimental  apparatus.  They  are  also  restricted  to 
ferfield  radiation,  and  theoretical  models  generally  break  down  when  the 


MTS  is  moving  at  a  velocity  close  to  the  speed  of  sound  in  the  medium.  The 
sophisticated  experimental  apparatus  used  in  our  investigation  yielded  very 
clean  thermoacoustic  signals.  The  theoretical  study  is  based  on  a  new 
approach  which  provides  information  about  the  nearfield  of  a  TS  and  about 
the  acoustic  signature  when  the  source  is  moving  at  transonic  velocity. 

The  experimental  apparatus  is  described  in  the  second  chapter  of  this 
study.  A  detailed  analysis  of  a  stationary  TS,  both  theoretical  and 
experimental,  is  presented  in  the  third  chapter.  The  theoretical  approach  is 
to  obtain  the  impulse  response  of  a  TS.  The  time  domain  approach  is  used 
because  it  can  be  very  easily  extended  to  the  case  of  a  TS  moving  at  any 
velocity,  including  at  transonic  velocity,  ft  also  yields  some  interesting 
results  about  the  nearfield.  Diffraction  effects  due  to  the  finite  width  of 
the  laser  beam  are  also  discussed  in  this  third  chapter,  and  the  comparison 
between  theory  and  experiment  is  shown  for  pressure  waveforms, 
directivity  patterns  and  spreading  curves. 

The  fourth  chapter  deals  with  thermoacoustic  radiation  by  a  moving 
source,  both  theoretically  and  experimentally.  In  addition  to  presenting 
pressure  waveforms,  directivity  patterns  and  spreading  curves,  for 
subsonic,  transonic  and  supersonic  source  velocities,  we  investigate  the 
dependence  of  the  thermoacoustic  sound  level  on  source  velocity. 

Chapter  five  is  an  analysis  of  the  Doppler  shift  of  a  moving  source.  The 
8im  is  to  understand  certain  ferfield  approximations  often  made  in  the  case 
of  a  transonic  source.  Specifically,  we  derive  in  this  chapter  an  expression 
for  the  Doppler  shift  which  is  valid  whether  the  angle  of  observation 
between  the  source  and  the  receiver  is  fixed  or  a  function  of  time. 
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The  mein  conclusions  of  this  study  ere  given  in  the  lest  chepter  elong 
with  idees  ebout  possible  theoreticel  end  experimentel  work  thet  would  be 
interesting  to  develop  in  the  future. 
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CHAPTER  II 


EXPERIMENTAL  APPARATUS 

The  experimental  apparatus  involved  in  our  study  of  moving 
thermoacoustic  sources  consisted  basically  of  a  high  power  laser,  a  means 
of  deflecting  the  beam  so  that  motion  of  the  source  over  a  water  surface 
could  be  achieved,  and  a  receiving  and  analyzing  system.  A  block  diagram  of 
the  experimental  system  is  shown  in  Fig.  1,  and  a  description  of  the  mein 
components  of  the  system  is  given  below.  More  information  on  the 
experimental  apparatus  can  be  found  in  reference  13. 

A. Laser 

The  laser  used  in  this  study  was  designed  and  built  by  Apollo  Lasers. 
It  is  the  same  laser  described  in  reference  14  .  Figures  2  and  3,  which  ere 
taken  from  reference  14,  show  the  components  of  the  modulated  laser 
system.  The  loser  head  assembly  consists  of  an  interchangeable  rod,  either 
(Neodymium)  Nd.gless  or  ruby.  Optical  pumping  of  the  rod  is  achieved  by  a 
fleshlamp  ,  and  the  wavelength  of  the  light  emitted  by  the  laser  is  1.06  ym 
(infrared)  for  Nd:g!ess  and  0.6943  ym  (red)  for  ruby  light.  The  laser  beam 
is  polarized  in  order  to  achieve  modulation  of  its  intensity.  The  polarizing 
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element  is  e  Brewster  stack  in  the  cese  of  the  Nd:g)ess  rod  end  a 
calcite  prism  analyzer  in  the  case  of  the  ruby  rod.  The  laser  cavity  is 
formed  os  usual  by  a  100%  reflective  back  mirror  end  an  output  front 
mirror  whose  coating  depends  on  the  optical  wavelength  being  emitted. 

The  laser  is  used  in  the  conventional  non-Q-switched  mode  so  that  its 
output  is  a  series  of  very  short  (20-100  ns)  Gaussian  pulses,  characteristic 
of  the  normal  modes  distributed  statistically  in  space  and  time.  The 
envelope  of  the  laser  pulse  is  denoted  l0(t),  and  the  equivalence  between  a 
modulated  uniform  pulse  and  a  modulated  series  of  dense  random  spikes  has 
been  established  in  reference  9 

In  order  to  achieve  periodic  heating  of  the  column  of  water  being 
illuminated  by  the  laser  beam,  one  has  to  modulate  the  intensity  of  the 
laser  beam  The  modulation  is  achieved  by  a  Pockels  cell  end  a  prism 
analyzer.  The  principle  of  the  modulation  of  the  laser  is  shown  in  Fig.  4  . 
When  a  voltage  is  applied  to  the  electrodes  of  the  Pockels  cell,  it  rotates 
the  plane  of  polarization  of  the  light  passing  through  it  by  an  amount  which 
is  nonlinearly  related  to  the  applied  voltage.  The  prism  analyzer  is  a 
polarizer  with  a  fixed  plane  of  polarization  and  the  relationship  between 
the  transmitted  light  intensity  l(t)  end  the  incident  light  intensity  l0(t)  is 
given  by 


l(t)  *  l0(t)  sin2(ifV/2V1/5)  ,  (2.1) 

where  V  is  the  voltage  applied  to  the  Pockels  cell  and  V)/2  is  the 


10 


I(t),  MODULATED  MODULATED  CONTROL 

INTENSITY  VOLTAGE 


(a)  NORMAL  MODE  LASER  PULSE 


(b)  MODULATION  CONTROL  SIGNAL 


0  10  t 

TIME  -  msec 


(c)  MODULATED  NORMAL  MODE  LASER  PULSE 


FIGURE  4 

REPRESENTATIVE  OPERATION  OF  MODULATED  LASER  SYSTEM 

ARL  -  UT 
AS-73-220 
TGm  •  DR 
3-19-73 
REV4-7-'3 


half-wave  voltage  of  the  modulation  system.  Typical  pe8k  values  of  V  are 
7.0  kV  when  using  the  Ndiglass  rod  and  4.5  kV  when  using  the  ruby  rod.  If 
the  voltage  V  is  chosen  to  be  a  triangular  function  of  time  with  an 
amplitude  of  V1/2  ut/tf  ,  as  indicated  in  Fig.  4(b),  the  resulting  modulation 
yields  a  transmitted  intensity  of  the  form 

l(t)  =  0.5  l0(t)  (  1  -  cos  w0t)  (2.2) 

The  amplitude  modulation  of  the  laser  intensity  is  expected  to  achieve  a 
periodic  heating  of  the  column  of  water  being  illuminated  by  the  laser 
beam,  and  therefore  periodic  changes  in  density,  and  subsequently  sound  at 
8  frequency  f 0=g)0/2tt.  The  normalized  envelope  l0(t)  of  the  laser  intensity 
in  the  experiment  is  adequately  described  by 

l0(t)  =  (10.8  t/tp)  exp(-5.0  t/Cp)  ,  (2.3) 

where  tp  is  the  laser  pulse  duration.  The  specifications  of  the  modulated 
laser  system  are  given  in  Table  I,  reprinted  from  reference  14  .  The 
average  laser  pulse  duration  varies  between  0.75  ms  and  1.25  ms.  The  laser 
C8n  deliver  up  to  5  joules,  end  in  most  experiments  the  laser  beam  was 
focused  so  that  the  spot  diameter  on  the  surface  of  the  water  was  8bout 
1  cm.  The  triangular  signal  used  to  monitor  the  modulation  of  the  laser 
intensity  was  provided  by  an  Exact  model  126  signal  generator,  and  the 
modulation  frequency  f0could  be  varied  between  5  kHz  and  80  kHz.  The 
repetition  rate  of  the  laser  pulse  was  limited  to  four  pulses  per  minute 
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TABLE  I 

MODULATED  LASER  SYSTEM  SPECIFICATIONS 


Optical  wavelength  (pm) 

Unmodulated  output 
energy  (J)  /  power  (kW) 

Modulated  output 
energy  (J)  /  power  (kW) 

Pulse  duration  (ms) 


RUBY 

Nd.GLASS 

0.6943 

1.06 

25  J 

30  J 

25  kW 

30  kW 

5  J 

5  J 

5  kW 

5  kW 

1  ms 

1  ms 

Modulation  frequency  (kHz)  5  to  80  5  to 


deam  divergence  (mrad) 


1.7 


2.2 


because  of  the  time  constant  needed  to  charge  the  capacitors  of  the  laser 
power  supply. 


B.  Experimental  set-up 

The  experiments  were  conducted  at  Applied  Research  Laboratories, 
The  University  of  Texas  at  Austin  (ARL:UT)  in  its  18.29  m  x  4.57  m  x  3.66  m 
(60  ft  x  15  ft  x  12  ft)  fresh  water  t8nk.  Figures  5  and  6  show  the  mechanical 
arrangements  for  the  experiment.  The  laser  is  mounted  on  a  table  near  the 
water  tank.  Shock  mounts  are  used  to  prevent  mechanical  vibrations  which 
occur  after  the  laser  discharge,  from  being  transmitted  to  the  water.  After 
deflection  of  90°  by  a  tot8l  internal  reflection  prism,  the  laser  beam 
passes  through  8  pipe  to  a  rotating  mirror.  There  8re  three  major  reasons 
for  using  a  pipe  on  the  light  path.  (1)  it  is  obviously  safer  as  far  as  eye 
damage  is  concerned;  (2)  the  rotating  mirror  can  be  mounted  directly  on  the 
pipe  rather  than  on  the  ceiling  so  that  alignment  problems  are  simplified; 
and  (3)  it  provides  a  rigid  support  where  optical  lenses  can  be  installed  in 
order  to  focus  the  laser  beam  to  the  desired  size  on  the  surface  of  the 
water. 

C.  Rotating  mirror 

The  rotating  mirror  was  procured  by  Lincoln  Laser  Company.  The 
specifications  are  given  in  Table  II.  Special  attention  was  given  to 
environmental  conditions  such  as  relative  humidity  (77%  in  the  tank  room) 
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FIGURE  6 

MECHANICAL  ARRANGEMENT  -  II 


TABLE  II 

SPECIFICATIONS  FOR  THE  ROTATING  MIRROR  SVSTEM 


Mirror  substrate  materiel. 
Number  of  facets: 
Circumscribed  circle  diameter: 
Facet  width: 

Facet  length: 

Surface  flatness. 

Surface  reflectivity: 

Speed  of  rotation: 

Linear  source  speed: 


aluminum  and  thorium  flouride 
6 

15.2  cm  (6  in.) 

3.2  cm  (1.25  in.) 

7.6  cm  (3  in.) 

1/4  optical  wavelength 
82  *  at  0.6943  pm 
88  %  at  1.06  pm 
150  -  3500  rpm  +  0.1  % 
x  Mach  0.1  to  Mach  2.6 


end  dust  accumulation.  The  rotational  speed  of  the  mirror  was  adjustable 
between  150  and  4500  rpm  so  that  source  velocities  on  the  water  surface 
could  be  varied  between  M  »  0.1  and  M  «  2.6  ,  where  M  is  the  Mach  number  of 
the  source  with  respect  to  the  sound  speed  in  water.  The  stability  of  the 
angular  velocity  of  the  rotating  mirror  was  given  to  be  ±  0.1  %  between  150 
and  3500  rpm,  i.e.,  between  M  »  0.1  and  M  »  2.  In  this  design  the  source 
velocity  on  the  water  surface  is  not  truly  constant  because  the  angle  the 
beam  makes  with  the  water  is  not  constant, i.e.,  the  laser  beam  is  not 
always  perpendicular  to  the  water  surface.  It  is  estimated  that  the  error 
Av/v  ,  where  v  is  the  source  velocity,  is  less  than  3  %  at  Mach  1  ,  and 
decreased  with  lower  Mach  number.  It  was  also  estimated  that  the 
maximum  tilt  angle  of  the  thermoacoustic  source  in  water  was  less  than  8° 
at  Mach  1  and  decreased  with  lower  Mach  number.  This  estimate  takes  into 
account  the  refraction  of  the  laser  beam  at  the  air-water  interface.  A 
detailed  discussion  is  given  in  Appendix  A,  where  it  is  shown  that  these 
secondary  effects  can  safely  be  ignored. 

D.  Synchronization 

The  firing  of  the  laser  must  be  synchronized  with  the  angular 
position  of  the  rotating  mirror  so  that  the  laser  beam  is  deflected  to  a 
known  and  repeatable  position  on  the  surface  of  the  water.  One  way  to  know 
the  angular  position  of  the  rotating  mirror  is  to  use  a  light  emitting  diode 
(LED)  and  a  phototrensistor,  both  combined  in  a  smell  reflective  transducer. 
When  the  light  beam  emitted  by  the  LED  is  perpendicular  to  a  facet  of  the 


mirror,  it  is  reflected  to  the  phototrensistor  and  an  electrical  impulse  of  a 
few  microseconds  duration  is  produced.  When  the  mirror  rotates  at 
constant  speed,  the  reflective  transducer  delivers  a  series  of  short 
impulses  giving  an  indication  of  the  angular  position  of  the  mirror.  Any  one 
of  these  impulses  may  be  delayed  in  time  and  used  to  trigger  the  firing  of 
the  laser.  A  digital  circuit  has  been  designed  and  built  in  order  to  compute 
the  required  time  delay  to  trigger  the  laser  at  the  appropriate  time.  A 
photograph  of  the  laser  rotation  firing  controller  circuit  is  shown  in  Fig.  7. 

The  manually  controlled  settings  in  the  firing  control  circuit  are  a 
set  of  three  dip  switches  (DS).  The  input  to  the  circuit  is  a  pair  of  analog 
signals  provided  by  two  LEO/phototransistor  reflective  transducers  located 
in  the  rotating  mirror  housing.  Outputs  of  the  firing  control  circuit  include 
signals  from  connectors  labeled  "T.P.A.",  "OUTPUT”,  end  "FIRE".  The  fire 
signal  is  a  pulse  of  approximative^  100  ys  duration  and  amplitude  of 
12  volts  to  trigger  the  discharge  of  the  capacitors  used  to  fire  the  loser. 
The  T.P.A  signal  is  a  series  of  impulses  delivered  by  the 
LED/phototransistor.  It  can  be  used  to  accurately  check  the  mirror  angular 
velocity  and  its  stability.  The  OUTPUT  signal  is  a  single  pulse  of  about  1  ms 
duration,  slightly  delayed  after  the  FIRE  pulse.  The  amount  of  the  time 
delay  is  manually  controlled  by  the  dip  switch  DS3  (see  Fig.  7  ).  The 
reason  for  this  time  delay  is  that  there  is  an  intrinsic  time  interval 
between  the  discharge  of  the  capacitor  bank  end  the  actual  lasing  of  the 
rod.  This  time  delay  is  about  180  ys  for  the  Nd.glass  rod  and  340  ys  for  the 
ruby  rod.  Calibration  curves  for  values  of  the  binary  number 
associated  with  the  dip  switch  DS3  controlling  the  time  delay  ,  indicate 
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that  NDS3  should  be  101  when  using  the  Nd:Gless  rod  end  1010  when  using  the 
ruby  rod.  The  OUTPUT  signal  is  intended  to  be  used  to  trigger  the  intensity 
modulator  external  signal  provided  by  the  Exact  signal  generator. 

The  firing  control  circuit  contains  two  other  dip  switches,  DS1  end 
DS2,  that  ere  manually  controllable.  These  are  indicated  in  Fig.  7.  The 
decimal  numbers  associated  with  0S1  and  DS2  are,  respectively,  Npg,  and 
Npsj.  These  integer  numbers  are  used  to  select  the  initial  angle  ^  between 
the  vertical  and  the  laser  beam  (see  Fig.  5)  at  the  beginning  of  the  laser 
pulse.  The  initial  angle  is  given  by 

120°  (Npsj/N,^)  -32°  +  (n  120°)  ,  (2.4) 

where  the  "plus"  sign  is  to  be  used  when  the  laser  beam  is  moving  away 
from  the  bench  and  similarly,  where  the  "minus"  sign  has  to  be  used  when 
the  beam  is  moving  towards  the  bench.  In  Eq.  (2.4),  n  is  a  positive  integer 
accounting  for  the  periodicity  of  the  facets  on  the  rotating  mirror.  As  an 
example,  if  one  wants  to  use  the  Nd:glass  rod  in  such  a  way  that  the  initial 
angle  ^  is  2ero,  Npgj  should  reed  "101",  NDS1  should  read  "100"  (i.e.,  4  in 
decimal),  8nd  Np52  should  read  "1111"  (i.e.,  15  in  decimal). 

E.  Receiving  system 

The  signal  emitted  by  the  thermoecoustic  source  is  received  by  a 
hydrophone,  amplified  end  filtered,  before  being  stored  end  analyzed  on  a 
digital  oscilloscope  which  is  interfaced  with  a  mainframe  computer  for 


more  advanced  signal  processing. 

The  hydrophone  used  in  most  experiments  described  in  this  study  is  a 
USRD  type  H-56  hydrophone.  This  hydrophone  is  very  sensitive  (~  -  163  dB 
re  1  V/yPa  at  60  kHz)  but  its  bandwidth  is  limited  to  frequencies  below 
100  kHz.  Therefore  in  some  cases  an  H-23  hydrophone  was  preferred. 
Calibration  curves  for  the  frequency  response  and  directivity  patterns  of 
the  hydrophone  (H-56)  at  several  frequencies,  are  given  in  Appendix  B. 

The  signal  received  by  the  hydrophone  was  amplified  by  two  battery 
operated  Burr-Brown  model  100  AC-decade  amplifiers,  with  a  calibrated 
total  gain  of  40.5  dB.  The  amplified  signal  was  bandpass  filtered  by  a 
Krohn-Hite  filter.  The  signal  was  then  transferred  to  a  Nicolet  4094  digital 
oscilloscope.  This  oscilloscope  provides  two  input  channels  (expandable  to 
four  channels)  with  simultaneous  12  bit  digitizers  that  sample  at  a  rate  of 
2  MHz.  The  data  is  stored  on  floppy  diskettes  in  the  external  disk  recorder. 
Special  diskettes  are  available  from  Nicolet  for  performing  valuable  signal 
processing  functions  on  the  oscilloscope.  For  example,  the  rms  value  of  a 
signal  over  a  specified  time  interval  can  be  calculated.  Similarly,  programs 
to  compute  Fourier  transforms,  digital  filtering,  integration,  etc.  are 
available.  In  general  the  sound  level  measured  from  an  experimental 
waveform  was  obtained  from  the  peek-to-peek  values  of  the  pressure.  Also, 
the  Nicolet  4094  was  interfaced  to  the  mainframe  CYBER-CDC  computer  at 
ARLUT  vie  an  RS-232  interface,  so  that  experimental  data  could  be 
transferred  to  the  computer  and  compered  with  theoretical  predictions. 
The  software  required  for  such  a  transfer  consists  essentially  of  two 
programs,  whose  binary  version  ere  called  “NICOLT"  end  "CNVRT".  The 


program  NICOLT  creates  a  file  on  the  CYBER  that  contains  integer  numbers 
representing  the  signal  being  transferred  from  the  Nicolet  to  the  CYBER. 
Program  CNVRT  is  used  to  convert  these  numbers  into  their  real  values. 
Note  that  NICOLT  can  also  be  used  to  transfer  waveforms  from  the  CYBER 
to  the  Nicolet.  More  information  on  the  procedure  for  data  transfer  is  given 
in  Appendix  C. 

In  general  the  sound  pressure  level  measured  from  an  experimenel 
waveform  was  obtained  by  taking  20  times  the  logarithm  based  10  of  the 
ratio  of  the  peek-to-peak  pressure  to  a  reference  pressure  of  1  yPa. 
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CHAPTER  HI 


STATIONARY  THERMOACOUSTIC  SOURCES15 

The  sound  field  redieted  by  e  stetionery  thermoecoustic  source  is 
enelyzed  in  this  chepter  by  introducing  the  impulse  response  of  the 
thermoecoustic  system.  This  epproech  hes  four  mein  edventeges:  (1)  it  cen 
eesily  be  extended  to  the  cese  of  e  TS  moving  et  eny  velocity  through  the 
medium,  end  this  will  be  done  in  Chepter  IV,  (2)  the  model  is  velid  in  the 
neerfield  of  the  source,  (3)  it  cen  be  used,  in  principle,  for  eny  spetiel  end 
temporel  leser  intensity  profiles,  end  (4)  trensforms  ere  not  needed  to  get 
time  informetion. 

In  the  first  section  of  this  chepter  we  describe  the  physicel  problem 
of  thermoecoustic  redietion  end  derive  the  impulse  response  of  e  TS.  The 
second,  third,  end  fourth  sections  deel,  respectively,  with  the  pressure 
weveform,  the  directivity,  end  the  spreeding  cherecteristics  of  the 
redietion  from  e  TS,  both  theoreticelly  end  experimentelly.  The  lest  section 
is  e  summery  of  the  importent  feetures  of  this  chepter. 

A.  Impulse  response 

Consider  the  thermoecoustic  source  shown  in  Fig.  8.  A  leser  beem  is 
shined  et  normel  incidence  in  the  +z  direction  into  en  optical ty 
ebsorbing  medium.  The  heet  produced  by  the  leser  beem  generetes  elmost 
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instantaneously  8  density  disturbance  due  to  thermal  expansion  which 
induces  a  pressure  fluctuation  which  is  observed  as  8n  acoustic  wave.  The 
objective  of  this  study  is  to  relate  in  the  most  general  wsy  the  spatial  and 
temporal  characteristics  of  the  laser  to  the  acoustic  pressure  observed  at 
any  point  in  the  field  . 

Larson16  h8s  shown  that  the  wave  equation  describing 
thermoacoustic  radiation  is,  in  the  C8se  of  an  inviscid  8nd 
non-heet-conducting  fluid,  given  by 
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V2p  -  c~2  pu  =  -  (p/Cp)  qt  ^  (3.1) 

where  the  subscript  t  denotes  a  time  derivative,  p  is  the  acoustic  pressure 
generated  by  the  thermal  expansion  mechanism,  c  is  the  small  signal  sound 
speed,  Cp  is  the  specific  heat  of  the  medium  at  constant  pressure,  p  is  the 
logarithmic  coefficient  of  thermal  expansion  of  the  medium,  and  q  is  the 
he8t  energy  per  unit  volume  per  unit  time  8dded  to  the  medium.  It  is  related 
to  the  laser  intensity  vector  I  by  the  relation 

q  =  -  V-  <l(x,y,z,t)  >  (3.2) 


A  derivation  of  the  wave  equation  for  thermal  generation  of  sound  in  a 
viscous  medium  is  given  in  Appendix  D.  For  practical  reasons,  we  consider 
the  medium  to  be  a  liquid  having  8  constant  optical  coefficient  of 
absorption  oc,  so  that  the  z-component  of  the  laser  beam  intensity  1 
decreases  with  depth  according  to  the  exponential  shading  exp(-ocz).  With  a 
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different  absorption  lew,  the  following  analysis  would  still  be  workable, 
although  more  complicated  and  not  necessarily  more  realistic  . 

When  Eqs.  (3.1)  and  (3.2)  are  combined  it  is  clear  that  the 
thermoacoustic  source  function  is  proportional  to  the  time  derivative  of 
the  laser  intensity.  It  is  therefore  convenient  to  define  the  impulse 
response  h(t)  of  the  thermoacoustic  system  in  the  following  manner. 


p(t)  =  h(t)  *  ft(t) 


(3.3) 


where  the  asterisk  denotes  convolution  over  time.  In  this  formulation  h(t) 
contains  ell  of  the  spatial  characteristics  of  the  system  at  any  arbitrary 
range  r,  and  1^(0  contains  all  of  the  temporal  information  of  the  source. 
The  spatial  characteristics  of  the  system  can  be  decomposed  into  its 
vertical  component  (x-z  plane)  and  its  horizontal  component  (x-y  plane).  Let 
us  denote  by  l\(t)  and  h0(t),  respectively,  the  vertical  and  horizontal 
impulse  responses  of  the  TS.  Then  h|_(t)  represents  the  response  of  an 
infinitely  narrow  TS  (a-»0  in  Fig.  8)  of  effective  length  L=l/cc.  Thus  hL(t)  is 
the  response  of  an  exponentially  tapered  semi-infinite  line  source  below  a 
pressure  release  boundary  such  as  the  air-water  interface.  On  the  other 
h8nd  he(t)  represents  the  response  of  an  infinitely  short  TS  (L-»0).  In  other 
words,  h0(t)  is  the  impulse  response  of  a  cross-section  of  the  source.  With 
these  notations  Eq.  (3.3)  can  be  rewritten  as 


p(t)  =  hL(t)  *  h0(t)  *  lt(t) 


(3.4) 


It  is  therefore  necessary  to  evaluate  hL(t)  and  h0(t)  in  order  to  find  the 
pressure  response  of  a  given  medium  excited  by  a  laser  beam  of  intensity 
l(t).  This  analysis  is  carried  out  below. 

1.  Vertical  impulse  response  hj_(t) 

The  vertical  impulse  response  is  obtained  by  imagining  that  the 
laser  beam  diameter  is  so  small  that  sound  coming  from  an  entire  beam 
cross  section  S  is  received  instantaneously.  This  has  the  effect  of  setting 
the  horizontal  impulse  response  in  Eq.  (3.4)  to  a  delta  function  in  time:  h0(t) 
->  S  S(t).  In  this  case  the  heat  q  is  related  to  the  laser  intensity  by 

q  =  -  d/dz  (  A  l(t)  e"  *  z  ]  S  8(x-xc)  8(y-yc)  ,  (3.5) 

where  A  is  the  coefficient  of  transmission  of  light  across  the  interface, 
and  where  xc  and  yc  denote  the  center  coordinates  of  the  beam.  The 
solution  of  Eq.  (3.1)  for  a  pressure  release  boundary  condition  at  the 
interface  is  then  obtained  by  integrating  the  Green's  function  of  the 
problem  over  the  source  volume  as  follows. 

CO 

p(t)  =  C1  (r)~l  e"**2  l^(t-r/c)dz 
o 

o 

-Ci  f  (rT1  e_ccz  It(t-rVc)  dz 


(3.6) 


where  the  coordinates  z,  r,  and  r*  are  shown  in  Fig.  8,  and  where 
C|*(ApccSc/4trcp).  The  second  integral  in  Eq.  (3.6)  represents  the 
contribution  of  the  mirror  image  above  the  pressure  release  boundary. 
From  Eqs.  (3.4)  and  (3.6),  then,  the  vertical  impulse  response  is  given  by 

oo 

hL(t)=(C|/S)  |  (r)"1  e-cCZ  8(t-r/c)  dz  -  mirror  image  ,  (3.7) 
o 

where  8(t)  is  the  Dirac  delta  function.  Shifting  the  integration  variable 
from  z  to  r,  yields 

1^(0  »K  J  (r)“*  e_a:z  8(r-ct)  (dz/dr)  dr  -  mirror  image  ,  (3.8) 

r 

where  K  -  (Apccc/4ircp).  Equation  (3.8)  provides  a  means  for  evaluating  the 
vertical  impulse  response  of  a  TS  once  the  relationship  between  the  depth  z 
of  the  receiver  and  its  radius  of  observation  r  is  known.  Assuming  that  the 
laser  beam  strikes  the  medium  at  normal  incidence,  one  may  divide  the 
source  into  four  regions.  These  are  indicated  in  Fig.  9.  For  each  region  the 
relationship  between  z  and  r  can  easily  be  found.  Equation  (3.8)  can  then  be 
written  in  terms  of  separate  integrals  over  the  four  regions  as  follows. 


© 


FIGURE  9 

STATIONARY  THERMO  ACOUSTIC  SOURCE  GEOMETRY  -  II 


A«L  UT 
AS-84  782 
YHB  GA 

9  5  84 


30 


hL<t) «  K  J  (1/r)  S(r-ct)  expl-ff(rocos0o-s)l  (-r/s)  dr 
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CD 


-  K  |  (1/r)  8(r-ct)  expl-oc(~rocos0o+s))  (r/s)  dr 


(3.9) 


where  s  «=  V  r2-(rQsin0o)2 . 


The  four  integrals  in  Eq.  (3.9)  represent  the  contributions  to  the 
acoustic  signal  from  the  four  source  regions  shown  in  Fig.  9  However,  it 
is  more  meaningful  to  distinguish  the  array  response  by  regions  in  time 
rather  than  in  space.  The  limits  of  integration  indicate  that  the  analysis 
of  the  Impulse  response  should  be  divided  into  five  time  regimes  as 
follows. 


o)  t<  rosin0D/c:  In  this  time  region,  the  impulse  response  reduces  to 


hL(t)  =  0  , 


(3.10) 


which  states  that  the  acoustic  emission  from  the  array  has  not  yet  had 
time  to  reach  the  receiver.  The  first  disturbance  emitted  by  the  TS  has  to 
travel  a  distance  r0sin90  before  reaching  the  observation  point.  Therefore 
the  impulse  response  is  identically  zero  for  any  time  t<r0sin90/c  . 

b)  t=r0s1n®0/c.  This  particular  instant  corresponds  to  the  time  of 
arrival  of  the  first  wavelet  at  the  receiver.  This  wavelet  originated  at  a 
distance  of  r0sin90  from  the  receiver,  that  is  to  say,  at  a  depth 
z0=r0cos90  below  the  pressure  release  boundary.  Consequently  the 
impulse  response  is  simply,  from  Eq.  (3.7), 


h^(t)  =  K  (r0sin90)-1  exp  (-ctr0cos90) 


(3.11) 


c)  rosin0o/c  <  t  <  r0/c:  From  Fig.  9  ,  it  is  obvious  that  the  first  wavelet 
which  will  reach  the  receiver  emanates  from  point  E.  An  instant  later, 
wavelets  emitted  from  points  B  and  C  will  arrive  simultaneously  at  the 
receiver.  The  progression  continues  with  wavelets  eventually  coming 
from  points  A  end  D.  These  are  all  the  wavelets  received  during  the 
interval  r0sin90/c<t<r0/c  which  resulted  from  sound  radiation  at  time 
t*0.  The  disturbances  coming  from  regions  2  and  3  of  the  TS  therefore  add 
by  pairs  at  the  observation  point.  The  only  non-zero  contribution  to  the 


impulse  response  is  given  by  the  first  two  integrals  in  Eq.  (3.9),  which 
may  be  evaluated  to  give 


2K  exp(-ccr0cos#0) 


d)  t=r0/c  :  There  are  two  points  on  the  array  that  8re  located  at  a  distance 
r0  from  the  receiver.  These  two  points  are  denoted  by  A  and  D  in  Fig.  9. 
However,  the  pressure  release  boundary  condition  at  the  free  surface  of  the 
liquid  implies  that  the  acoustic  contribution  from  point  A  is  null,  so  that  in 
fact  the  acoustic  response  at  time  t*r0/c  is  simply  given  by  the  wavelet 
coming  from  D,  i.e.,  from  the  point  of  the  source  that  is  located  8t  a  depth 
z=2rocos0o.  Equation  (3.7)  yields  an  expression  for  the  impulse  response  at 
t=r0/c: 

hL(t)  =  Kr0_1  exp(-2ccrocos0o)  (3.13) 

e)  t>r0/c:  The  disturbances  emitted  by  regions  1  and  4  of  the  TS  will  arrive 
by  pairs  at  the  observation  point.  In  this  case,  however,  the  boundary 
condition  at  the  air-water  interface  implies  that  signals  coming  from 
region  1  (the  mirror  image  region)  are  inverted.  Hence  the  only  non-zero 
contribution  to  the  impulse  response  is  given  by  the  last  two  integrals  in 
Eq.  (3.9).  These  may  be  evaluated  to  give 


The  vertical  impulse  response  of  a  TS  is  therefore  given  by  Eq.  (3.15a-e). 


0 

K  e"r  O^sin^)"1 

hL(t)  -  <  2Ke“ry'1  cosh(ccy) 
K  e_2r  {r0V] 
-2Ky~1e“cf^  sinh(D 


for  1  <  t| 

(3.15e) 

for  t  -  ti 

»  / 

(3.15b) 

for  tj  <  t  <  t0 

(3.15c) 

fort-t0 

(3.15d) 

for  t  >  t0 

(3.15e) 

where  t0  -  r0/c  ,  tj  -  rosin0o/c  ,  r  -  ccrocos0o  ,  and  p  -  c*v  t2-tj2  . 

Figure  10  shows  how  the  shape  of  h|_(t)  changes  os  the  nondimensionel 
parameter  r  varies  from  0.02  to  100.  r  con  be  regarded  os  the  ratio  of 
the  depth  of  the  receiver  to  the  effective  length  L  ■  cT1  of  the  TS.  It  is 
assumed  in  Fig.  10  that  the  laser  is  shined  into  fresh  water  where  c*1486 
m/s.  Note  that  no  forfield  assumptions  were  required  to  derive  the 
analytical  expressions  above.  Figure  11  is  a  three  dimensional 
representation  of  the  impulse  response  h^t)  os  a  function  of  time  and  r. 


FIGURE  10 

IMPULSE  RESPONSE  h.  (t) 


Figures  10  and  11  show  a  curious  behavior  for  hj_(t)  when  tj<t<t0 
During  this  time  interval  the  shape  of  hj_(t)  can  change  drastically  from  a 
monotonically  decreasing  function  having  a  maximum  at  time  t^*r0sind0/c, 
to  a  monotonically  increasing  function  having  a  maximum  at  time  t0=r0/c. 
This  can  be  explained  by  looking  at  the  two  different  effects  determining 
the  shape  of  the  impulse  response:  the  exponential  shading  along  the  axis  of 
the  array  which  is  due  to  the  absorption  of  light  in  the  medium,  and  the 
spherical  spreading  associated  with  each  wavelet  radiated  by  the  TS.  For 
simplicity  let  us  assume  first  that  the  exponential  decay  along  the  depth  of 
the  array  is  much  more  important  than  the  spherical  spreading.  This  is 
expected  to  be  the  case  in  the  farfield  because  the  distance  from  each  point 
from  the  TS  to  the  receiver  is  almost  the  same.  Also  let  zA,zB,zc  and  ^  be 
the  depths  of  the  elementary  sources  located  at  points  A,B,C  end  D  on  the 
source  as  shown  in  Fig.  9.  It  can  be  shown  that,  for  any  positive  tf. 


e  *  e  -tfz»  <  e  "*zb  ♦  e 


(3.16) 


Therefore  the  strength  of  the  acoustic  disturbance  coming  from  points  B 
and  C  is  always  less  important  than  the  strength  of  the  disturbance 
originated  from  points  A  and  D.  Since  the  latter  arrives  at  the  receiver  at 
time  t0,  the  impulse  response  h^t)  is  expected  to  increase  gradually  from 
ti  to  t0 

Let  us  now  assume  that  the  spherical  spreading  loss  associated  with 
each  wavelet  coming  from  the  array  is  much  more  important  then  the 
exponential  shading  along  the  z-exis  of  the  array.  In  this  case,  all  of  the 


elementary  sources  will  hove  approximately  the  same  strength. 
This  case  applies  to  the  neerfield  of  a  TS.  The  strongest  disturbance  is 
therefore  the  one  which  suffers  the  least  spherical  spreading  loss. 
Consequently  it  originates  from  the  elementary  source  on  the  array  which 
is  the  closest  to  the  receiver.  Since  the  minimum  distance  between  the 
source  end  the  receiver  is  rosin0o,  the  vertical  impulse  response  hj_(t) 
exhibits  a  peak  at  time  t^=r0sinfi0/c  and  then  decreases  gradually  until 

Wc 

It  is  also  interesting  to  note  that  the  impulse  response  l\(t)  of  8 
stationary  thermoacoustic  source  is  very  closely  related  to  the  half-order 
derivative  operator.  This  is  discussed  in  Appendix  E. 

2.  Horizontal  impulse  response  ha(t) 

The  horizontal  impulse  response  is  obtained  by  imagining  that  the 
optical  absorption  by  the  fluid  is  so  high  that  the  laser  beam  penetrates 
only  a  small  amount  compared  to  tne  beam  diameter.  Under  this  condition, 
the  vertical  impulse  response  is  set  to  a  delta  function  :  hj_(t)  -»  (K/c)  S(t), 
end  the  horizontal  impulse  response  may  be  calculated.  Note  from  Eq.  (3  3) 
that  the  time  derivative  operates  on  the  intensity,  and  from  Eq.  (3.7),  thst 
the  mirror  image  of  the  source  over  the  interface  has  already  been  t8ken 
into  account  in  the  vertical  impulse  response,  including  the  constant  of 
proportionality.  This  means  that  it  would  be  redundant  to  again  consider 
these  features  in  ha(t).  The  horizontal  impulse  response  is  thus  simply 


(3.17) 


(3.18) 


Let  us  first  examine  the  case  of  a  uniform  beam  of  radius  a.  In  this 
case  l|y|(x,y)  is  equal  to  unity,  and  following  Morse’s  analysis17  of  the 

impulse  response  of  a  piston,  one  finds,  if  r0»a  sinG0  , 
a  sin0o 

hg(t)«=  j  2  V  (a  sin  0Q)2-y2  8  [t-(r+y)/cj  dy  ,  (3.19) 

-a  sine0 

or 


h0(t)  =  2  c 


(a  sin  9q)2  -  (rQ-ct)2 


(3.20) 


if  (rQ-e  sin  #0)/c  <  t  <  (rQ+e  sin60  )/c,  end  zero  otherwise.  Similarly,  in 
the  case  of  a  Gaussian  beam  where  a  denotes  the  1/e  radius,  one  would  find, 
provided  that  r0»e  sin$0  , 


h0(t)  *  ir1/2  c  o  sin0o  exp  (  -{(r0-ct)/(e  sin0o)]2] 


(3.21) 


if  (rQ-e  sin$0)/c  <  t  <  (r0+e  sin$0)/c,  end  zero  otherwise.  It  con  be  verified 
thet 

Lim  (  h0(t  j  «  if  a2  8(t-rQ/c)  (3.22) 

•”*0*0**° 

Equations  (3.20)  end  (3.21)  provide  expressions  for  the  horizontol 
impulse  response  in  analytical  form  for  either  a  uniform  or  a  Gaussian 
intensity  shading  across  the  laser  beam.  The  spatial  extent  of  the  source 
induces  a  frequency  filtering  which  may  be  found  by  taking  the  Fourier 
transform  of  ha(t).  For  a  uniform  intensity  distribution  across  a  section  of 
the  laser  beam,  the  frequency  filtering  has  a  Bessel  characteristic  of  the 
form  2Jj(ko  sin0o)/(ka  sin$0)  end  for  the  case  of  the  Gaussian  distribution, 
it  takes  the  Gaussian  form  exp  i-(ko  sin0o/2)2l.  The  frequency  filtering 
associated  with  diffraction  due  to  finite  laser  beamwidth  is  therefore  a 
low  pass  filtering,  or  a  smoothing  in  the  time  domain. 

B.  Pressure  waveform 

In  this  section  we  present  results  obtained  using  the  impulse 
response  technique  previously  described  to  predict  the  pressure  waveform 
radiated  by  a  TS.  The  double  convolution  of  Eq.  (3.4)  is  evaluated 
numerically.  It  is  then  used  to  evaluate  the  case  of  an  unmodulated  TS.  The 
results  in  the  ferfield  limit  ere  compared  with  analytical  predictions 


derived  by  Lyamshev  and  Sedov18  in  their  review  paper  on  thermoacoustic 
sources.  The  analytical  predictions  are  8lso  compared  to  experimental 
results. 

Consider  a  thermoacoustic  source  of  the  sort  discussed  previously. 
For  simplicity  we  assume  that  the  laser  intensity  is  not  amplitude 
modulated  and  that  the  finite  beamwidth  effects  are  negligible,  i.e.,  that 
h8(t)  -♦  S  8(t).  To  facilitate  comparisons  with  experimental  data,  we  will 
assume  that  the  laser  pulse  shape  matches  the  measured  pulse  shape,  which 
is  adequately  described  by  Eq.  (2.3).  Using  Eqs.  (3.22)  and  (3.15)  in  Eq.  (3.4) 
leads  to  8  straightforward  numerical  computation  of  the  pressure  received 
at  a  specified  location  as  a  function  of  time. 

Figure  12  shows  the  evolution  of  the  predicted  pressure  waveform 
received  in  fresh  water  at  the  observation  point  as  the  nondimensionel 
parameter  r  changes  from  200  to  0.2.  The  position  of  the  receiver  is  kept 
constant  (r0  ■  4  m,  $0  =  60°)  so  that  the  variation  of  r  corresponds  in  fact 
to  a  variation  of  the  effective  length  L  ■  cT1  of  the  source  from  1  cm 
(short  source)  to  10  m  (long  source).  The  pulse  duration  was  set  to  250  ys. 
The  sharp  peak  in  the  acoustic  response  observed  in  Fig.  12(a)  is  due  to  the 
fact  that  the  laser  pulse  described  by  Eq.  (3.22)  has  a  discontinuous  time 
derivative  at  t=0.  The  negative  pressure  pulse  which  appears  for  all  values 
of  T  in  Fig.  12(e)-(f)  after  the  retarded  time  r0/c  can  be  explained  by  the 
effect  of  the  pressure  release  interface  between  air  and  water.  As  shown 
in  the  figures,  the  pressure  waveform  broadens  and  changes  substantially 
in  shape  as  the  source  length  increases.  In  the  case  of  a  very  long  source 
the  pressure  waveform  is  seen  to  begin  at  time  tj=  r0  sin0o/c  rather  then 
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the  time  t0  =  r0/c  as  it  would  in  the  case  of  a  short  source.  This  result 
agrees  with  the  analysis  presented  in  Section  III.A.1. 

Lyamshev  end  Sedov18  have  shown  that  in  the  farfield,  the  pressure 
waveform  radiated  by  a  TS  is  proportional  to  the  second  time  derivative  of 
the  laser  pulse  intensity  provided  the  length  of  the  array  is  very  small 
compered  to  a  typical  acoustic  wavelength.19  The  reason  the  pressure 
response  is  proportional  to  the  second  time  derivative  of  the  laser 
intensity  is  explainable  in  terms  of  the  impulse  response  approach.  As 
shown  by  Fig.  10,  the  impulse  response  of  a  very  short  source  tends  to 
behave  like  the  time  derivative  of  a  delta  function  centered  at  time  t0.  The 
physical  cause  of  this  behavior  is  the  pressure  release  boundary  between 
the  air  and  water.  From  Eq.  (3.3)  it  then  follows  that  the  pressure  is 
proportional  to 

$t  (t-r0/c)  *  lt(t)  *  ltt(t-r0/c)  (3.23) 

Figure  13(a)  shows  a  comparison  between  the  impulse  response 
approach  8nd  Lyamshev  snd  Sedov’s  predictions  for  the  sound  pressure  of  a 
short  thermoacoustic  source.  The  comparison  was  made  for  a  modulation 
frequency  of  5  kHz,  a  laser  pulse  duration  of  1  ms,  end  for  r=lOO.  The  laser 
intensity  was  described  by  Eqs.  (2.2)  end  (2.3).  As  shown  in  the  figure,  the 
numerical  and  analytical  predictions  for  the  short  source  case  are  in 
excellent  agreement. 

The  case  of  a  long  and  narrow  TS  has  also  been  studied  by  Lyamshev 
and  Sedov18  in  the  special  case  of  farfield  radiation.  Their  results  show 


that  under  these  conditions  the  pressure  waveform  appears  to  be  simply 
proportional  to  8n  inverted  laser  pulse.  20  As  indicated  by  Fig.  13(b), 
the  convolution  technique  produces,  as  predicted,  a  pressure  proportional 
to  8n  inverted  laser  pulse. 

An  experimental  investigation  of  thermoacoustic  sources  was 
carried  out  in  the  fresh  water  tank  at  ARL:UT  .  The  lasing  element  of  the 
optical  system  was  a  NeodymiunrGlass  rod  producing  infrared  light  at  1.06 
pm  in  the  conventional  non-Q-switched  mode.  At  this  optical  frequency  the 
absorption  of  light  in  water  W8S  measured  to  be  13.7  Np/m.  Determinations 
of  the  optical  coefficient  of  absorption  cc  are  reported  in  Appendix  F.  More 
information  on  the  laser  apparatus  is  given  in  Chapter  II. 

Figure  14(a)  shows  a  typical  acoustic  signal  recorded  on  the  digital 
oscilloscope.  This  waveform  was  obtained  under  the  following  conditions: 
r0  =  0.5  m  ;  90  =  75°  ;  a  =  0.5  cm  ,  tf  =  13.7  Np/m  ;  f0  =  35  kHz  ( f0=u0/2ff); 
tp  =  1.2  ms,  r  =1.77.  Figure  14(b)  shows  the  waveform  as  predicted  by  the 
theory  described  previously.  In  the  computer  prediction,  a  Gaussian 
intensity  distribution  across  the  laser  beam  was  assumed.  The  beam  radius 
was  very  small  compered  to  the  main  acoustic  wavelength  X0=c/f0,  and  thus 
the  laser  beamwidth  had  little  effect  on  the  resulting  pressure  waveform. 
Figure  14  shows  that  reasonable  agreement  between  the  theory  and  the 
experiment  is  obtained.  A  better  fit  could  easily  be  obtained  by  altering  the 
envelope  constants  for  the  laser  intensity  (see  Eq.  (2.3)).  This  did  not  seem 
appropriate,  since  these  constants  have  been  determined  by  averaging 
the  envelope  over  several  laser  pulses.  The  discrepancy  is  attributed  to  the 
fact  that  the  laser  pulse  envelope  was  not  very  repeatable.  It  was  also 
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noted  that  the  rise-time  of  the  experimental  pressure  waveform  was  often 
shorter  than  the  predicted  one.  It  shows  again  the  inherent  non-stationerity 
of  the  laser  pulse. 

The  acoustic  peek  signal  shown  in  Fig.  14(a)  was  measured  to  be 
122.8  dB  re  1  yPa.  This  level  is  obtained  by  taking  20  times  the  logarithm 
based  10  of  the  peak  to  peak  pressure  divided  by  the  reference  pressure.  It 
also  takes  into  account  the  directional  characteristics  of  the  hydrophone 
(see  Appendix  B),  the  gain  of  the  amplifier,  end  the  loss  through  the 
bandpass  filter  (30-40  kHz).  The  predicted  peak  level  can  be  found  from 
Fig.  14(b)  to  be  123.1  dB  re  1  yPa.  This  remarkable  agreement  between 
theory  end  experiment  is  certainly  fortuitous  because  it  was  found  that  the 
average  difference  between  experimental  levels  and  theoretical  levels  is 
8bout  +4.6  dB.  In  general,  our  numerical  predictions  underestimate  the 
acoustic  levels  obtained  experimentally,  but  are  within  a  reasonable  margin 
if  we  consider  that  the  intrinsic  properties  of  the  laser  pulse,  such  as  its 
shape,  its  energy  and  its  duration,  vary  significantly  from  pulse  to  pulse. 

C.  Directivity  patterns 

The  impulse  response  described  previously  provides  a  means  of  finding 
the  directivity  characteristics  of  a  thermoacoustic  source.  This  is  due  to 
the  fact  that  all  the  directional  characteristics  of  the  system  are 
embedded  in  the  impulse  response  of  the  system,  and  therefore  one  can 
relate  the  directivity  pattern  to  the  impulse  response  by  a  simple  Fourier 
transform.  Let  us  denote  by  H(u)  the  frequency  response  of  the  optoacoustic 


system.  Then,  using  FT  to  denote  the  Fourier  transform. 


H(y)  =  FT  l  h(t)]  «  FT  [  h8(t)  *  hL(t)  ]  =  He(y)  HL(y)  ,  (3.24) 

where  H0(y)  end  Hj_(y)  denote,  respectively,  the  Fourier  transforms  of  the 
horizontal  and  vertical  impulse  responses  h0(t)  and  h^t).  The  directivity 
pattern  at  a  given  frequency  y  is  then  simply 

I!  H  (*;y)  II 

D(0)=  20  Log,  0 -  '  <325> 

II  H  (0N;y)  II 

where  0^  yields  the  .ximum  value  of  H(0)  at  the  frequency  y,  and  the 
double  vertical  bars  indicate  the  modulus  of  a  complex  quantity.  The  Fourier 
transform  of  h^t)  is  obtained  directly  from  Eq.  (3.21),  and  in  the  case  of  a 
Gaussian  beam,  ft  gives 

H0(y)  »  S  exp  ( -  (ka  sin  0o)2/4  J  ,  (3.26) 

where  k  *y/c  is  the  acoustic  wave  number.  Similarly  one  can  find  the 
Fourier  transform  of  h^(t)  from  Eq.  (3.15)  .  It  is  convenient  to  use  the 
identity21 

t-r/c 

\  f(0)  l(t-0)2-(r/c)2]‘,/2  d0  ■  J  f|t-(r  cosh  u)/cl  du  (3.27) 

-oo  o 

By  way  of  Eq.  (3.27),  the  Fourier  transform  of  hi  (0  can  be  expressed  as 


Hl(cj)  =  Kj  \  cosht'f  sinh  u)  exp(-ir)  cosh  u)  du 
0 

CO 

'  <2  \  exp(-^  sinh  u)  expHq  cosh  u)  du  ,  (3.28) 

“o 

where  Kj  =  Apcr  e“r/2ircp 

<2  ■  Apec  sinh(D  /  2ffcp 
u0  *  Jin  I(1+cos0o)/sin0o] 

T  =  ccr0sine0 
r  =  tfrocos0o 
r)  =  krosin0 

end  we  sre  considering  only  angles  0  such  that  0  <  9  <  ir/2.  This 
expression  is  not  edsy  to  integrate  analytically  but  may  be  easily 
evaluated  numerically.  In  most  practical  ceses  the  beam  radius  is  very 
small  compared  to  an  acoustic  wavelength  so  that  h^t)  -»  S  S(t)  or 
He(u)  -» S  is  a  good  approximation. 

Figures  15(a)-(c)  show  the  result  of  the  numerical  computation  of  the 
directivity  pattern  of  a  TS  at  three  radii  of  observation:  r0  of  1  m  ,  2  m, 
and  3  m  (i.e.,  for  values  of  r/cosfi0  of  5,  10,  and  15).  The  modulation 
frequency  was  set  to  10  kHz  and  the  optical  absorption  coefficient  was 
5  Np/m  so  that  the  effective  length  of  the  array  was  20  cm.  Figure  15(e) 
shows  that  sidelobes  exist  in  the  nearfield  of  a  TS.  As  the  receiver  moves 
into  the  farfield,  the  sidelobe  structure  tends  to  disappear  and  eventually 
the  computed  directivity  reduces  to  the  usual  farfield  analytical  form 
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described  by  Lyamshev  and  Sedov:10 


where  d  »  k/d  is  a  measure  of  the  length  L  of  the  array  in  terms  of  an 
acoustic  wavelength  X. 

Figure  16  shows  the  predicted  directivity  patterns  for  different 
wavelengths  when  the  radius  of  observation  and  the  source  length  ere  held 
constant.  The  figure  shows  that  by  decreasing  the  number  of  wavelengths 
along  the  length  of  the  source  (i.e.,  for  d  decreasing) ,  one  decreases  the 
number  of  sidelobes  and  increases  the  angle  of  tilt  of  the  main  lobe  away 
from  the  interface.  Also  the  half-power  beamwidth  angle  increases 
drastically  so  that  the  end  result  is  a  dipole  type  of  directivity  as 
expected  from  a  very  small  source  (relative  to  a  wavelength)  radiating  on 
a  pressure  release  interface. 

Figure  17  shows  the  computed  directivity  when  the  length  of  the 
source  varies  while  the  distance  of  observation  r0  and  the  modulation 
frequency  f0  are  kept  constant.  Since  we  have  assumed  an  exponential  law 
of  absorption  of  light  in  water,  with  an  optical  coefficient  of  absorption 
oc  (see  Eq.  (3.5)),  the  quantity  1/ct  represents  the  depth  at  which  the 
intensity  is  attenuated  by  a  factor  1/e.  Again  the  number  of  sidelobes  and 
the  angle  of  maximum  emission  decrease  with  d,  whereas  the  half-power 
beamwidth  increases  when  d  decreases.  It  is  also  interesting  to 
note  that  the  neerfield  directivity  is  not  simply  dictated  by  the 
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FIGURE  17 

DIRECTIVITY  OF  A  STATIONARY  THERMOACOUSTIC  SOURCE  -  III 
(rQ  ■  1  m,  f„  *  10  kHz) 


nondimensional  parameter  d=k/cc,  as  in  the  fsrfield  case,  but 
independently  by  k  and  c( . 

Figure  18  shows  the  nearfield  effects  on  the  half-power  beamwidth 
as  a  function  of  distance  r0  and  plotted  for  several  array  lengths.  As 
mentioned  previously,  fyp  decreases  with  6.  For  large  values  of  r0,  0^ 
reduces,  as  expected,  to  its  farfield  estimate,  which  can  easily  be 
obtained  from  Eq.  (3.29).  In  the  farfield,  0^,  reduces  to 

0HP  =  cos“,[(/2-1)/d)  -  cos_1[(/2+1)/(5}  (3.30) 

It  is  also  interesting  to  note  that,  for  a  given  source  length,  0^  goes 
through  a  minimum  at  a  certain  distance  r0  opt.  This  may  be  of  practical 
interest  for  applications  where  angular  resolution  is  more  important  than 
long  range  propagation. 

Figure  19  shows  how  the  nearfield  effects  the  angle  of  maximum 
radiation  0max,  for  different  array  lengths,  as  a  function  of  distance. 
Again,  for  large  values  of  r0  ,  0max  reduces  to  its  farfield  limit,  which 
can  be  evaluated  from  Eq.  (3.29)  to  be 

0max  ■  cos_1(l/<5)  (3.31) 

Note  that  Figures  18  and  19  can  be  used  to  define  8  farfield  criterion. 

Experimental  data  were  obtained  for  8  modulation  frequency  of 
35  kHz,  with  o (  ■  13.7  Np/m,  at  a  radius  of  observation  of  0.50  m.  (Seme 
experiment  as  the  one  reported  in  Fig.  14).  Infrared  light  of  optical 


wavelength  of  1.06  ym  was  shined  into  fresh  water  (c* I486  m/s). 
The  pulse  length  was  of  the  order  of  1  ms.  The  receiver  was  an  H-56 
hydrophone  of  sensitivity  -171  dB  re  1  V/yPa  at  35  kHz.  The  vertical 
directivity  pattern  of  the  transducer  H-56  (see  Appendix  B)  was  taken 
into  account  in  the  analysis  of  the  directivity  of  the  TS.  The  signal  was 
amplified  and  filtered  before  being  stored  in  a  digital  oscilloscope.  For 
each  waveform,  the  rms  pressure  was  computed  numerically  end  yielded  a 
value  for  the  sound  pressure  level  in  decibels.  The  rms  pressures  were 
averaged  over  several  pulses.  The  comparison  between  predictions  based 
on  Eqs.  (3.24H3.28)  and  experiment  is  shown  in  Fig.  20.  Good  agreement 
between  theory  and  experiment  shows  again  that  the  impulse  response 
approach  is  very  well  suited  for  describing  the  neerfield  of  a 
thermoecoustic  source. 

D.  Spreading  curves 

It  is  known  that,  in  general,  the  sound  level  dependence  on  distance, 
or  spreading  curve,  is  just  a  function  of  the  geometrical  dimensions  of 
the  acoustic  source  in  terms  of  the  acoustic  wavelength.  Thus,  the  sound 
level  dependence  on  distance,  denoted  by  Ly(r),  can  be  found  at  a  given 
frequency  u  by  taking  the  Fourier  transform  of  the  impulse  response  of 
the  TS.  This  can  be  expressed  in  a  fashion  similar  to  Eq.  (3.25), 

il  H(u,r)  II 


L^r)  ■  20  Log,0 


II  H(u,r)  II  r.r 


(3.32) 


where  H(u,r)  is  the  Fourier  transform  of  1^(0,  which  is  valid  in  the 
nearfield. 

The  approach  for  finding  Lw(r)  is  thus  similar  to  the  one  described  in 
the  previous  section  for  finding  the  neerfield  directivity  Du(fl). 
Experimental  results  presented  in  Fig.  21  show  the  spreading  characteristic 
of  a  stationary  thermoacoustic  source.  The  parameters  associated  with 
Fig.  21  8re  0O  s  84°  ,  f0  =  35  kHz,  end  oC=  13.7  Np/m.  In  both  theory  and 
experiment,  the  sound  level  seems  to  decrease  at  a  rate  of  3  dB  per 
doubling  of  distance  until  a  range  of  about  1.5  m,  after  which  the  slope 
increases  to  a  value  of  6  dB  per  doubling  of  distance.  The  reason  for  the 
change  of  slope  is  that,  in  the  nearfield,  the  TS  appears  to  be  essentially  a 
cylindrical  source,  whereas  in  the  ferfield,  spherical  spreading  obviously 
predominates.  The  transition  between  cylindrical  and  spherical  spreading 
may  also  be  used  as  a  practical  estimate  of  the  Rayleigh  distance  of  the 
thermoecoustic  source.  Another  way  to  see  the  transition  is  to  look  at  the 
directivity  patterns  in  Fig.  15.  As  r0  increases,  the  major  lobe  moves 
toward  84°.  Thus  the  tendency  toward  a  6  dB  slope  is  pertly  neutralized  by 
the  progression  up  the  major  lobe  toward  its  peak.  This  argument  implies 
that  at  same  angles  the  slope  might  exceed  6  dB  in  the  neerfield. 

E.  Conclusions 

The  sound  field  radiated  by  a  thermoacoustic  source  has  been 
investigated  both  theoretically  and  experimentally.  The  theory  is  based  on  a 
time  domain  approach.  It  shows  that  the  acoustic  pressure  radiated  by  the 
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FIGURE  21 

SPREADING  OF  A  STATIONARY  THERMO  ACOUSTIC  SOURCE 
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thermal  mechanism  is  a  convolution  of  the  impulse  response  of  the  system 
with  the  thermoacoustic  source  strength  (time  derivative  of  the  laser 
intensity).  Although  the  convolution  has  to  be  performed  numerically  in 
most  practical  cases,  it  offers  several  advantages:  (1)  it  can  be  used  to 
investigate  the  case  of  a  moving  thermoacoustic  source,  (2)  it  is  valid  even 
in  the  neerfield  of  the  source,  (3)  although  we  have  restricted  our  analysis 
to  the  most  applicable  case  of  a  laser  beam  of  Gaussian  cross-section 
intensity  distribution  with  an  exponential  shading  along  the  beam 
penetration  axis,  the  model  presented  here  can  easily  be  extended  to  any 
laser  intensity  temporal  or  spatial  profile,  and  (4)  this  approach  does 
not  require  transforms  to  get  time  information.  Furthermore,  the  impulse 
response  of  a  thermoacoustic  array  may  provide  insight  into  several 
thermoacoustic  transient  phenomena  such  as  extremely  short  laser  pulses, 
shock  waves  induced  by  a  laser  source  moving  at  transonic  velocity, 
pressure  waveforms  radiated  by  lightning,  or  even  cosmic  rays  entering 
the  ocean. 

The  validity  of  the  theoretical  model  W8S  tested  by  an  experimental 
investigation  of  thermoacoustic  pressure  waveforms  and  directivity 
patterns.  Results  seem  to  indicate  that  in  the  nearfield,  the  acoustic  levels 
may  be  quite  different  from  those  predicted  by  a  farfield  theory,  and  that 
sidelobes  cannot  be  neglected  in  a  nearfield  directivity  pattern.  In  general, 
experimental  results  were  in  support  of  the  theoretical  predictions  in  both 
farfield  and  nearfield  cases. 


CHAPTER  IV 


MOVING  THERMOACOUSTIC  SOURCES22 


In  the  previous  chapter,  we  analyzed  in  detail  the  acoustic  field 
radiated  by  a  stationary  thermoacoustic  source.  The  objective  of  this 
chapter  is  to  extend  these  results  to  the  case  of  a  moving  source.  As 
mentioned  earlier,  the  major  reason  for  studying  moving  thermoacoustic 
sources  (MTS)  is  that  they  are  expected  to  produce  substantially  higher 
acoustic  levels  when  the  source  is  moving  at  a  velocity  close  to  the  speed 
of  sound  in  water. 

In  the  first  section  of  this  chapter,  we  present  the  theory  that  will 
be  used  to  predict  the  properties  of  the  acoustic  field  of  on  MTS,  even  in  the 
neorfield  of  the  source.  The  theory  is  based  on  the  impulse  response 
described  in  Section  III.A. 

Then,  in  the  second  section  of  this  chapter,  we  proceed  to  analyze 
successively  pressure  waveforms,  directivity  patterns,  sound  pressure 
level  dependence  on  source  velocity,  and  spreading  curves.  Comparisons 
between  theory  end  experiments  are  made  for  subsonic,  transonic,  and 
supersonic  source  velocities. 

The  lest  section  of  this  chapter  summarizes  the  main  conclusions  of 
this  chapter. 


A.  Theory 


1.  Pseudo-convolution 

The  pressure  signal  radiated  by  the  moving  thermoecoustic  source  (MTS) 
may  be  constructed  in  a  three-step  procedure  as  follows.  First,  the  motion 
of  the  source  is  decomposed  in  time.  As  in  8  movie,  a  "picture"  is  t8ken  at 
constant  time  intervals  At,  showing  the  laser  beam  at  different  positions 
along  its  path  during  the  laser  pulse  duration  tp  .  The  geometry  associated 
with  an  MTS  is  shown  in  Fig.  22  .  Second,  for  each  of  these  positions 
(defined  by  the  subscript  i  in  Fig.  22)  the  impulse  response  hj(t)  and  the 
corresponding  elementary  pressure  response  Pj(t)  are  evaluated 
numerically  from  the  analysis  presented  in  Section  III.A.  Third,  the  total 
pressure  pj(t)  received  at  the  hydrophone  is  obtained  by  adding  all  the 
elementary  acoustic  responses  Pj(t)  with  the  suitable  time  delays 
corresponding  to  the  motion  of  the  source.  Note  that  this  approach  is  in 
principle  also  suitable  for  a  source  moving  with  a  nonuniform  velocity 
along  its  path.  However  in  the  following  analysis  it  is  assumed  for 
simplicity  that  the  source  is  moving  at  a  constant  velocity  v.  This  model 
also  assumes  that  the  laser  beam  remains  always  perpendicular  to  the 
water  surface.  This  assumption  is  discussed  in  Appendix  A. 

The  first  step  in  the  analysis  is  to  discretize  in  time  and  space 
the  motion  of  the  source  so  that  the  problem  can  be  analyzed 
numerically.  Given  that  tp  is  the  laser  pulse  duration  and  At  the  time 
increment,  the  motion  of  the  source  is  represented  by  N  pictures  taken 
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every  At  seconds,  such  that  N=tp/At  .  For  a  source  moving  rectilinearly  at 
constant  velocity  v,  i  At  represents  the  time  taken  for  the  laser  to  move 
from  the  origin  to  the  position  shown  on  the  i-th  picture.  The  second  step  is 
to  compute  for  each  of  these  positions  :he  elementary  acoustic  pressure 
radiated  by  the  source  when  it  W8S  at  the  position  shown  on  the  i-th 
picture,  that  is  to  S8y,  after  a  time  delay  i  At.  The  laser  source  intensity  at 
that  instant  is  simply  !(i  At)  8nd  the  impulse  response  hj(t.)  C8n  be 
evaluated  from  Eq.  (3.15).  However  it  is  important  to  realize  that  the 
impulse  response  is  a  function  of  the  position  of  the  source  because  the 
distance  and  the  angle  between  the  source  and  the  receiver  are  time 
dependent  in  the  C8se  of  8  moving  source.  Let  us  denote  by  r^  and  these 
coordinates.  It  can  easily  be  shown  that 

1/2 

rj  =  [  r02  -  2  r0  (v  i  At)cos80  sin$0  +  (v  i  At)2  ]  (4.1) 

and 

0j  =  cos"1  (r0  cos  0Q  i  q)  ,  (4.2) 

where  r0,  f?0  and  4>0  refer  to  values  at  the  origin  of  the  motion  of  the 
source  (t=0).  The  impulse  response  h-j(t)  is  therefore  given  by  Eq.  (3.15) 
where  the  coordinates  r0  and  0O  have  been  replaced  by  rj  and  using 
Eqs.  (4.1)  and  (4.2). 

The  elementary  acoustic  pressure  Pj(t)  is  then  computed  as  follows. 
Equation  (3  3)  shows  that  the  optoecoustic  source  strength  is  proportional 
to  the  time  derivative  of  the  laser  intensity.  The  i-th  picture  therefore 
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shows  a  line  source  excited  by  an  impulse  of  strength  l‘(i  At)  during  a  time 
At ,  and  thus  the  elementary  acoustic  pressure  pj(t)  is  given  by 

Pi<t)  =  l'(  i  At)  fyU)  At  ,  (4.3) 

where  the  prime  denotes  a  time  derivative. 

The  third  step  is  to  add  ell  the  elementary  acoustic  pressures  Pj(t) 
radiated  during  the  motion  of  the  source,  taking  into  account  the  time  delay 
i At  taken  by  the  source  to  travel  from  the  origin  (t=0)  to  the  i-th  position. 
This  gives  the  total  acoustic  pressure  pj(t)  radiated  by  the  thermoacoustic 
source  during  its  motion: 

N 

PjU)  =  2  Pt<t-iAt)  At  (4.4) 

i=0 

Combining  Eqs.(4.3)  snd  (4.4)  gives 

N 

PT(t)  =  2  1  ('At)  hj(t-iAt)  At  (4.5) 

i*0 

Equation  (4.5)  shows  that  the  total  pressure  received  at  the  observation 
point  is  a  convolution  type  summation  in  the  time  domain  between  the 
"unsteady"  (changing  shape  with  time  and  distance)  impulse  response  ty(t) 
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and  the  optoacoustic  source  strength  dl/dt  .  Note  that  in  the  case  of  a 
stationary  source,  the  shape  of  the  impulse  response  remains  constant  8nd 
the  acoustic  pressure  pj(t)  is  a  true  convolution  between  the  impulse 
response  end  the  time  derivative  of  the  laser  intensity. 

In  the  above  results  laser  beemwidth  effects  have  been  neglected. 
They  become,  however,  very  important  when  the  acoustic  wavelength  is  of 
the  order  of  the  laser  beam  diameter,  and  this  is  always  the  case  when  the 
source  is  moving  at  velocities  close  to  transonic.  In  order  to  have  a  model 
valid  for  any  source  velocity,  it  is  therefore  necessary  to  take  into 
account  laser  beamwidth  effects.  This  is  the  object  of  the  next  section. 

2.  Laser  beamwidth  effects 

It  was  shown  in  Chapter  III  (see  Eq.  (3.4))  th8t  laser  beamwidth 
effects  can  be  accounted  for  by  8  convolution  involving  the  impulse 
response  ha(t)  of  8  cross  section  of  the  laser  beam.  This  impulse  response 
is  defined  as  the  time  response  at  the  receiver  to  an  impulse  of  heat 
applied  at  a  specified  depth  and  over  an  area  corresponding  to  the  laser 
beemwidth.  If  r0  »  a  ,  where  a  is  the  laser  beam  radius,  then 

h0(t)  =  ir,/2  c  a  sin0o  exp  {  -[  (r0-ct)/(a  sin80)  ] 2  }  ,  (3.21) 

where  we  have  assumed  that  the  intensity  distribution  across  a  section  of 
the  leser  beam  is  Gaussian.  This  last  equation  has  been  derived  previously 
in  Chapter  III.  The  effect  of  the  finite  beamwidth  of  the  laser  may  be 
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thought  of  as  introducing  many  infinitesimally  thin  line  sources  into  each 
"picture"  of  the  source  at  time  tj.  This  effect  may  be  included 
mathematically  by  convolving  the  cross-section  impulse  response  with 
Pj(t).  The  result  is  an  induced  smoothing  of  the  pressure  waveform  that  is 
8  low-pass  filtering  in  the  frequency  domain.  The  frequency  filtering  H0(u) 
associated  with  laser  beamwidth  effects  may  be  obtained  by  taking  the 
Fourier  transform  of  h0(t).  It  is  necessary  to  note,  however,  that  in  order 
to  account  for  the  motion  of  the  source,  the  main  wave  number  k0=u0/c  of 
the  acoustic  radiation  must  be  replaced  by  its  Doppler  shifted  version 
k(j  =  k0/D  =  k0/|1-Msin$0cos4»0l,  where  D  is  known  as  the  Doppler  factor. 
(The  sin#0cos$0  term  in  the  Doppler  shift  correction  is  due  to  the  fact  that 
the  receiver  sees  8n  apparent  velocity  that  is  different  from  the  actual 
velocity  of  the  source  which  has  Mach  number  M=v/c.  This  was  first  noticed 
by  Doppler23  in  1842).  The  result  is 

H0(u)  =  it  82  sin20o  exp  (  -(k^  8  sin#0/2)  (4.6) 

It  should  be  understood  that  in  the  case  of  a  moving  source,  the 
coordinates  of  the  source  r0  and  are  time  dependent,  so  that  in  fact 
Eq.  (4.6)  is  an  approximation  to  the  real  frequency  filtering  associated  with 
finite  laser  beamwidth  effects.  Since  these  effects  are  usually  small,  it 
seems  unnecessary  to  get  an  exact  and  much  more  complicated  expression 
for  He(<j),  and  therefore  Eq.  (4.6)  is  expected  to  yield  good  results  provided 
that  finite  laser  beamwidth  effects  are  relatively  small.  When  the  source 
is  moving  at  velocities  close  to  the  speed  of  sound  in  the  fluid,  the  Doppler 
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shifted  wavelength  characteristic  of  the  acoustic  radiation  shrinks  to  a 
V8lue  very  small  compared  to  the  beam  radius  a  ,  thus  generating  very 
significant  diffraction  effects,  in  this  case  the  approximation  that 
diffraction  effects  are  small  breaks  down.  The  finite  beamwidth  effects 
near  transonic  velocity  will  be  discussed  later  in  more  detail. 

The  time  domain  approach  used  to  predict  the  pressure  waveform 
radiated  by  an  MTS  is  perfectly  suitable  for  a  numerical  computation  and 
the  theoretical  predictions  obtained  with  this  model  will  be  discussed  in 
the  next  section. 

3.  Numerical  predictions 

Some  typical  predicted  pressure  waveforms  radiated  by  8  moving 
thermoacoustic  source  are  presented  in  Fig.  23.  These  waveforms  were 
obtained  by  computer  simulation,  using  the  impulse  response  approach 
described  8bove.  The  computer  program  used  for  these  simulations  is  given 
in  Appendix  G.  These  predictions  were  made  for  some  realistic  values  of  the 
important  parameters:  r0  =  4  m  ,  0O  =  75°  ,  *0  =  0°  ,  oc  =  13.7  Np/m  (Nd:Gl8ss 
lasing  wavelength  of  1.06  pm  in  fresh  water),  8=0.5  cm,  f0=7  kHz,  tp=1  ms, 
and  c  =  1486  m/s.  The  Mach  number  of  the  source  M  =  v/c  was  varied  from  0 
to  1.5. 

Figure  23(a)  shows  the  pressure  waveform  observed  at  the  receiver 
when  the  source  is  stationary.  There  is  a  periodic  nature  to  the  response 
due  to  the  modulation  of  the  intensity  of  the  laser  pulse  end  8  roughly 
exponential  decay  in  the  amplitude  due  to  the  decay  in  the  laser  intensity. 
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NUMERICAL  PREDICTIONS  OF  MTS  PRESSURE  WAVEFORMS 

rQ  =  4  m,  0Q  =  75°, <*>Q  =  0°,  a=  13.7  Np/m,  fn  =  7  kHz 


Figure  23(b)  shows  the  acoustic  response  when  the  source  is  moving 
at  M=0.5.  Clearly  there  is  a  time  compression  of  the  acoustic  signal 
observed  at  the  receiver,  i.e.,  an  increase  in  amplitude  and  frequency.  This 
is  a  typical  effect  of  the  Doppler  shift  which  occurs  for  moving  sources. 
One  important  advantage  of  the  time  domain  approach  used  in  the  numerical 
computation  of  the  pressure  waveform  is  that  the  Doppler  shift  is 
implicitly  taken  into  account  by  the  time  delays  of  the  pseudo-convolution 
(see  Eq.  (4  5));  therefore  the  Doppler  shift  appears  naturally  and  it  is  valid 
even  if  the  source  is  moving  at  transonic  velocity.  A  detailed  analysis  of 
the  Doppler  shift  associated  with  a  moving  source  of  finite  length  may  be 
found  in  reference  24.  To  a  first  approximation  the  Doppler  factor  D  seen 
by  the  receiver  in  Fig.  23(b)  is  (l-Msin#0cos*0)  *  0.517  so  that  the  Doppler 
shifted  frequency  in  Fig.  23(b)  is  about  fd  =  f0/D  =  13.5  kHz.  This  is  in  good 
agreement  with  the  Doppler  shifted  frequency  of  13.6  kHz  measured  from 
Fig.  23(b). 

Figure  23(c)  shows  the  pressure  waveform  observed  at  the  receiver 
when  the  source  is  moving  at  8  velocity  such  that  Msin0ocos*o*  1  , 
i.e.,  so  that  the  receiver  perceives  the  source  to  be  moving  at  Mach  one.  In 
this  case,  the  Doppler  factor  D  goes  to  zero  and  subsequently  the  Doppler 
shifted  frequency  goes  to  infinity  so  that  the  Doppler  shifted  wavelength 
goes  to  zero.  This  means  that  in  the  transonic  regime,  diffraction  effects 
due  to  the  finite  width  of  the  laser  beam  will  be  important.  Fig.  23(c) 
shows  that  the  received  pulse  is  extremely  short  (of  the  order  of  the 
transit  time  2e/c  across  the  laser  beam  diameter)  and  there  is  therefore 
no  periodic  structure  in  the  pressure  waveform  but  only  a  strong  inverted 


pulse  whose  shape  depends  mainly  on  the  intensity  distribution  across  a 
section  of  the  laser  beam.  This  is  discussed  in  more  detail  in  Appendix  H. 
Note  also  that,  as  predicted  by  Eq.  (4.6),  laser  beamwidth  effects  limit  the 
amplitude  of  the  pressure  response  when  the  source  is  moving  at  transonic 
velocity.  Previously  reported  results  which  describe  the  pressure  field  of 
an  MTS  in  terms  of  the  pressure  field  of  a  stationary  source  compensated 
by  a  Doppler  shift  D  ■  (1-Msin0ocos4o)  predict  an  infinite  amplitude  for 
a  source  moving  at  Mach  one  because  they  hsve  not  included  the  fact  that 
the  angle  from  the  source  to  the  receiver  is  location  dependent,  and  hence 
have  used  a  poor  approximation  to  the  Doppler  factor. 

Finally  it  is  interesting  to  note  in  Fig.  23(d)  that  the  time  axis  seems 
to  be  inverted.  This  effect  is  sometimes  referred  to  as  time  inversion  end 
it  occurs  when  a  source  moves  toward  the  receiver  faster  than  the 
disturbances  it  generates.  In  such  an  event,  the  first  sound  wavelets  to  be 
observed  by  the  receiver  will  be  those  excited  last.  This  effect  has  been 
known  for  more  th8n  a  century.  Rayleigh25,  for  instance,  explains  how  a 
musical  piece  played  at  Mach  two,  can  be  heard  "in  tune,  but  played 
backwards".  Note  also  that  the  predicted  Doppler  shifted  frequency  in 
Fig.  23(d)  is  16.0  kHz  which  is  in  good  agreement  with  the  Doppler  shifted 
frequency  predicted  by  f<j  =  f0/D  =  15.6  kHz. 

B.  Experimental  results 

The  experimental  results  presented  in  this  section  are  grouped  into 
four  categories:  pressure  waveforms,  directivity  patterns,  sound  level 


dependence  on  source  velocity,  and  spreading  curves.  The  experimental 
procedure  was  described  in  Chapter  II. 

1.  Pressure  waveforms 

The  pressure  waveform  radiated  by  a  thermoacoustic  source  moving 
at  Mach  1.6  is  shown  in  Fig.  24.  In  both  cases  the  pressures  ere  normalized 
to  one  at  their  peak  values.  The  values  of  the  parameters  relevant  to 
Fig.  24  ere  :  r0  =  4  m  ,  0O  =  80°  ,  $0  =  0°  ,  f0  =  35  kHz  ,  and  tp  =  0.8  ms. 
There  is  fairly  good  agreement  between  the  theory  (Fig.  24(e))  and  the 
experiment  (Fig.  24(b)).  Note  that  time  inversion  does  occur  at  such  a 
supersonic  source  velocity,  8nd  it  is  very  clearly  observed  experimentally. 
The  slight  discrepancy  (93  ys)  between  the  predicted  and  experimental 
values  of  the  arrival  time  of  the  acoustic  pulse  at  the  receiver  is 
attributed  to  the  fsct  that  the  time  delay  between  the  discharge  of  the 
laser  fleshlamp  and  the  actual  lasing  of  the  rod  (about  180  ys)  was  not 
perfectly  constant  so  that  the  initial  coordinates  r0  and  90  of  the 
thermoacoustic  source  were  known  within  only  a  certain  accuracy. 

Next  we  turn  our  attention  to  the  pressure  waveform  when  the 
source  velocity  is  such  that  M  sin0ocos$o  =  1  .  Figure  25(e)  shows  the 
results  of  three  different  theoretical  prediction  schemes.  The  solid  line  in 
Fig.  25(a)  represents  the  prediction  obtained  from  the  numerical  analysis 
presented  in  the  first  section  of  this  chapter.  The  long  dashed  line  in  the 
same  figure  was  obtained  by  using  a  simplified  approach  which  is 
described  in  Appendix  H  .  In  this  approach,  the  pressure  waveform  can  be 
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LASER-INDUCED  MACH  WAVE 
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approximated  by  the  second  time  derivative  of  the  impulse  response  of  a 
cross  section  of  the  laser  beam.  The  third  curve  (short  dashed  line)  in 
Fig.  25(a)  represents  the  analytical  prediction  of  Lyamshev  and  Sedov  given 
by  Eq.  (5)  in  Ref.  26.  The  values  of  the  parameters  relevant  to  Fig.  25  ere 
r0  =  3.50  m  ,  60  =  80°  ,  ^  0°  ,  li  =  1.03  ,  f0  =  20  KHz  ,  and  tp  *  0.8  ms. 

The  three  curves  in  Fig.  25  have  the  same  main  features:  an 
exponential  growth,  followed  by  a  strong  inverted  pulse  and  an  exponential 
decay.  The  exponential  growth  is  characteristic  of  the  arrival  at  the 
receiver  of  wavelets  coming  successively  from  a  depth  ^j=rocos0o  to  a  zero 
depth  at  the  air-water  boundary.  The  inverted  nature  of  the  strong  pulse  is 
characteristic  of  a  pressure  release  boundary  condition  at  the  air-water 
interface.  The  exponential  decay  is  the  result  of  the  exponential  lew  of 
absorption  of  light  in  water.  The  last  wavelets  arriving  at  the  receiver 
come  from  the  furthest  points  on  the  column  of  water  being  illuminated  by 
the  laser  and  therefore  they  have  very  low  intensity. 

Figure  25(b)  shows  an  acoustic  pulse  recorded  experimentally  in  a 
Mach  wave.  There  is  an  overall  agreement  with  the  theoretical  predictions 
given  in  Fig.  25(e)  although  the  result  shown  in  Fig.  25(b)  has  a  higher  peak 
after  the  inverted  pulse  end  shows  more  high  frequency  content.  The  pulse 
was  recorded  by  an  H-23  hydrophone,  amplified  by  40.5  dB  and  filtered 
between  1  and  200  kHz.  The  beam  radius  on  the  surface  of  the  water  was 
about  7  mm,  and  the  corresponding  peak  sound  level  was  measured  to  be 
about  142  dB  re  1  yPa  at  3.50  m  from  the  source.  The  asymmetry  in  the 
experimental  waveform  is  attributed  to  nonlinear  effects.  It  has  been 
shown27  that  nonlinear  effects  associated  with  the  temperature 


dependence  of  the  thermodynamic  coefficients  (p,cp)  cause  an  increase  in 
the  tail  of  the  received  signal.  This  is  due  to  the  fact  that  the  heating  of  the 
medium  being  illuminated  increases  the  coefficient  of  thermal  expansion 
and,  consequently,  the  efficiency  of  the  conversion  of  electromagnetic 
energy  into  acoustic  energy. 

2.  Directivity  patterns 

Next  consider  the  directivity  pattern  of  an  MTS  in  the  x-z  plane.  It  is 
important  to  note  th8t  the  motion  of  the  MTS  induces  8  variation  of  the 
angle  of  observation  between  the  source  and  the  receiver,  and  hence  e 
variation  of  the  Doppler  shifted  frequency,  so  that  the  received  frequency 
is  in  fact  a  function  of  source  location.  The  directivity  pattern  of  8n  MTS  is 
not  therefore  given  8t  a  single  frequency  but  rather  in  a  frequency  interval. 
The  smaller  the  ratio  Mctp/r0  ,  the  smaller  the  frequency  interval.  It  is 
therefore  necessary  to  define  the  concept  of  directivity  for  a  moving 
source.  In  our  study  a  directivity  pattern  is  a  plot  of  the  sound  pressure 
level  8s  a  function  of  the  angle  of  observation  between  the  source  end  the 
receiver,  at  the  origin  (time  t=0)  of  the  laser  pulse,  in  other  words,  all  the 
parameters  such  8S  modulation  frequency, laser  pulse  duration,  and  initial 
distance  of  observation,  are  kept  constant  in  a  directivity  pattern.  The  only 
parameter  allowed  to  change  is  the  initial  angle  of  observation  $0  ,  in  the 
case  of  a  directivity  pattern  in  the  vertical  plane,  and  4>0  in  the  case  of  a 
directivity  pattern  in  the  horizontal  plane.  The  frequency  content  of  the 
received  signal  is  therefore  allowed  to  change  (sometimes  drastically)  in 


the  same  directivity  pattern.  The  sound  pressure  level  is  obtained  by  taking 
20  times  the  logarithm  based  10  of  the  ratio  of  the  peak-to-peek  pressure 
signal  received  at  the  hydrophone  to  a  reference  value  of  1  yPa.  It  should  be 
noted  that  the  shape  of  the  directivity  patterns  does  not  change 
significantly  when  rms  pressures  are  used  instead  of  pe8k-to-peak 
pressures.  A  directivity  pattern  defined  in  such  a  way  is  very  practical 
from  an  experimental  standpoint,  since  it  gives  an  estimate  of  the  sound 
pressure  level  versus  angle  under  the  same  operating  conditions  of  the 
laser 

The  vertical  directivity  of  thermoacoustic  source  is  shown  in  Fig.  26. 
The  theoretical  curves  were  obtained  by  running  the  numerical  program 
described  previously,  over  several  values  of  the  angle  fl0.  The  other 
parameters  were  r0=  3  m  ,  ^  «*  0°  ,  CC  s  13.7  Np/m,  and  tp  =  0.8  ms.  The 
circles  represent  experimental  measurements.  The  left-h8nd  directivity 
pattern  in  Fig.  26  wes  obtained  with  a  stationary  source  modulated  at  35 
kHz  while  the  right-hand  one  wes  obtained  for  a  source  moving  at  Mach  1.6 
and  a  modulation  frequency  of  25  kHz.  The  experiment  confirms  that  even 
though  the  modulation  frequency  wes  lower  for  the  moving  source,  its 
directivity  turns  out  to  be  sharper  than  that  of  the  stationary  source, 
because  the  received  Doppler  shifted  frequency  is  greater  th8n  f0  by  the 
Doppler  factor.  The  half-power  beamwidth  W8S  measured  to  be  11°  for  the 
stationary  source  and  7°  for  the  moving  source.  This  agrees  well  with  the 
predicted  respective  vslues  of  11°  snd  6°.  Although  the  source  appears  at 
first  glance  to  be  a  broadside  array  whose  main  lobe  should  be  at  90  »  90°, 
the  presence  of  the  air-water  interface  produces  a  tilt  in  the  lobe  so  that 
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there  is  no  pressure  at  the  boundary.  The  tilt  angle  for  the  main  lobe  W8S 
measured  to  be  84°  for  the  stationary  source  and  86°  for  the  supersonic 
source.  Again  this  confirms  the  validity  of  the  model  which  predicts  83°  in 
the  first  case  8nd  85°  in  the  second. 

Two  experimental  directivity  patterns  in  the  vertical  x-z  plane  of  8n 
MTS  moving  at  Mach  one  are  shown  in  Fig.  27.  The  laser  intensity  was 
modulated  at  25  kHz  and  35  kHz  respectively  for  the  two  directivity 
patterns,  and  in  both  cases  the  initial  coordinates  of  the  MTS  were  r0  =  4  m 
and  <t>0  =  0°.  The  two  curves  are  very  similar  which  confirms  the  fact  that 
the  modulation  frequency  does  not  play  8n  important  role  in  the  directivity 
pattern  of  an  MTS  moving  at  transonic  velocity.  Theoretical  curves  8re  not 
shown  in  Fig.  27  because  the  numerical  codes  seem  to  yield  unstable  results 
near  M8ch  one.  The  reason  for  this  discrepancy  has  not  been  determined  yet. 
It  should  be  emphasized  th8t  experimental  results  for  a  source  moving  at 
Mach  1  8re  highly  dependent  upon  the  frequency  response  of  the  receiving 
hardware  (hydrophone-emplifier-filter)  which  should  optimally  have  a  very 
large  bandwith.  However,  the  H-56  hydrophone  used  in  this  experiment  has 
an  upper  frequency  limit  of  about  70  kHz  (see  calibration  curves  in 
Appendix  B).  Moreover,  the  correction  for  the  vertical  sensitivity  of  the 
hydrophone  is  very  difficult  to  assess  when  the  received  signal  is 
broadband.  Thus  quantitative  measurements  around  Mach  1  are  to  be  taken 
cautiously. 

The  directivity  pattern  in  a  plane  parallel  to  the  interface  between 
air  and  water  (x-y  plane)  will  be  referred  to  as  a  horizontal  directivity 
pattern.  Figure  28  shows  the  horizontal  directivity  pattern  of  an  MTS  for 
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three  different  source  velocities,  but  under  similar  conditions:  r0  *  2  m  , 
f0  =  25kHz  and  B0  *  80°.  When  the  source  is  stationary,  the  measured 
acoustic  radiation  is  fairly  omnidirectional,  as  predicted  by  the  theory 
given  by  the  solid  line  in  Fig.  28(a).  The  results  are  normalized  by  the 
pressure  signal  received  at  40=90°.  The  theoretical  curves  were  obtained 
by  running  the  numerical  program  described  previously,  over  several  values 
of  the  angle  ♦<,.  Fig.  28(b),  whose  scale  differs  from  that  of  Fig.  28(a), 
shows  the  horizontal  directivity  pattern  of  an  MTS  moving  at  M  =  0.4 
towards  the  receiver,  and  Fig.  28(c)  applies  to  the  case  of  an  MTS  moving  at 
M  *  0.5  away  from  the  receiver.  In  both  cases  there  is  reasonable 
agreement  between  theory  and  experiment  and  the  results  show  that  indeed 
acoustic  radiation  is  stronger  in  the  direction  of  motion  of  the  source.  The 
theoretical  model  seems,  however,  to  overestimate  the  acoustic  level  in 
the  direction  opposite  to  the  motion  of  the  source  and  this  discrepancy  is 
not  fully  understood. 

The  horizontal  directivity  pattern  of  8  transonic  TS  was  investigated 
experimentally  end  the  result  is  given  in  Fig.  29.  It  shows  that  more  than  30 
dB  in  acoustic  amplitude  can  be  gained  in  the  Mach  cone  where  a  very  strong 
broadband  signal  is  emitted  by  the  moving  source,  as  indicated  by  the  upper 
time  waveform  in  Fig.  29.  Both  time  waveforms  in  Fig.  29  are  normalized  to 
unity.  When  the  Doppler  shifted  wavelength  X<j*Xol1-Msin0ocos*ol  is  of  the 
order  of  the  laser  beam  diameter,  diffraction  effects  are  obviously 
important  and  a  strong  diffraction  loss  is  expected  to  reduce  the  sound 
level.  This  happens  when 
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that  is  to  say,  when  32°  ,  and  this  may  be  an  explanation  for  the 
sudden  drop  in  level  around  $0  =  30°.  The  second  waveform  shown  in 
Fig.  29  indicates  that  for  angles  *0  >  ♦jjm  (i.e.,  well  outside  the  Mach 
cone),  the  received  pressure  waveform  is  not  a  strong  transient  but  rather 
a  waveform  modulated  at  its  Doppler  shifted  frequency.  Hence  there  are 
two  angular  regions  associated  with  transonic  source  motion.  In  one  region 
the  waveform  is  a  Mach  wave  and  in  the  other  region  the  laser  modulation 
waveform  is  dominant. 

Directivity  patterns  were  also  obtained  for  the  case  of  a  long  and 
narrow  thermoacoustic  source  created  by  illuminating  the  water  with  the 
ruby  laser  (a>1.5  Np/m).  Experimental  date  are  given  in  Appendix  I. 

3.  Sound  level  versus  Msch  number 

In  another  set  of  measurements  we  investigated  the  dependence  of 
the  sound  pressure  level  on  the  source  velocity.  The  results  are  shown  in 
Fig.  30,  where  the  dashed  line  represents  the  theoretical  predictions 
obtained  by  numerical  analysis  of  the  pressure  waveforms  computed  for 
several  values  of  the  Mach  number  as  seen  by  the  receiver,  M  *  M  sine0 
cos40.  The  general  shape  of  the  predicted  curve  follows  the  trend  of  the 
experimental  data  with,  however,  an  unexplained  discrepancy  in  the 
amplitude  of  the  sharp  peak  which  occurs  around  M=1.  It  is  possible  that  in 
the  transonic  regime,  some  nonlinear  effects,  which  have  not  been  included 
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in  the  present  model,  become  of  significant  importance.  The  results  shown 
in  Fig.  30  are  normalized  so  that  both  theory  and  experiment  have  a  value  of 
0  dB  at  M=0  In  fact  the  absolute  level  measured  experimentally  at  M=0  was 
of  116.7  dB,  whereas  the  theory  predicts  a  level  of  109.4  dB.  This  significant 
difference  is  so  far  unexplained. 

The  physical  interpretation  of  the  two  symmetric  dips  around  M  =  1  is 
that  diffraction  effects  due  to  the  finite  width  of  the  laser  beam  become 
important  and  tend  to  lower  significantly  the  acoustic  radiation.  In  the  case 
of  monochromatic  radiation  we  have  seen  that  the  frequency  filtering 
associated  with  finite  beamwidth  effects  is  given  by  Eq.  (4.6)  which  shows 
that  a  diffraction  loss  of  3  dB  is  expected  to  occur  when 

(  kd  8  sin0o/2  )2  =  Ln  (  10  3/20  )  ,  (4.8) 

that  is  to  say  ,  when  M  =  1  ±  5.35  0  ,  where  Q  ■  f0  /  (  c/a  sin0o)  is  a 
nondimensional  frequency.  In  our  experiments  Q  is  0.0828  so  that  the  3  dB 
points  8re  expected  to  occur  around  M  =  0.6  and  M  *  1.4 .  As  indicated  by 
Fig.  30,  this  is  confirmed  experimentally.  When  M  approaches  unity,  the 
acoustic  signal  becomes  so  short  that  the  monochromatic  assumption  is  not 
valid  any  more  and  thus  Eq.  (4.6)  cannot  be  used  to  predict  the  diffraction 
loss.  In  fact  at  M  =  1,  all  of  the  acoustic  disturbances  emitted  during  the 
motion  of  the  source  seem  to  add  coherently  ,  so  that  a  strong  transient  is 
observed  at  the  receiver.  A  gain  of  more  than  25  dB  has  been  observed  over 
the  amplitude  of  a  stationary  laser  source  of  similar  characteristics. 

The  three  normalized  waveforms  shown  in  the  windows  of  Fig.  30 
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indicate  the  typical  behavior  of  the  pressure  waveform  for  the  three 
regions  of  the  curve.  In  the  first  region  the  thermoecoustic  signal  is 
modulated  at  its  Doppler  frequency.  When  M  is  further  increased, 
diffraction  effects  reduce  the  acoustic  level  until  the  monochromatic 
assumption  does  not  hold  and  a  very  strong  transient  is  generated  as 
indicated  in  the  second  region.  In  the  third  region,  when  diffraction  effects 
become  negligible,  one  C8n  observe  time-inverted  thermoecoustic  signals. 

4.  Spreading  curves 

The  spreading  curves  shown  in  Fig.  31  show  the  sound  level 
dependence  on  the  distance  r0  between  the  thermoecoustic  source  and  the 
receiver.  These  curves  are  plotted  for  three  different  source  velocities. 
M=0  ,  M=1  ,  and  M=1.5.  The  other  parameters  of  interest  are  f0=35  kHz  , 
0O=84°  ,  and  40=0°  .  The  predicted  spreading  curves,  indicated  by  a  solid 
line  in  Fig.  31 ,  were  obtained  as  follows. 

In  the  case  of  a  TS,  the  predicted  spreading  curve  was  obtained  from 
Eq.  (3.32).  For  the  case  of  an  MTS,  Eq.  (3.32)  does  not  apply  because  the 
received  frequency  is  actually  time  dependent.  In  Fig.  31,  the  theoretical 
predictions  for  transonic  and  supersonic  regimes  are  based  on  the 
assumption  that  the  Doppler  shifted  wavelength  characteristic  of  the 
acoustic  radiation  is  small  compered  to  the  distance  of  observation,  so 
that  spherical  spreading  (6  dB  per  doubling  distance)  dominates  over 
cylindrical  spreading  (3  dB  per  doubling  distance).  Experimental  data 
support  this  assumption  very  well,  as  indicated  by  Fig.  31. 
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C.  Conclusions 


The  acoustic  characteristics  of  a  moving  thermoacoustic  source  created 
by  moving  a  modulated  laser  across  a  fluid-air  boundary  have  been  examined 
by  imagining  the  source  to  consist  of  many  line  sources  each  of  which  are 
excited  at  successive  times  corresponding  to  the  motion  of  the  laser  source 
along  its  path.  In  this  manner  it  is  possible  to  construct  the  acoustic 
pressure  received  by  a  hydrophone  by  adding  the  sound  contribution  of  each 
line  source  with  the  appropriate  delay.  The  result  is  a  pseudo-convolution 
between  the  impulse  response  of  the  system  and  the  acoustic  source 
strength.  This  analysis  lends  itself  well  to  numerical  simulation.  It  is  valid 
in  the  nearfield  of  the  thermoacoustic  source  and,  in  principle,  for  any 
source  velocities.  In  the  case  of  a  stationary  source,  the  impulse  response 
of  the  system  is  fixed  in  time,  so  the  pressure  calculation  reduces  to  a  true 
convolution.  Further,  it  is  possible  to  include  the  effect  of  the  finite 
beamwidth  of  the  laser  beam  by  convolving  the  pressure  calculated 
assuming  an  infinitesimal  beam  with  the  impulse  response  of  the  system  to 
laser  heat  applied  across  the  beam  at  a  fixed  depth  in  the  fluid.  Although 
this  impulse  response  is  a  function  of  the  location  and  velocity  of  the 
source,  the  corrections  produced  by  the  nonstationarity  of  the  source  are 
small  and  may  be  ignored. 

Experimental  results  have  been  presented  which  indicate  excellent 
agreement  between  the  measured  and  predicted  waveforms  including  the 
subsonic,  transonic,  and  supersonic  cases.  It  has  been  shown  that  the 
waveforms  in  these  three  velocity  regions  differ  markedly  from  one  another. 


In  the  subsonic  region  the  waveform  consists  of  a  modulated  signal  with  an 
envelope  that  first  grows  rapidly  and  then  decays  slowly.  As  the  source 
velocity  increases,  the  amplitude  and  frequency  of  the  waveform  increase 
and  the  total  time  of  duration  decreases.  In  the  supersonic  region,  the 
pressure  waveform  is  very  similar  to  that  for  the  subsonic  region,  but  the 
signal  is  time  inverted.  This  time  inversion  phenomenon  is  caused  by 
moving  the  source  toward  the  receiver  at  a  speed  faster  than  the  speed  of 
sound.  In  the  transonic  region  the  waveform  shows  almost  no  frequency 
modulation,  but  is  a  very  fast  negative  pulse  with  smaller  leading  and 
lagging  positive  exponential  tails.  Previous  analyses18'26  in  which  the 
pressure  field  radiated  by  an  MTS  is  replaced  by  the  pressure  field  of  an 
equivalent  stationary  array,  corrected  by  a  Doppler  shift,  have  failed  at 
transonic  source  velocities  because  of  neglecting  the  fact  that  the  Mach 
number  of  the  source  as  seen  by  the  receiver  depends  on  angle  and  hence 
location  of  the  source. 

Experimental  measurements  and  theoretical  predictions  for  the 
directivity  patterns  produced  by  a  moving  laser  source  were  shown  to  be  in 
good  agreement.  These  results  were  presented  for  both  the  horizontal  and 
vertical  planes.  The  results  also  show  that  near  transonic  source  velocity 
the  modulation  frequency  of  the  laser  is  not  a  significant  determinant  of  the 
directivity  pattern.  Below  and  above  transonic  source  velocity,  the 
modulation  frequency  has  a  direct  effect  on  the  directivity  pattern;  a 
narrower  beam  is  produced  as  the  modulation  frequency  increases. 

Experimental  results  and  theoretical  predictions  were  also  presented 
for  the  sound  pressure  level  as  a  function  of  the  source  velocity  as  seen  by  a 


hydrophone  at  a  fixed  position.  These  results  show  the  largest  discrepancy 
of  all  of  the  theory  and  experiment  comparisons.  Although  the  curves  have  a 
similar  shape  there  is  a  notable  and  large  disagreement  near  Mach  one. 
Experimentally,  large  acoustic  signals  were  observed  near  transonic  source 
velocity,  with  amplification  of  almost  20  dB  above  the  levels  in  the 
subsonic  and  supersonic  region,  and  more  than  30  dB  in  the  horizontal 
directivity  pattern  in  the  direction  of  motion  of  the  source.  Theoretically, 
the  predicted  amplification  at  transonic  velocity  is  much  smaller  than  the 
one  measured  experimentally.  This  suggests  that  there  may  be  nonlinear 
mechanisms  which  become  important  near  Mach  one  and  which  have  not  been 
included  in  the  physical  model.  More  work  is  necessary  to  explain  the  exact 
nature  of  this  discrepancy. 

Finally,  experimental  results  for  the  spreading  of  the  moving  and 
stationary  laser  source  were  compared  to  theoretical  predictions.  It  was 
shown  that  in  the  nearfield,  the  laser  source  exhibits  cylindrical  spreading, 
then  switches  over  to  spherical  spreading  in  the  farfield.  These  results 
agree  well  with  physical  intuition. 
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CHAPTER  V 


THE  DOPPLER  SHIFT  OF  AN  ACOUSTIC  SOURCE  MOVING 
AT  TRANSONIC  VELOCITY24 

In  order  to  improve  the  intrinsically  low  efficiency  of  the  conversion 
by  thermal  expansion,  of  light  into  sound,  Soviet  ph  sicists  have 
studied28-32  closely  the  sound  field  radiated  by  a  thermoacoustic 
source  moving  at  transonic  velocity.  They  have  shown  experimentally  that 
indeed  a  substantial  gain  in  sound  level  occurs  when  the  source  is  moving 
trensonicelly.  This  has  also  been  confirmed  experimentally  in  the  previous 
chapter.  It  is  a  striking  fact,  however,  that  most  theories  which  rely  on  the 
Doppler  shift  to  account  for  the  motion  of  the  source  break  down  at 
transonic  source  velocity.  The  reason  for  this  is  that,  according  to 
classical  theory,  the  Doppler  shift  associated  with  an  acoustic  source 
moving  at  constant  velocity  induces  a  straining  of  the  time  coordinate,  8s 
perceived  at  the  receiver.  This  time  straining  is  described  by 

f  =  t  ( 1  -  M  cos  0  )  ,  (5^ 

where  t'  denotes  the  time  coordinate  in  the  fixed  frame  of  reference 
(receiver),  t  is  the  time  coordinate  in  the  moving  frame  of  reference 
(source),  M  is  the  Mach  number  of  the  source,  and  O  is  the  angle  between  the 
direction  of  motion  of  the  source  and  the  direction  of  observation.  At 


transonic  velocities,  t*  can  assume  a  value  of  zero,  in  which  case  all  the 
acoustic  energy  emitted  by  the  source  is  received  in  an  infinitesimal 
increment  of  time.  This  mplies  an  infinite  instantaneous  power  which  is 
not  physically  realizable.  The  object  of  this  chapter  is  to  use  the  method  of 
characteristics  to  find  the  exact  Doppler  shift  of  a  transonic  source.  Two 
cases  are  considered:  a  point  source  and  a  line  source  of  finite  length,  both 
moving  rectilineerly  towards  or  away  from  a  point  receiver. 

A.  Point  source 

1.  Characteristics 

Consider  a  point  source  moving  at  constant  velocity  in  a  straight 
line,  as  shown  in  Fig.  32.  Let  us  assume  that  the  source  is  emitting  a 
rectangular  pulse  of  duration  tp.  This  simplification  is  justified  by  the  fact 
that  we  are  primarily  concerned  with  time  scales  and  not  wave  shapes.  The 
present  analysis  is  not,  however,  restricted  to  transient  sources.  It  is 
understood  th8t  tp  is  just  a  characteristic  time  of  the  acoustic  emission 
The  Doppler  factor  D  may  then  be  described  as  the  ratio  of  the  pulse 
duration  tp  recorded  at  the  receiver  to  the  pulse  duration  tp  emitted  by  the 
source. 

D  =  tp  /  tp  (5.2) 

Note  that  tD,  and  therefore  D,  are  allowed  to  take  on  negative  values  which 


is  whet  will  happen  et  supersonic  source  velocities.  It  is  en  effect  of  the 
time  reversal  phenemenon  which  occurs  when  a  source  is  moving  faster 
than  the  disturbances  it  generates. 

Another  quantity  relevant  to  the  motion  of  en  acoustic  source  is  the 
Doppler  gain  6,  defined  by 


G  =  1 1  /  D I 


(5.3) 


The  gain  can  be  explained  in  terms  of  conservation  of  energy.  The  motion  of 
the  source  induces  a  change  in  the  time  scales,  but  the  net  energy  flux  et 
the  receiver  is  unchanged.  The  average  power  in  the  received  signal  must 
change  accordingly  then.  Reducing  the  received  pulse  duration  therefore 
increases  the  average  received  power  by  the  factor  G.  It  is  important  to 
note  that,  if  the  source  is  directional,  the  motion  of  the  source  induces  e 
redistribution  of  energy  as  a  function  of  angle  of  observation,  and  thus  the 
gain  G  is  not  just  the  gain  in  amplitude  of  the  acoustic  pressure.  This 
dependence  of  the  source  directivity  properties  as  seen  by  the  receiver,  on 
source  motion,  makes  the  problem  much  more  complicated  and  it  is  beyond 
the  scope  of  the  analysis  presented  in  this  chapter.  Time  inversion  has  no 
effect  on  the  gain  G,  since  only  its  amplitude  has  a  significant  meaning.  This 
is  why  an  absolute  value  was  introduced  in  Eq.  (5.3). 

The  problem  is  therefore  to  find  the  pulse  duration  tp  recorded  at 
the  receiver.  A  convenient  approach  for  solving  this  problem  is  to  use  the 
method  of  characteristics  which  consists  of  finding  the  time  of  arrival  of 
each  wavelet  coming  from  the  source  during  its  motion  along  the  x-exis. 


Let  us  define  tp(x)  es  the  time  of  arrival  at  the  receiver  of  a  wavelet 
emitted  when  the  source  was  at  a  distance  x  from  the  origin.  The  zero  of 
time  is  chosen  for  convenience  to  coincide  with  the  zero  of  space,  i.e. ,  t=0 
when  the  source  is  at  x»0.  Then  tR(x)  can  be  expressed  as  follows: 

tR(x)  -  (x/v)  ♦  (rx/c)  ,  (5.4) 

where  c  is  the  sound  speed  in  the  medium  and  rx  is  the  distance  between  the 
source  and  the  receiver  when  the  source  is  at  a  distance  x  from  the  origin. 
It  con  easily  be  shown  that 

tp(x;M)  =  (x/Mc)  +  c~M r02-2r0x  cos0o+x2),/2  ,  (5.5) 

where  the  subscript  o  refers  to  sn  initial  value  (x-O;  t*0),  and  M=v/c  is  the 
Mach  number  of  the  source.  Note  that  cos0o  corresponds  to  sin#0  cos*0  in 
the  notation  of  Chapter  IV.  The  slope  of  the  characteristics  tR(x;M)  in  the 
x-t  plane  may  be  found  by  differentiating  Eq.  (5.5)  with  respect  to  x.  It  is 
found  that 

d/dx  { tR(x)  J  =  ( 1  -  M  cos0x  )/v  ,  (5.6) 

where  0X  indicates  that  the  angle  of  observation  0  depends  on  the  source 
location  x.  If  M  is  a  constant  along  the  source  path,  fl  cos0x  is  8  nonlinear 
function  of  the  spatial  coordinate  x.  This  means  that,  in  general,  the 


characteristics  are  curved. 

Subsonic,  transonic,  and  supersonic  regimes  may  be  defined  from  the 
slopes  of  the  characteristics.  Equation  (5.6)  shows  that  a  negative  slope  can 
occur  only  if  M  >  l/cos0x  ,  i.e.,  only  at  supersonic  velocities.  If  the  source 
is  moving  at  supersonic  velocity,  but  has  already  passed  over  the  receiver 
and  is  moving  away  from  it,  cos0x  is  negative  and  therefore  the  slope  is 
positive.  The  transonic  regime  occurs  only  for  11  «  l/cos0x  ,  i.e.,  when  the 
slope  is  zero.  In  the  Soviet  literature,  cos0*l/M  is  sometimes  referred  to 
as  the  Cerenkov  direction.  The  subsonic  regime  occurs  for  positive  slopes 
of  the  characteristics.  The  characteristics  tp(x;M)  are  plotted  in  Fig.  33  for 
various  values  of  the  Mach  number  M. 

In  the  limit  of  x  »  r0 , 0X  approaches  if ,  and  from  Eq.  (5.6),  the  slope 
of  the  characteristics  goes  to  a  maximum  value  of  (1/v  ♦  1/c).  The 
interpretation  of  this  result  is  that  both  subsonically  and  supersonically 
moving  point  sources  do  not  induce  time-inversion  on  the  receiver 
waveform,  once  they  are  moving  8wey  from  the  receiver.  Regardless  of 
the  source  speed,  supersonic  or  subsonic,  wavelets  emitted  at  a  given 
time  will  arrive  at  the  receiver  before  those  emitted  a  short  time  later. 
For  supersonically  moving  sources,  the  time  gap  between  the  two  wavelets 
is  accentuated  because  the  later  wavelet  has  so  much  farther  to  travel. 

It  therefore  seems  appropriate  to  make  a  distinction  between  a 
supersonic  source  and  a  source  that  is  seen  as  supersonic  by  the  receiver. 
In  the  first  case,  the  Mach  number  of  the  source  is  such  that  M  cos0o  >  1  , 
whereas  in  the  second  case,  time  inversion  occurs  and  the  slope  of  the 
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characteristic  is  negative.  This  distinction  is  emphasized  in  Fig.  33  by 
indicating  the  region  of  time  inversion. 

Also  of  interest  is  the  locus  of  the  points  of  zero  slope  in  the 
characteristics,  i.e.,  the  transonic  region.  From  Eq.  (5.5),  we  can  set  d/dx 
(tp(x)l  =  0  ,  and  obtain  the  following  equation  for  the  points  of  zero  slope: 

xt  -  r0  cos0q  [  1  -  V  (cos~20o  -  1  )/(M2-1)  ]  (5.7) 

The  transonic  region  is  indicated  by  the  dashed  line  on  Fig.  33. 

Another  interesting  parameter  is  the  distance  Xq  at  which  the 
characteristics  cross  the  line  tR  *  r0/c  .  Wavelets  emitted  from  this 
location  will  arrive  at  the  receiver  simultaneously  with  those  emitted  at 
x  ■  0,  t  *  0.  Clearly,  Xc  exists  only  for  a  source  moving  at  supersonic  speed, 
and  Xc  is  always  greater  then  x*  it  is  possible  to  solve  for  xc  by  setting 
tR(xc)  *  r0/c  in  Eq.  (5.5).  It  is  found  that,  provided  ri  >  1/cos0o  , 

xc«2r0f1  { ( M  cos  0O  -1 )/(  (5.8) 

2.  Pulse  duration 

The  x-t  plane  is  a  convenient  way  to  find  the  received  pulse  duration  and 
therefore  the  corresponding  Doppler  shift.  First  let  us  define  xN  8S  the 
distance  traveled  by  the  source  during  the  pulse  duration,  x^  ■  v  tn  .  Since 


tg(x)  is  the  time  of  arrive!  of  a  wavelet  coming  from  x,  the  set  of  all  tp(x) 
for  0  <  x  <  xn  ,  written  as  t^Oixix^),  subsequently  represents  the  time  gap 
between  the  arrivals  of  all  wavelets  emitted  during  the  motion  of  the 
source.  The  recorded  pulse  duration  tp  is  the  difference  between  the 
maximum  and  the  minimum  values  of  the  set.  In  other  words,  tp  is  the 
projection  on  the  time  axis  of  the  characteristic  for  8  given  Mach  number, 
over  all  x,  such  that  o  <  x  <  x^  . 

tp  ■  Pro3t  t  tR(  0  i  x  <  xN  )  1  (5.9) 

From  Eqs.  (5.2)  and  (5.3),  it  is  clear  that  in  order  for  the  Doppler  gain 
to  be  infinite,  it  is  necessary  for  tp  to  be  zero.  This  requires  ell  of  the 
wavelets  emitted  by  the  source  to  be  received  simultaneously,  i.e.,  that  the 
characteristic  be  a  straight  line  at  tp  =  r0/c  .  From  Eq.  (5.6),  we  see  th8t 
the  characteristic  cannot  be  a  horizontal  line  unless  M  cos0x  is  equal  to 
unity.  For  8  source  moving  at  uniform  Mach  number,  this  requires  Ox  to  be  a 
constant,  and  the  only  geometry  which  permits  0X  to  be  constant  places  the 
receiver  directly  on  the  path  of  the  source.  If  the  source  is  not  constrained 
to  move  at  a  constant  Mach  number,  other  geometries  which  permit  tp  to  be 
zero  are  possible,  but  this  would  violate  the  assumption  of  a  source  moving 
at  constant  velocity,  and  the  preceding  analysis  would  be  invalid.  Hence, 
unless  the  receiver  is  directly  on  the  path  line  of  the  source,  cp  end  the 
Doppler  gain  are  always  finite  If  the  receiver  is  on  the  source  path  line,  it 
is  necessary  to  include  diffraction  effects  induced  by  the  finite  dimension 
of  any  real  physical  source,  in  order  to  prove  that,  in  reality,  the  received 


pulse  duration  cannot  be  less  than  the  transit  time  across  the  source  and, 
hence,  that  the  Doppler  gain  is  still  bounded. 

If  we  restrict  the  analysis  to  the  case  of  a  source  moving  at  constant 
velocity,  it  is  straightforward  to  find  an  analytical  expression  for  the 
received  pulse  duration  end  the  corresponding  Doppler  factor. 

First,  consider  subsonic  motion  of  the  point  source.  As  mentioned 
earlier,  when  li  <  l/cos0o  ,  the  slope  of  the  characteristic  is  always 
positive.  This  means  that  the  order  of  arrival  of  the  wavelets  at  the 
receiver  remains  unchanged  as  compared  to  their  order  of  emission.  In  this 
case,  the  "projection  operator"  defined  by  Eq.  (5.9)  is  simply 


where 


tp  -  tR(xN)  -  tR(0) 


xN  =  M  c  tp 


(5.10) 


(5.11) 


Combining  with  Eq.  (5.5)  yields 


tp  *  (xN/v)  ♦  c"1  y  r02-2  roxNcos0o+xN2  -  (r0/c) 


(5.12) 


Substituting  Eq.  (5.12)  into  Eq.  (5.2)  yields  an  expression  for  the  Doppler 


D  3  1  ♦  (M/xn)V  r02  -  2  r0xN  cos0o  ♦  xN2  -  (r0/ctp) 


(5.13) 


Defining  a  nondimensionel  distance  Rg  =  r0/ctp  (this  quantity  can  also  be 


regarded  as  8n  inverse  nondimensional  pulse  duration,  or  a  nondimensional 
frequency),  we  finally  rewrite  the  Doppler  factor  in  the  form 


D  -  1  -  R0  l  1  -  V  1  -  2(M/R0)  cos0q  +  (M/R0)2  ]  (5.14) 

Equation  (5.14)  gives  a  practical  estimate  of  the  Doppler  factor  of  a  moving 
source,  when  one  wants  to  take  into  account  the  change  in  angle  of 
observation  as  the  source  is  moving.  As  an  example,  the  received  pulse 
duration  was  measured  in  the  study  of  moving  thermoacoustic  sources,  and 
it  was  compered  with  the  predicted  value  of  D  tp  (see  Eq.  (5.2)).  The 
Doppler  shift  D  W8S  computed  first  assuming  that  the  angle  of  observation 
0  was  constant  during  the  motion  of  the  source,  and  second,  in  the  more 
accurate  way  described  by  Eq.  (5.14).  Figure  34  shows  that  Eq.  (5.14)  fits 
the  data  better  than  the  simple  estimate  D  ■  1  -  (1  cos  0O  .  The  experimental 
data  were  obtained  from  the  study  of  the  horizontal  directivity  described 
in  Section  IV.B.  The  Msch  number  of  the  source  was  0.4  ,  the  radius  of 
observation  was  2  m,  the  depth  of  the  receiver  was  0.35  m  (0o=8O°),  and 
the  modulation  frequency  was  25  kHz.  This  simple  experiment  shows  that 
Eq.  (5.14)  is  useful  in  some  practical  cases. 

Next,  consider  supersonic  motion  of  the  source.  In  this  case,  the 
problem  is  more  complicated  since  the  slope  of  the  characteristics  is 
negative  until  x  =  x^,  then  it  becomes  positive  (see  Fig.  33).  Therefore  the 
projection  operation  defined  by  Eq.  (5.9)  can  be  expressed  analytically  as 
follows. 
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If  XNiXt  , 

cp  =  tR^N*  “ 

(5.15a) 

If  X^  i  X  i  xc  , 

tp  =  tR(xt)  -  tR(0) 

(5.15b) 

If  xN  >  xc  , 

cp  =  lR<xN*  "  ^*0 

(5.15c) 

where  x^,  x^,  and  xc  are  defined  by  Eqs.  (5.11),  (5.7),  and  (5.8).  The  Doppler 
shift  can  be  expressed  explicitly  by  making  use  of  Eqs.  (5.5)  end  (5.7).  The 
results  are  summarized  below  in  nondimension8l  form. 


For  xN  <  x^ , 

D0  =  1  -  R0  [1  -  V7  -  2(M/R0)  COS  0q  ♦  (M/R0)2  1 
For  xt  <  x^  s  xc 

D]  =  R0  1  (cos0o/M)  -  (sinGo/MVT?7! )  -  1  ♦  sinO^ \  1 
For  xn  >  xc 

D2  s  Do  '  D1 


(5.16a) 


(5.16b) 

(5  16c) 


Notice  that  Eq.  (5.14)  and  Eq.(5.16e)  ere  identical.  By  examining  the  farfield 
limit  of  the  solutions  given  in  Eqs.  (5.14)  and  (5.16),  it  is  possible  to 
compere  the  results  above  to  the  classical  solution.  It  is  necessary  to  take 
the  farfield  limit  because  there  is  an  implicit  assumption  in  the  classical 
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solution  that  the  angle  of  observation  O  is  independent  of  source  location. 
This  assumption  is  only  valid  if  r0  »  c  tp,  i.e.,  in  the  ferfield.  in  the 
ferfield  limit,  Eq.  (5.16a)  applies  and  hence  Eq.  (5.14)  describes  the  Doppler 
shift  for  any  Mach  number  and  any  source  location. 

Expanding  the  square  root  in  Eq.  (5.14)  in  terms  of  the  small 
parameter  (M/R0)  and  retaining  only  first  order  terms,  yields 

D  =  1  -  M  cos  0O  ,  (5.17) 

which  is  the  standard  form  for  the  Doppler  factor.  Note  that,  for  M  = 
1/cos0o  ,  the  gain  becomes  infinite.  With  the  exception  of  the  previously 
discussed  case  of  a  source  moving  at  Mach  1  on  a  straight  line  towards  the 
receiver,  the  infinite  gain  disappears  when  higher  order  terms  are  retained 
in  Eq.  (5.14).  Hence,  in  the  ferfield,  the  analysis  presented  here  recovers  the 
standard  form  but,  with  only  one  exception,  the  analysis  guarantees  finite 
Doppler  gain. 

It  is  also  interesting  to  note  thst,  at  least  for  subsonic  velocities,^ 
the  received  pulse  duration  can  be  expressed  in  the  following  manner 

*N 

tp  =  \  (dtR/dx)dx  (5.18) 

0 

Substituting  Eq.  (5.6)  yields 

*N 

D=  1  -  (xm)'1  [ 


M  cos  0-  dx 


(5.19) 


which  clearly  shows  that  the  global  Doppler  shift  (over  the  whole  path  of 
the  source)  is  an  average  of  the  Doppler  shifts  for  each  position  of  the 
source  along  its  path.  The  integral  in  Eq,  (5.19)  can  be  evaluated  by 
expressing  cosOx  in  terms  of  x.  It  can  easily  be  shown  from  the  geometry 


r„  cos  0„  -  x 


r0*  -  2  x  rQ  COS0Q  ♦  x* 


Note  that  when  the  source  has  passed  over  the  receiver,  cos0o  <  0  so  that 
Eq.  (5.20)  still  applies.  Combining  Eqs.  (5.19)  and  (5.20)  leads  as  expected  to 
Eq.  (5.14),  which  in  turn  reduces  to  the  standard  form  of  the  Doppler  factor, 
Eq.  (5.17),  provided  that  r0  »  ctp  . 

Equations  (5.2),  (5.14),  end  (5.15)  have  been  used  together  with 
Eqs.  (5.5),  (5.7),  (5.8),  end  (5.11)  to  compute  the  Doppler  factor  end  therefore 
the  gain  of  a  moving  point  source.  Figure  35  shows  a  three-dimensional  plot 
of  the  gain  6  versus  Mach  number  end  nondimensional  pulse  duration  T0  * 
ctp/r0  =  (R0)_1.  The  initial  angle  0O  was  set  to  45°.  It  can  be  seen  from 
Fig.  35  that  for  any  Mach  number,  the  acoustic  response  observed  at  the 
receiver  is  always  finite.  Also  of  interest  is  the  fact  that  the  maximum 
gain  in  Fig.  35  occurs  at  a  Mach  number  which  depends  on  the 
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nondimensional  parameter  R0,  whereas  one  would  have  expected  from  the 
classical  theory  an  infinite  response  occuring  always  8t  M  ■  1/cos90.  The 
object  of  the  next  paragraph  is  to  find,  for  a  given  value  of  R0,  the  Mach 
number  which  will  minimize  the  Doppler  factor,  or  maximize  the  Doppler 
gain. 

3.  Maximum  gain  end  optimum  Mach  number 

It  has  been  shown  previously  that  the  gain  of  a  point  source  moving  at 
any  constant  velocity  is  always  bounded,  except  for  9=0.  This  result  leads 
to  questions  about  whether  it  is  possible  to  predict  the  value  of  the 
maximum  gain  and  the  value  of  Mach  number  at  which  the  peek  occurs.  These 
answers  are  easily  found  numerically,  but  in  this  Haragraph,  we  instead 
develop  an  analytical  scheme  for  prediction  purposes. 

First,  consider  the  case  of  a  source  moving  8t  a  subsonic  velocity. 
The  maximum  gain  6*  may  be  found  by  differentiating  Eq.  (5.14)  with 
respect  to  M,  to  find  a  value  of  the  M8Ch  number  M*,  which  minimizes  the 
Doppler  factor  to  a  value  D*.  From  Eq.  (5.14)  it  is  shown  that 


M*  =  R0  cos  90 


(5.21) 


Substituting  M*  into  Eq.  (5.14)  yields 


D*  =  1  -  Rn  ( 1  -  sin  9n  ) 


(5.22) 


C 1  -  R0  ( 1  -  sin  CC0  )  r1 


(5.23) 


Note  that  the  assumption  M  <  l/cos80  can  be  combined  with  Eq.  (5.23)  to 
show  that  this  estimation  scheme  is  valid  only  if 

R0  <  t/cos2e0  (5.24) 

Next,  consider  the  case  of  a  source  moving  at  M  >  1/cos0o  .  In  this 
case,  the  differentiation  procedure  discussed  above  is  not  analytically 
simple  end  it  seems  more  appropriate  to  apply  some  physical  insight.  From 
Eqs.  (5.2)  and  (5.3)  the  maximum  gain  occurs  when  the  received  pulse 
duration  is  minimum.  The  received  pulse  duration  is  minimized  if  xjj  <  x^  , 
but  then  this  implies  that  the  source  at  no  time  appears  to  the  receiver  to 
be  moving  transonicelly.  Since  we  expect  the  gain  to  be  maximized  close  to 
the  transonic  region,  we  therefore  require  x^  >  x^  .  Given  this  constraint, 
tp  is  minimized  at  a  constant  value  ,  provided  x^  is  between  x^  and  xc. 
Although  (p  does  not  vary  if  x^  is  in  the  range  xt  <  x^  i  %  ,  the  pulse 
duration  in  the  frame  of  the  source,  tp  ,  increases  8s  x^  increases.  Because 
we  wish  to  minimize  the  Doppler  factor,  Eq.  (5.2)  suggests  that  we  want  to 
maximize  tp  within  the  range  that  produces  the  minimum  acceptable  tp.  We 
therefore  choose  x^  =  xc.  By  recalling  that  x^'Mc  tp  ,  and  using  Eq.  (5.8),  we 
get 


M*  =  R0  cosOq  ±  V  (Rocos0o)2  -  2R0  +1 


(5.25) 


The  negative  root  in  Eq.  (5.25)  h8S  to  be  discorded  because  Eq.  (5.8) 
indicates  that,  for  any  Mach  number,  xc  >  2  r0  cos0o  ,  so  M*  must  be  larger 
than  R0  cos0o.  By  substituting  Eq.  (5.25)  into  Eq.  (5.15b),  it  is  found  that  D* 
is  given  by 

D*  =  R0  { (cos0o/M*)  -  (sin0o/M*^M52:l)  -  1  ♦  (M*sin0o/VM**-l) )  . 

(5.26) 

The  maximum  gain  G*  is  given  by  1  /  |D*|  .  Equation  (5.25)  con  be  applied 
only  if  the  argument  of  the  square  root  is  greater  or  equal  to  zero.  It  can  be 
shown  that  this  implies 

R0  i  (  1  +  sin0o  )/cos20o  (5.27) 

If  inequality  (5.27)  is  not  satisfied,  one  has  to  use  a  different  prediction 
scheme. 

Table  III  shows  a  comparison  between  the  estimated  values  of  G*  and 
M*  supplied  by  Eq.  (5.25)  end  the  inverse  of  D*  in  Eq.  (5.26),  and  numerically 
determined  values.  For  values  of  R0  such  that 


1/cos20o  iR0i  (l+sin0o)/cos20o 


(5.28) 
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we  have  used  the  prediction  scheme  based  on  Eqs.  (5.21)  to  (5.23).  Although 
Eq.  (5.24)  is  violated  for  the  range  of  values  of  R0  (1.07  <  R0  i  1.35),  this 
method  still  yields  a  reasonable  estimate  for  the  optimum  gain  and  Mach 
number  of  a  moving  source.  This  is  due  to  the  fact  that  M*  is  a  slowly 
varying  function  of  R0.  As  indicated  in  Table  III,  good  agreement  exists  for 
a  wide  range  of  nondimensional  pulse  durations.  Hence  it  is  possible  to 
estimate  the  maximum  gain  of  a  moving  point  source  ,  based  on  a  pulse 
duration  analysis. 

B.  Line  source 

The  analysis  of  the  Doppler  shift  of  a  point  source  may  be  extended  to 
the  more  realistic  case  of  an  acoustic  array  of  finite  length  L  moving  at 
constant  velocity  v  =  Me  in  the  +x  direction.  The  approach  is  similar  to  the 
one  used  above.  The  method  of  characteristics  is  used  to  study  the  pulse 
duration  recorded  at  the  fixed  point  receiver  and  consequently,  the  Doppler 
shift  and  the  Doppler  gain  induced  by  the  motion  of  the  source. 

I.  Characteristics 


Consider  a  line  source  of  length  L  which  moves  at  constant  velocity 
perpendicular  to  its  axis  past  a  fixed  point  receiver  as  shown  in  Fig.  36.  For 
the  sake  of  simplicity,  it  is  assumed  as  in  Section  V.A  that  the  source  is 
emitting  a  rectangular  pulse  of  duration  tp.  It  is  convenient  to  define  the 
leading  edge  (first  wavelet)  end  the  trailing  edge  (last  wavelet)  of  the 
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SOURCE  c£>  M  =3 


FIGURE  36 

MOVING  LINE  SOURCE  GEOMETRY 


pulse  being  received  at  the  observation  point.  Let  tm(x)  be  the  time  of 
arrival  of  the  first  wavelet  emitted  by  the  source  when  it  was  at  a 
distance  x  from  the  origin;  and  similarly  let  t^(x)  be  the  time  of  arrival  of 
the  lest  wavelet  emitted  when  the  source  was  at  location  x.  Again  we 
choose  the  origin  of  time  to  coincide  with  the  origin  of  the  source  emission 
at  x  ■  0  .  The  first  wavelet  emanates  from  the  source  location  which  is 
closest  to  the  receiver.  Let  us  denote  by  rm(x)  the  corresponding  minimum 
distance.  Similarly  the  last  wavelet  comes  from  the  point  on  the  line 
source  farthest  from  the  receiver.  The  corresponding  distance  will  be 
denoted  rM(x).  Both  rm(x)  and  rj-j(x)  are  shown  in  Fig.  36.  Using  the 
definitions  of  rm(x)  and  rM(x),  we  have 


tm(x)  =  (x/v)  ♦  ( rm(x)/c  ) 


(5.29) 


tM(x)  =  (x/v)  +  {  rM(x)/c  ) 


(5.30) 


Let  us  define  the  time  t(x)  as  a  characteristic  time  of  the  "vertical 
diffraction"  across  the  line  source: 


t(x)  =  tM(x)  -  tm(x) 


(531) 


Combining  Eq.  (5.31)  with  Eqs.  (5.29)  and  (5.30)  yields 


t(x)  =  (  Tf-|(x)  -  rm(x)  ]  /  c 


(5.32) 


This  vertical  diffraction  time  represents  the  time  gap  between  the  arrivals 
of  the  wavelets  coming  from  the  closest  and  farthest  points  of  the  line 
source,  when  it  is  at  a  distance  x  from  the  origin.  From  elementary 
geometrical  considerations,  it  can  be  seen  thst  t(x)  increases  as  the  source 
nears  the  receiver,  and  decreases  as  the  source  moves  away  from  the 
receiver.  Note  that  a  minimum  pulse  duration  that  C8n  be  observed  8t  the 
receiver  is  always  greater  th8n  zero,  whether  the  source  moves  or  is 
stationery.  This  implies  that  the  Doppler  gain  will  always  be  finite,  even  if 
the  receiver  is  in  the  path  of  motion  of  the  source. 

If  the  line  source  is  viewed  8S  an  array  of  an  infinite  number  of  point 
sources,  then  the  received  pulse  duration  may  be  considered  as  having  two 
contributing  factors:  the  pulse  duration  associated  with  a  moving  point 
source,  and  the  vertical  diffraction  time.  The  time  coordinate  straining  due 
to  the  motion  of  a  point  source  was  discussed  in  Section  V.A.  It  remains 
then  to  use  these  results  to  determine  the  characteristics  of  a  line  source. 
The  line  source  geometry  may  be  divided  into  three  cases,  each  of  which  is 
treated  separately  below,  the  short  array,  the  array  of  length  comparable 
to  the  receiver  depth,  and  the  long  array. 

First  consider  8  short  source  as  defined  by  L/h  <  1  ,  where  h  is  the 
depth  of  the  receiver.  By  referring  to  Fig.  36,  one  can  see  that  the  nearest 
point  to  the  receiver  is  the  point  at  the  deep  end  of  the  source,  and  that  the 
farthest  point  from  the  receiver  is  the  point  located  at  z  ■  0  .  If  Gm(x)  and 
0j-j(x)  are  the  angles  formed  between  the  direction  of  motion  and  the 


vectors  to  the  receiver  from  the  nearest  end  farthest  points  on  the  source, 
respectively,  when  it  is  located  at  position  x,  then 


i 

d/dx  lrm(x)]  =  -  cos  Om(x)  ,  (5.33) 

i  and 

d/dx  Irj-j(x))  ■  -  cos  9^(x)  (5.34) 

Equations  (5.33)  and  (5.34)  ere  valid  even  if  the  source  has  passed  over  the 
receiver.  Combining  Eqs.  (5.33)  and  (5.34)  with  Eqs.  (5.29)  and  (5.30)  yields 


d/dx(tm)  *  [  1  -  M  cos  0m(x)  1  /  v 


and 


i 

i 


d/dx  (tfl)  =  1 1  -  M  cos  Ofl(x)  1  /  v 


(5.35) 


(5.36) 


From  Eqs.  (5.35)  and  (5.36),  the  characteristics  of  a  short  array  may  be 
drawn.  As  in  the  point  source  case,  the  characteristics  are  in  general 
curved. 

Next,  consider  an  array  whose  length  is  comparable  to  the  depth  of 
the  receiver.  Such  a  source  is  defined  by  1  i  L/h  <  2  .  In  this  case,  the 
shortest  distance  between  the  srray  and  the  receiver  is  the  distance  from 
the  array  point  located  at  the  same  depth  as  the  receiver.  Let  us  denote  this 
distance  by  r  j^(x)  end  the  angle  to  the  receiver  as  Oj^(x).  Then 
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rm(x)  *  rjte) 


(5.37) 


Under  these  conditions. 


d/dx  (rm)  =  -  cos  0jjx)  ,  (5.38) 


and 


d/dx  (r^j)  =  -  cos  0j-j(x)  (5.39) 

This  leads  to  the  following  slopes  for  the  characteristics: 

d/dt  (tm)  =  ( 1  -  M  cos  0±(x) )  /  v  ,  (5.40) 


and 

d/dx  (tpj)  =  { 1  -  M  cos  ©|-|(x) )  /  v  (5.41) 


Note  that  Eqs.  (5.40)  and  (5.41)  are  identical  to  Eqs.  (5.35)  and  (5.36)  with 
the  transformation  0m  -»  0j. 

Finally  consider  the  case  of  a  long  line  source  (  L/h  i  2  ) .  Now  the 
shortest  distance  to  the  receiver  is  that  from  the  source  point  located  at 
the  seme  depth,  and  the  longest  distance  is  from  the  deepest  source  point. 
Retaining  the  definition  of  ©j-|  as  the  angle  associated  with  the  vector  from 

118 


.'Iv'lsV 


the  farthest  source  point  to  the  receiver,  Eqs.  (5.40)  and  (5.41)  are  still 
valid. 

Figure  37  shows  the  characteristics  of  a  line  source  by  showing  the 
characteristics  associated  with  motion  of  the  nearest  and  farthest  points 
from  the  receiver.  The  characteristics  for  each  point  resemble  Fig.  33 
closely,  8nd  it  is  possible  for  each  source  location  on  the  line  source  to 
define  xt  and  xc.  For  a  given  Mach  number,  the  time  separation  between  the 
two  characteristics  corresponding  to  the  two  location  is  t(x). 

An  interesting  phenomenon  occurs  when  1/cos^  <  M  <  l/cosO^  .  In  this 
case,  time  reversal  occurs  on  wavelets  coming  from  the  point  source 
closest  to  the  receiver  (i.e.,  the  slope  of  t^x)  is  negative),  while  no  time 
reversal  affects  the  wavelets  coming  from  the  farthest  point  of  the  line 
source  (i.e.,  the  slope  of  tj^(x)  is  still  positive).  This  "partial  time 
inversion"  region  is  the  smooth  transition  from  apparently  subsonic  to 
apparently  supersonic  behavior,  as  seen  by  the  receiver.  It  is  therefore  the 
transonic  region.  This  phenomenon  can  also  be  understood  as  a  phase 
inversion  affecting  only  wavelets  coming  from  a  certain  part  of  the  line 
source.  This  suggests  the  possibility  of  reshuffling  the  wavelets  between 
the  source  end  the  receiver.  If  a  source  is  emitting  a  signal  which  varies 
with  its  depth,  for  example,  the  sequence  A,  B,  C,  D,  E,  then  it  seems 
possible  to  receive  the  sequence  A,  D,  C,  B,  E.  This  may  possibly  have  some 
applications  in  signal  encryption. 


FIGURE  37 

CHARACTERISTICS  FOR  A  MOVING  LINE  SOURCE 


ARL  UT 
AS-84-694-P 
YHB - GA 
8-6-84 


120 


2.  Pulse  duration 


As  discussed  in  the  case  of  a  point  source,  it  is  convenient  to  describe 
the  emission  at  the  source  8S  a  rectangular  pulse  of  duration  cp  during 
which  the  source  travels  a  distance  Xfj  ■  v  tp  .  Since  t(x)  has  been  defined 
as  the  time  duration  of  the  disturbance  generated  by  the  source  when  it 
was  at  a  distance  x  from  the  origin,  t(Oix<x^)  subsequently  represents  the 
time  gap  between  the  arrivals  of  ell  wavelets  emitted  during  the  motion  of 
the  source.  Therefore  the  received  pulse  duration  tp  is  the  projection  on 
the  time  axis  of  t(0<x<xN).  As  usual,  the  Doppler  factor  is  the  straining  of 
the  time  coordinate  due  to  the  motion  of  the  source.  It  is  given  by 


tp  ♦  t(0) 


where  t(O)  is  the  vertical  diffraction  time  at  the  origin.  Note  that  Eq.  (5.42) 
differs  from  Eq.  (5.2)  for  the  point  source  because  the  finite  extent  of  the 
line  source  will  cause  the  received  pulse  to  spread  in  time. 

The  received  pulse  duration  can  be  expressed  in  analytical  form  by 
examining  the  characteristics  to  determine  their  projection  on  the  time 
axis.  As  in  the  point  source  case,  a  distinction  must  be  made  between  a  line 
source  which  appears  to  the  receiver  to  be  moving  subsonically  and  a  line 
source  which  seems  to  be  moving  supersonically. 

If  M  <  1/cos0m(O)  then  the  entire  line  source  appears  to  the  receiver 


to  be  moving  subsonically.  Another  way  of  stating  this  is  to  say  that  the 
characteristic  for  each  point  on  the  line  source  h8s  a  positive  slope.  For  8 
moving  line  source  of  this  type,  the  earliest  arrival  of  sound  comes  from 
the  first  wavelet  emitted  by  the  point  source  nearest  to  the  receiver.  The 
latest,  sound  comes  from  the  last  wavelet  emitted  by  the  point  source 
farthest  from  the  receiver.  Hence 

tp  =  -  tm(0)  (5.43) 

If  l/cos0m(O)  <  M  <  l/cos0^j(O)  ,  the  characteristic  t^(x)  has  always  a 
positive  slope,  whereas  the  slope  of  tm(x)  is  negative  for  x  <  xtm  and 
positive  for  x  >  xtm.  Here  x^m  denotes  the  location  at  which  the  nearest 
point  source  8ppe8rs  to  be  moving  8t  transonic  velocity.  It  may  be  found 
from  Eq.  (5.7).  It  is  therefore  necessary  to  consider  two  cases  in  order  to 
evaluate  the  received  pulse  duration  tp.  The  cases  depend  on  the  distance 
Xjij  traveled  by  the  source.  If  x^  <  x^m  ,  the  first  wavelet  that  reaches  the 
receiver  arrives  at  time  tm(x^)  (see  Fig  37).  Similarly,  the  last  wavelet 
that  reaches  the  receiver  arrives  at  time  t^(xM) .  However,  if  xN  >  x^m  ,  the 
time  of  arrival  of  the  first  wavelet  is  tm(xtm)  and  the  last  wavelet  still 
reaches  the  receiver  at  t^(x^).  The  received  pulse  duration  is  therefore 
given  by  Eqs.  (5.44)  8nd  (5.45). 


XN  <  xtm  - 


tD  =  "  WXN^ 


(5.44) 


If  XN2Xlm, 


tp  =  -  tm(Xtm) 


(5.45) 


if  II  >  1/cos0j-j(O) ,  then  the  entire  line  source  appears  to  the  receiver 
to  be  moving  supersonically  at  time  t=0.  This  implies  that  the 
characteristic  of  each  point  source  starts  out  with  a  negative  slope.  Two 
different  esses  depending  on  x^,  must  be  investigated.  As  shown  in  Fig.  37, 
if  xN  <  x^m,  the  time  of  arrival  at  the  receiver  of  the  first  wavelet  is  given 
by  tm(xN),  whereas  the  time  of  arrival  of  the  last  wavelet  is  tM(0). 
However,  if  xN  >  xtm  ,  Fig.  37  shows  that  the  times  of  arrival  of  the  first 
end  lest  wavelets  ere,  respectively,  tm(x^m)  end  ttt(x^).  The  received  pulse 
duration  is  therefore  given  by  Eqs.  (5.46)  and  (5.47). 


If  Xfyj  <  x^m  , 


tp =  *ii^)  ~ 


(5.46) 


If  Xfy|  >  X^m  , 


tp  =  tM(XN)  -  tm(X|m) 


(5.47) 


Note  that  in  ell  of  these  cases,  subsonic  and  supersonic,  the  received  pulse 
duration  is  greater  then  the  vertical  diffraction  time  at  the  origin,  t(O). 
From  Eq.  (5.42)  and  the  definition  of  the  Doppler  gain,  Eq.  (5.3),  we  may 
determine  an  upper  bound  for  the  gain  of 


(5.48) 


6  <  1  ♦  {  Cp  /  t(0)  ] 

Hence,  e  line  source  moving  at  any  Mach  number  will  always  produce  a 
finite  gain  at  the  receiver.  For  a  point  source  t(O)  =  0,  and  we  recover  the 
result  that  6  <  <»  ,  i.e.,  that  the  Doppler  gain  is  finite. 

A  computer  program  has  been  written  to  evaluate  the  gain  associated 
with  an  acoustic  line  source  of  finite  length  moving  at  constant  velocity. 
Figures  38(8)  and  38(b)  are  three  dimensional  plots  showing  the  Doppler 
gain  as  a  function  of  both  Mach  number  and  pulse  duration.  The  length  of  the 
array  was  set  to  L  =  0.01  m  in  Fig.  38(a),  and  L  *  0.1  m  in  Fig.  38(b).  In  each 
case,  the  receiver's  coordinates  were  x  *  4  m,  y  =  0  m,  and  z  ■  4  m,  (see 
Fig.  36  ).  The  pulse  duration  varied  between  tp  *  0.02  ms  end  tp  ■  2  ms;  in 
nondimensionel  form,  this  means  that  T0  *  ctp/r0  varied  from  0.005  to 
0.525. 

As  expected  from  the  previous  discussion,  the  gain  associated  with 
the  motion  of  the  source  is  always  finite.  It  is  also  interesting  to  note  that 
the  finite  length  of  the  array  introduces  some  "damping"  into  the  response. 
Figure  38(e)  shows  a  noticeably  sharper  peak  in  the  gain  than  Fig.  38(b) 
which  is  at  the  same  scale.  This  csn  be  explained  by  noticing  th8t  the  longer 
the  source,  the  longer  its  "vertical  diffraction  time",  and  therefore  the 
more  time  spreading  in  the  response,  and  the  lower  the  gain. 


C.  Conclusions 


It  has  been  shown  that  the  Doppler  gain  associated  with  a  point 
source  moving  at  constant  velocity  in  a  line  past  a  point  receiver  is  always 
finite,  even  when  the  source  is  moving  at  transonic  velocity.  The  only 
exception  is  when  the  source  is  moving  straight  at  the  receiver.  In  this 
case,  it  is  necessary  to  include  the  diffraction  effects  induced  by  the  finite 
dimensions  of  any  reel  source,  in  order  to  argue  that  the  Doppler  gain  is 
still  bounded. 

A  simple  analytical  prediction  scheme  has  been  used  to  find  the 
maximum  gain  8nd  the  optimum  Mach  number  for  a  moving  point  source.  The 
results  are  valid  for  the  case  of  a  source  passing  over  the  receiver.  In 
addition,  simple  expressions  have  been  developed  by  examining  the  source 
characteristics  for  the  Doppler  shifts  associated  with  a  source  moving  at 
constant  velocity. 

The  results  for  a  moving  point  source  were  extended  to  include  the 
case  of  a  moving  line  source.  The  analysis  of  the  Doppler  shift  for  a  line 
source  of  finite  length,  moving  at  constant  velocity  in  a  straight  line,  ere 
very  much  like  those  of  a  point  source.  The  Doppler  gain  is  finite  for  all 
Mach  numbers,  including  transonic, end  the  Mach  number  which  produces  the 
maximum  gain  is  a  function  of  the  nondimensional  distance  Ro=r0/ctp.  It 
was  also  shown  that  the  Doppler  gain  of  an  acoustic  array  moving  at 
constant  velocity  decreases  with  the  array  length.  Finally,  it  was  shown 
that  even  in  the  case  of  a  line  source  moving  at  transonic  velocity  directly 
toward  a  point  receiver,  the  Doppler  gain  is  always  finite. 
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Although  the  present  analysis  is  carried  out  for  the  case  of  a 
transient  pulse  emission  of  finite  duration  at  the  source,  the  results  can 
easily  be  extended  to  the  case  of  any  arbitrary  source  function.  In  the 
analysis,  the  pulse  duration  plays  the  role  of  a  characteristic  time,  such  os 
the  period  for  the  case  of  monochromatic  signals. 
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CHAPTER  VI 


CONCLUSIONS  AND  SUGGESTIONS 

The  thermoacoustic  mechanism  of  sound  generation  by  a  laser  pulse 
beamed  into  water  has  been  investigated,  both  theoretically  end 
experimentally,  for  both  stationary  and  moving  sources.  The  theoretical 
model  is  based  on  the  linear  Westervelt-Larson  equation8  describing 
thermoacoustic  sound  generation  for  an  inviscid  medium  containing  a  heat 
source  such  as  a  laser  beam  in  water.  A  time  domain  analysis  is  developed 
in  order  to  find  the  pressure  field  radiated  by  the  thermoecoustic  source. 
The  acoustic  pressure  is  expressed  as  a  convolution  type  summation 
between  the  impulse  response  of  the  system  and  the  thermoecoustic  source 
strength  (time  derivative  of  the  laser  intensity).  Although  this  approach 
requires  in  most  esses  the  use  of  a  computer,  it  offers  several  advantages: 
(1)  It  is  valid  in  the  neerfield  of  the  source  and  this  has  been  verified 
experimentally.  (2)  It  is  also  valid  for  any  source  velocities,  because  the 
Doppler  shift  associated  with  the  motion  of  the  source  appears  naturally 
within  the  "pseudo-convolution".  (3)  Although  the  analysis  presented  here  is 
restricted  to  the  common  case  of  a  laser  beam  with  a  Gaussian  cross 
section  intensity  distribution,  with  an  exponential  shading  along  the  beam 
penetration  axis,  the  model  can  be  extended  to  any  laser  intensity  profile. 
(4)  The  impulse  response  of  a  thermoecoustic  source  is  a  useful  tool  in 
analyzing  the  acoustic  response  to  thermoecoustic  transient  phenomena 
such  as  extremely  short  laser  pulses  or  Mach  waves  produced  by  a  laser 


beam  moving  at  high  velocity.  It  turns  out  that  these  transient  signals  are 
among  the  most  promising  ways  to  generate  reasonably  high 
thermoacoustic  levels. 

Experimental  results  were  obtained  with  a  laser  system  providing  up 
to  5  joules  of  energy  over  a  pulse  duration  of  approximately  1  ms,  during 
which  the  intensity  was  modulated  at  a  single  frequency  between  5  and  80 
kHz.  Most  of  the  experimental  results  were  obtained  with  a 
Neodymium-Glass  laser  (optical  wavelength  1.06  pm),  so  that  the 
thermoacoustic  source  length  in  water  was  fairly  short  (~0.1m).  Some 
experimental  results  are  also  given  for  the  case  of  a  long  Mm) 
thermoecoustic  source  which  was  obtained  by  using  a  ruby  laser  (optical 
wavelength  0.6943  pm). 

The  main  experimental  results  can  be  grouped  into  four  categories: 
pressure  waveforms,  directivity  patterns,  sound  level  dependence  on 
source  velocity,  8nd  spreading  curves.  Source  velocities  up  to  Mach  2  were 
investigated  with  special  attention  to  the  case  of  transonic  velocity.  In 
general,  the  experimental  results  ere  in  good  agreement  with  theoretical 
predictions,  and  this  confirms  the  validity  of  the  Westervelt-Larson  model 
even  for  supersonic  velocities. 

Finally,  the  Doppler  shift  of  an  acoustic  source  moving  at  transonic 
velocity  was  analyzed  theoretically  by  taking  into  account  the  time 
dependence  of  the  angle  of  observation  between  the  source  and  the  receiver. 
It  was  shown  that  this  time  dependence  removes  the  singularity  which 
occurs  in  the  standard  theory,  when  the  source  is  moving  at  transonic 
velocity.  It  was  also  shown  that,  in  the  case  of  a  line  source,  the  transition 


between  subsonic  and  supersonic  regimes  is  a  region  where  only  wavelets 
coming  from  a  certain  part  of  the  source  are  time-inverted,  whereas  the 
rest  of  the  wavelets  are  not.  This  partial  time  inversion  suggests  the 
possibility  of  reshuffling  the  wavelets  in  8  signal  8nd  this  m8y  have  some 
applications  in  signal  encryption.  Future  work  on  this  topic  should  also 
include  diffraction  effects  due  to  the  finite  size  of  the  source,  and 
non-rectilinear  motion  of  the  source  to  analyze  p8th  curvature,  source 
acceleration,  caustics  and  focusing  effects. 

The  results  reported  in  this  study  show  that,  in  general,  the 
experimental  results  are  in  good  agreement  with  the  theoretical 
predictions.  The  major  discrepancy  was  found  to  be  that  the  theoretical 
model  underestimates  by  a  significant  amount  the  measured  acoustic  level 
when  the  source  is  moving  at  transonic  velocity.  However,  the  predicted 
waveform  is  in  reasonable  agreement  with  the  measured  waveform.  It  was 
found  experimentally  that  by  moving  the  source  at  transonic  velocity,  a  gain 
of  about  20  dB  could  be  expected,  and  a  very  strong  broadband  pressure 
transient  was  generated.  A  sound  pressure  level  of  152  dB  re  1  yPa  at  1  m 
was  measured  for  such  a  M8Ch  wave.  This  high  level  shows  that  the  low 
frequency  components  of  the  broadband  signal  can  propagate  at  fairly  large 
ranges  before  being  absorbed.  The  theoretical  model  predicts  that  in  the 
vicinity  of  Mach  one,  diffraction  effects  become  predominant  and  actually 
tend  to  ruin  the  efficiency  of  the  process  Although  diffraction  effects 
were  clearly  observed  around  Mach  one,  it  is  not  fully  understood  why  the 
theoretical  model  underestimates  the  measured  level.  A  possible 
explanation^?  is  that  nonlinear  effects,  such  as  the  temperature 


dependence  of  the  thermodynamic  coefficient  p,  become  important  and 
cannot  be  neglected. 

The  surprisingly  high  thermoacoustic  levels  measured  at  transonic 
source  velocity  indicate  that,  at  least  for  sonar  applications,  future  work 
should  focus  on  laser  induced  Mach  waves.  It  h8S  been  suggested12,26,34 
that  the  acoustic  disturbance  generated  by  a  transonic  thermoacoustic 
source  is  identical  to  the  one  generated  by  8  stationary  source  of  the  same 
optical  energy  but  delivered  over  an  extremely  short  period  of  time. 
Although  this  is  approximately  true  for  an  observer  very  far  from  the 
source,  there  is  still  8  major  advantage  in  using  a  moving  source  instead  of 
a  stationary  source:  The  region  of  the  medium  which  absorbs  the  heat 
produced  by  the  laser  beam  is  spread  over  the  laser  path  length  whereas  in 
the  case  of  a  stationary  source  it  is  localized  to  the  laser  beam  cross 
section.  This  means  that  the  maximum  energy  density  acceptable  before 
surface  evaporation  or  explosive  boiling  occurs  in  the  medium  is  much 
higher  in  the  case  of  a  moving  source  than  for  a  stationary  source.  In  f8ct 
it  is  increased  by  a  factor  vtp/2a  ,  where  v  is  the  source  velocity,  tp  the 
laser  pulse  duration,  and  2a  the  beam  diameter. 

There  are  several  reasons  for  choosing  the  thermoacoustic 
mechanism  of  laser  sound  generation  over  other  mechanisms  such  as 
surface  evaporation,  explosive  boiling,  or  optical  breakdown.  Although 
these  other  mechanisms  are  intrinsically  more  efficient  than  the 
thermoacoustic  mechanism,2-35'3^37  they  cannot  produce  narrow 
directional  sound  beams.  Moreover  these  energy  conversion  mechanisms  are 
still  not  very  well  understood.  It  seems  therefore  safer  to  use  the 


thermoacoustic  process  and  to  push  its  upper  limit  by  scanning  the  laser 
beam  at  high  velocities.  The  strongest  acoustic  signal  that  can  be  generated 
by  a  moving  thermoacoustic  source  occurs  when  the  laser  beam  is  moving 
8t  transonic  velocity.  However  it  should  be  emphasized  that  in  this  case  the 
recorded  acoustic  signal  has  8  very  short  duration  so  that  the  information 
content  in  the  signal  is  limited.  Also  the  high  frequencies  of  such  8  pulse 
8re  attenuated  very  quickly  in  the  water  as  the  signal  propagates.  This 
problem  can  be  overcome  by  find? ng  judicious  W8ys  of  depositing  the  he8t 
source  in  the  medium. 

A  possible  design  which  has  been  investigated  in  detail  both 
theoretically  and  experimentally  by  Soviet  physicists,'58  consists  of  using  a 
stationary  laser  source  delivering  a  periodic  train  of  very  short  laser 
pulses,  separated  in  time  by  an  arbitrary  delay  T0.  This  is  achieved  for 
instance  by  Q-switching  the  laser  source.  In  this  design,  a  strong  end 
directional  acoustic  disturbance  of  m8in  frequency  component  1/T0  ,  is 
expected  to  propagate  at  long  distances  under  water.  As  shown  in 
reference  38,  this  design  increases  by  a  factor  of  about  four  the 
thermoacoustic  intensity  over  that  of  a  standard  monochromatic 
modulation  system.  However  the,  upper  limit  of  the  energy  density 
acceptable  to  stay  within  the  range  of  the  thermal  mechanism  of  sound 
generation  is  limited  by  the  fact  that  the  source  is  stationary. 

A  more  sophisticated  design  has  recently  been  suggested  by  Pierce  end 
Hsieh.^  it  consists  of  splitting  the  laser  beam  into  several  beams, 
separated  by  an  arbitrary  distance  \0  on  the  surface  of  the  water,  and  all 
moving  at  a  slightly  supersonic  velocity.  The  theoretical  analysis,  which  is 


based  on  a  very  general  acoustic,  pumping  principle/*'  predicts  that  this 
configuration  allows  a  very  strong  quasi-plane  wave  of  main  frequency 
component  c/X0  ,  to  propagate  very  f8r  under  water  in  a  very  collimated 
beam.  The  motion  of  the  source  would  allow  very  high  electromagnetic 
energy  levels  to  be  deposited  into  the  water,  within  the  range  of  the 
thermoacoustic  sound  generation  mechanism,  and  without  surface 
evaporation  or  explosive  boiling  deteriorating  the  directional  properties  of 
the  sound  beam.  At  sufficiently  high  energy  levels,  which  would  produce 
other  mechanisms  if  the  source  were  stationary,  nonlinear  effects  due  to 
the  temperature  dependence  of  the  thermodynamic  parameters  (p,cp)  could 
even  enhance  the  thermoecoustic  conversion  efficiency.27  The  experimental 
investigation  of  such  a  design  seems  to  be  the  logical  extension  of  the 
work  on  laser  thermoacoustic  sources.  If  experimental  results  confirm  the 
usefulness  of  this  design,  it  would  certainly  place  thermoacoustics  within 
the  range  of  some  sonar  applications. 

Laser  performances  are  of  course  a  determinant  factor  in  the 
capabilities  of  a  thermoacoustic  sonar.  Fortunately,  laser  technology  has 
improved  considerably  over  the  last  two  decades,  and  it  is  expected  to  do  so 
in  the  near  future.  It  is  now  possible  to  find  pulsed  solid  state  lasers  that 
emit  a  single  pulse  of  several  hundred  joules  over  several  milliseconds. 

Another  type  of  laser  that  may  be  suitable  for  some  underwater 
applications  of  thermoacoustics,  is  the  fairly  new  excimer  laser. This 
laser  delivers  very  short  pulses  (4-50  ns)  in  the  ultraviolet  spectrum 
(0.18  to  0.3  pm)  at  a  very  high  repetition  rate  so  that  high  mean  powers  are 


achievable  The  optical  frequency  can  be  downshifted  to  wavelengths 
suitable  for  thermoecoustic  purposes  (0.5  to  1  ym)  by  a  very  efficient 
method  known  as  the  stimulated  Raman  scattering.  Moreover  these  lasers 
are  reliable,  easy  to  operate,  very  pure  and  they  can  produce  a  very  small 
beem  diameter  so  that  very  high  energy  densities  can  be  produced  fairly 
easily.  However,  the  pulse  duration  is  too  short  to  take  advantage  of  any 
motion  of  the  source. 

Infrared  lasers  such  as  the  CO2  laser,  or  the  NH3  laser,  are  8mong 
the  most  efficient  CW-lesers.  In  certain  applications,  the  CO2  laser  can 
deliver  up  to  500  kW  continuous.42  Unfortunately  the  spectral  region  of 
these  lasers  is  too  f8r  away  from  the  region  of  interest  for  thermoecoustic 
applications.  (Stimulated  Raman  scattering  may  be  used  only  to  increase  the 
lasing  wavelength,  not  to  decrease  it.) 

However  recent  developments  on  chemically  pumped  lasers  seem  to 
indicate43  that  8  continuous  high-power  laser  of  wavelength  1.3  ym  could 
be  developed  ,  and  this  laser  would  certainly  be  perfectly  suitable  for 
underwater  applications  of  thermoacoustics. 

In  summary,  it  seems  that  the  laser  technology  is  improving  so 
rapidly  that  in  a  not  too  distant  future,  one  could  conceivably  design  a 
l8ser-induced  sonar  based  on  the  principles  of  thermoecoustic  generation 
of  sound. 
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APPENDIX  A 


NONUNIFORMITY  OF  THE  LASER  SCANNING  VELOCITY 
AND  OBLIQUE  INCIDENCE  OF  THE  LASER  BEAM 

In  this  study,  it  is  assumed  that  the  laser  beam  remains 
perpendicular  to  the  surface  of  the  water  during  the  laser  pulse  duration. 
This  is  just  an  approximation  to  the  real  physical  situation  encountered  in 
the  experiment,  because,  as  indicated  in  Fig.  1,  the  rotating  mirror  imposes 
during  the  laser  pulse  duration  a  deflection  angle  <Kt)  which  is  in  general 
not  zero.  This  has  two  consequences:  first,  the  laser  scanning  velocity  v 
will  not  be  rigorously  constant  in  time;  and  second,  the  laser  beam  will  be 
refracted  as  it  enters  the  water  because  the  speed  of  light  is  different  in 
air  and  water.  The  objective  of  this  appendix  is  to  quantify  the  importance 
of  these  two  effects. 

It  can  be  shown  from  elementary  geometrical  considerations  that 
the  velocity  v  of  the  laser  spot  on  the  water  surface,  is  given  by 


v(t)  = 


cos^(t) 


where  h  is  the  distance  between  the  rotating  mirror  and  the  water  surface, 
Q  is  the  angular  velocity  of  the  rotating  mirror,  and  <|>(t)  is  the  angle 
between  the  laser  beam  and  the  vertical  axis  as  indicated  in  Fig.  1.  The 
relative  error  Av/v  can  therefore  be  written  as 


Av  /  v  =  I  -  cos"2  Wt) 


(A2) 


which  shows  that  the  maximum  error  occurs  when  Wt)  is  maximum  If  we 
assume  that  the  laser  beam  is  perpendicular  to  the  water  surface  at  the 
mid-point  in  the  laser  pulse  duration  tp,  so  that  Wo)  *  Wcp),  then  it  can 
be  shown  that  Wt)  is  maximum  at  t=0,  and  it  takes  the  value 

Wo)  =  tan"1  (  M  c  tp  /  2h)  ,  (A3) 

where  M=v/c  is  the  Mach  number  of  the  thermoacoustic  source,  and  c  is  the 
sound  speed  in  W8ter.  Equations  (A2)  and  (A3)  can  be  combined  to  find  the 
maximum  relative  error  in  the  nonuniformity  of  the  laser  spot  on  the 
surface  of  the  W8ter. 

The  second  effect  induced  by  the  non-perpendicul8rity  of  the  laser 
beam  on  the  W8ter  surface  is  a  refraction  of  the  beam  as  it  penetrates  into 
the  water.  This  refraction  effect  is  described  by  Snell's  l8w: 

sin  W)  /  c5|r  =  sin  *«)  /  cwster  _  (A4) 


where  Wt)  is  the  refracted  angle  between  the  laser  beam  and  the  vertical, 
in  the  water;  Wt)  is  the  incident  angle  of  the  laser  beam  on  the  water 
surface;  and  cW0^er  and  c0jr  are,  respectively,  the  speed  of  light  in  water 
and  in  air.  The  index  of  refraction  n,  defined  as  the  ratio  of  c0jf-  over 
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cwetep  is  for  both  optical  wavelengths  (0.6943  ym  and  1.06  ym)  used 
in  the  experiment.  Therefore  the  tilt  angle  Wt)  of  the  laser  beam  in  W8ter 
is  given  by 


sin  t(t)  =  0.75  sin  Wt)  ,  (A5) 

8nd  the  maximum  tilt  angle  occurs  when  Wt)  is  maximum,  i.e.,  at  t=0. 
Combining  Eqs.  (A3)  and  (A5)  gives  an  estimate  of  the  maximum  tilt  angle  of 
the  thermoacoustic  source,  as  a  function  of  Mach  number: 

♦(0)  =  sin-1  {  0.75  sin  [  tan"*  (  M  c  cp  /  2  h  ) )  }  .  (A6) 

The  following  table  gives  some  numerical  estimates  of  WO),  WO),  in 

degrees,  and  Av/v  as  a  function  of  Mach  number  of  the  laser  spot  on  the 
surface  of  the  water.  It  is  assumed  that  the  speed  of  sound  in  water  is 
c  =  1486  m/s  ,  the  height  between  the  rotating  mirror  and  the  water 

surface  is  h  =  4.06  m  ,  and  that  the  laser  pulse  duration  is  tp  =  1  ms. 
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APPENDIX  B 


FREQUENCY  RESPONSE  AND  BEAM  PATTERNS  OF  THE  HYDROPHONE 

The  hydrophone  used  in  the  study  of  thermoacoustic  sources  was,  in 
most  cases,  the  USRD  type  H-56  hydrophone  (serial  number  1 14).  This 
hydrophone  is  a  high  sensitivity  and  low  self-noise  device,  so  that  weak 
acoustic  signals  can  be  detected.  The  most  important  electroacoustic 
characteristics  of  the  H-56  can  be  found  in  reference  44. 

Figure  39  shows  the  sensitivity  of  the  H-56  as  a  function  of 
frequency,  from  10  kHz  to  90  kHz.  Figure  40  (a)-(g)  represents  the 
sensitivity  of  the  H-56  as  a  function  of  angle  in  the  x-z  plane  (see  Fig.  22), 
for  several  frequencies:  10  kHz,  20  kHz,  30  kHz,  40  kHz,  50  kHz,  75  kHz,  and 
100  kHz. 

The  sensitivity  of  the  receiving  hydrophone  as  a  function  of  (received) 
frequency  and  angle  was  taken  into  account  in  the  experimental  analysis  of 
thermoacoustic  sources. 


DIRECTIVITY  PATTERN 

FOR:  HYDROPHONE 

PATTERN  AXIS:  VERTICAL  DATE:  JAN  29 
PROJ:  TR  225  HYDR:  H56-114 

FREQ.  Hz:  .  20000  SENS: 

XMIT  GATE  MSEC:  .  5  REC  GATE  MSEC:  .  4 

DRIVE:  FILTER:  5-100K 


FIGURE  40  (b) 
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DIRECTIVITY  PATTERN 


FOR:  HYDROPHONE 

PATTERN  AXIS:  VERTICAL  DATE:  jan  29 


PROJ:  IR  225 
FREQ.  Hz:  .40000 

XMIT  GATE  MSEC:  .  5 
DRIVE: 


HYDR:  H56-114 

SENS: 

REC  GATE  MSEC:  .4 
FILTER:  5-100K 


FIGURE  40  (d) 
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DIRECTIVITY  PATTERN 

FOR:  HYDROPHONE 

PATTERN  AXIS:  VERTICAL  DATE:  JAN  29 
PROJ:  TR  225  HYDR:  H56-114 

FREQ.  Hz:  .  50000  SENS: 

XMIT  GATE  MSEC:  -5  REC  GATE  MSEC:  .4 

DRIVE:  FILTER:  5-10QK 
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FIGURE  40  (e) 


ARL:UT 
AS-85-661 
YHB - GA 
6-17-85 


DIRECTIVITY  PATTERN 

FOR:  HYDROPHONE 

PATTERN  AXIS:  VERTICAL  DATE:  JAN  29 
PROJ:  TR  225  HYDR:  H56-114 

FREQ.  Hz:  .  100000  SENS: 

XMIT  GATE  MSEC:  .  5  REC  GATE  MSEC:  .  4 

DRIVE:  FILTER:  5-100K 

E  (Vrms):  .10.0 
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FIGURE  40  (g) 
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APPENDIX  C 


INTERFACE  BETWEEN  THE  OSCILLOSCOPE  AND  THE  COMPUTER 

The  Nicoiet  4094  digital  oscilloscope  was  interfaced  to  the 
CDC-CVBER  computer  so  that  recorded  waveforms  could  be  transferred  to 
the  mainframe  computer  and  compared  with  theoretical  predictions.  This 
appendix  summarizes  the  procedure  for  the  data  transfer. 

The  output  labeled  "TERMINAL"  on  the  back  panel  of  the  Nicoiet  should 
be  connected  to  a  CYBER  terminal.  The  switch  labeled  "GPIB-RS  232", 
located  also  on  the  back  panel  of  the  Nicoiet,  should  be  placed  in  the 
position  "RS  232".  The  actual  transfer  of  data  is  done  via  that  RS  232 
interface.  The  software  for  the  data  transfer  is  a  program  called  "NICL" 
(binary  version.  "NICOLT”),  which  was  written  by  D.  Gallant-Offer.  This 
program  must  be  run  interactively  from  the  CYBER  terminal  connected  to  the 
oscilloscope.  To  run  the  program,  type: 

/GET,NICOLT/UN=GQ 

/NICOLT 

The  message  "ENTER  COMMAND"  will  appear.  There  are  six  possible 
commands. 


Dli  transfer  waveform  from  diskette  to  Nicolet's  memory 

MD  transfer  waveform  from  Nicolet's  memory  to  diskette 

MC  transfer  waveform  from  Nicolet's  memory  to  CYBER 

DC  transfer  waveform  from  diskette  to  CYBER 

WF  display  data  about  waveforms  in  the  Nicolet's  memory 

QU  quit  the  program 

An  invalid  command  causes  the  message  "INVALID  COMMAND"  to  appear 
before  it  again  prompts  for  a  command.  A  carriage  return  without  command 
aborts  the  program  with  an  end-of-file  error.  Each  command  prompts  the 
information  it  needs.  Typing  while  the  Nicolet  is  communicating  with  the 
CYBER  usuelly  aborts  the  program  and  requires  that  the  Nicolet  be  reset.  (To 
reset  the  Nicolet,  switch  back  and  forth  the  RS  232/6PIB  switch  on  the  back 
panel  of  the  Nicolet).  While  the  program  is  running,  it  will  display  error 
codes  from  the  Nicolet.  An  error  code  of  zero  means  that  no  error  occurred. 
The  explanation  for  the  other  error  codes  is  given  on  pp.  13-60  (Table  13-25) 
in  the  Nicolet  4094  waveform  analysis  manual. 

The  following  messages  appear  before  the  transfer  of  data. 


Message 


Explanation 


ENTER  COMMAND 
WAVEFORM  NUMBER 

RECORD  NUMBER 
MEMORY  SECTION 
0 
1 

2 

3 

4 

5 

6 


enter  one  of  the  6  commands  (usually  MC) 

enter  the  number  of  the  waveform  to  be  transferred 

(usually  1). 

enter  diskette  record  number  to  read  from  or  save  to 

enter  the  segment  of  memory  to  read  from  or  save  to 

all  of  memory  (ALL) 

first  quarter  (Q1) 

second  quarter  (Q2) 

third  quarter  (Q3) 

fourth  quarter  (Q4) 

first  half  (HI) 

second  half  (H2) 


STARTING  POINT  first  point  to  be  transferred  from  the  waveform 

(1  to  15871) 

TOTAL  NUMBER  OF  POINTS  total  number  of  points  to  be  transferred 


FILENAME 


(1  to  15871) 

name  of  file  to  write  data  to. 


The  output  of  "NICQLT"  is  a  file  on  the  CYBER  in  free  format,  called 
“filename".  It  contains  the  values  of  the  data  points  of  the  waveform,  in 
integer  values.  In  order  to  convert  this  set  of  data  into  their  real  values,  a 
program  called  “CONVERT"  (“CVRT"  in  binary)  has  been  written.  Note  that 
both  "CONVERT"  and  "NICL"  have  to  be  compiled  in  Fortran  5.  The  input  file 
(tape2)  of  "CONVERT"  should  be  "filename".  The  output  file  (tape3)  of 
"CONVERT”  is  "PEXP".  It  contains  the  waveform  transferred  from  the  Nicolet, 
in  the  following  format: 

WRITE(3,100)  Tlf1E(l),VOLT(l) 

IOO  F0RMAT(2(5X,E  1 2.4)) 

where  TIME(I)  and  VOLT(I)  are  the  coordinates  of  the  i-th  data  point.  A 
listing  of  "NICL"  and  "CONVERT"  is  given  below. 


PROGRAM  NICL1TAPE2.  PSAVE.  TAPE3-PSAVS) 

C 

C  THIS  PROM  AH  ALLOWS  THE  TRANSFER  OF  DATA  FROM  A  NICOLET  DIGITAL 
C  OSC  ILL  I  SCOPE  4094  WITH  ONE  DiBK  DRIVE  TO  A  CVBER. 

C 

C  A  PROCEDURE  IS  CREATED  WHICH  WILL  SAVE  ALL  LOCAL  WAVEFORH  FILES  AS 
C  PERMANENT  WAVEFORM  FILES  WHEN  IT  IS  EXECUTED. 

C 

C  BEE  NOS  S.  I.  VOL.  4.  CH.  12  FOR  AN  EXPLANATION  OF  TRANSPARENT  OUTPUT 
C  1:0)  AND  TERMINAL  REDEFINITION  (HALF!. 

C 

LOGICAL  CONT 

CHARACTER  CHD#2,  0PER*4 

WRITE <3.  '!*.  PROC.  SPROC.  *./.  *CTIHE.  *)') 

C 

C  INITIATE  HALF  DUPLEX  TO  PREVENT  T)C  CVBER  FROM  ECHOING  THE  NICOLET  *S 
C  OUTPUT. 

C 

HALF  -  0*00144041400000000000* 

WRITE!#,  '!A10>'>  HALF 
C 

C  TURN  THE  NICOLET  ON  WITH  A  CTL-A.  CHANGE  THE  COMMAND  DELILMITER  TO  A  ' 
C  AND  CHANGE  THE  DEFAULT  SPACING  OF  DATA. 

C 

WRITE!#.  *  I  *:  GOA*.  *>C.  13.33l*>'» 

READ  #.  LERR1 

WRITE!#.  '!*C.4.  1.32!*)') 

REAS  #.LERR2 
WRITE!#.  *!*«!*»*» 

READ  LERR3 

WRITE!#.  33)  LERR1.LERR2.  LERR3 
38  FORMAT ! 'ERROR  CODES  *.313) 

CONT  -  .  TRUE. 

30  WRITE  I#.  '! ‘ENTER  COHtMND*  > ' ) 

READ!#.  13)  CHD 
IS  FORMAT!  A2) 

IF  !CMD.  EG.  'DM')  THEN 
OPER  -  'R.  O.  ' 

CALL  DI8HM!0PER.  HEH8EC) 

ELSE  IF  !CMD.  EG.  'DC')  THEN 
CALL  DI8HC 

ELSE  IF  ICMD.  EG.  'HD')  THEN 
OPER  -  '8.  O.  ' 

CALL  D1BHMI0PER.  HEHSEC) 

ELSE  IF  !CMD.  EG.  *HC'>  THEN 

'WRITE!#.  'CENTER  WAVEFORH  NUMBER*)') 

READ!#.  33)  NUAVEN 
33  FORMAT!  12) 

CALL  HEHCINWAVEN) 

ELSE  IF  <CHD.  EG.  'WF')  THEN 
CALL  WAVEF 

ELSE  IF  !CMO.  EG.  <0U')  THEN 
CONT  «  .FALSE. 

ELSE 

WRITE!*. 'I -INVALID  COMMAND’) ') 

END  IF 

IF  ICQNT )  GOTO  30 

WRITEI3.  'I'REVERT.  PSAVE.  *)') 

END 

C 

C  DISH  TO  MEMORY  OR  MEMORY  TO  DISK  TRANSFER 
C 

SUBROUTINE  DISKMtOPER.  HEHSEC) 

CHARACTER  0PER*4 

WRITE!#.  '! -RECORD  NUMBER*)') 

READ  #.  NREC 

WRITE!#.  '! “MEMORY  SECTION*)') 

READ  #.  HEHSEC 

WRITE!#, '<•: G3A>  *)') 

WRITE!#.  40)  OPER,  NREC.  HEHSEC 
40  FORMAT! A4.  12,  *,  *.  12.  *!•) 

READ  #, LERR 
WRITE!#. ' ! *0 ! • ) ' ) 

READ  #.  LERR2 
WRITE!#,  13)  LERR.LERR2 
13  FORMAT ! -ERROR  CODES  ',213) 

RETURN 

END 
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r.r. 


ICHORV  TO  CYBER  TRANSFER 

SUBROUTINE  HEHC  (NWAVEN) 

DIMENSION  I  DATA  <  16000).  IWDAT  (39.  7) 
CHARACTER  ruHM,  F*22 

REAO  PARAMETERS  PROM  USER 

WRITE!*.  '(  "START  I  NO  POINT’)') 

REAS  *.  I START 

WRITE!*.  'I ’TOTAL  NUMBER  OF  POINTS*)') 
Ae AD  *.  I TOTAL 
WRITE!*.  '(’FILENAME*)') 

READ(*.  MA4) ')  FLNH 
WRITE!*.  '(*:09A>  •>') 


READ  SATA  ABOUT  NAVES  FROM  THE  NICOLET. 

WRITE!*.  '(’C.4.0!*)') 

READ  *.  LCS1 

WRITE!*.  '(’C.O.  1. 13!’)*) 

READ  *.LCS2 
WRITE!*.  '<’H!*)'> 

READ  *.  LUR1 

READ!*.  MIS)')  ITOT 

ITEM  -  I  TOT/10  *  I  CHAR  (  '0*  ) 

I0NE  -  MODI  ITOT.  10)  *  ICHAR('O') 

F  -  '(  '//CHARI  I  TEN) //CHAR!  IONE)//  '(12. 16.  212.  311  )*  13)  ' 

READ!*. F)  ((IWDAT(1.K>.K*1.7>.  I-I.  IT0T).LNR2 
WRITE**.  '(’C.  0. 0!  ’)  ' ) 

READ  •*  LC83 

WRITE!*.  '(’C.4.1.32!*)') 

REAO  *.LC84 
I  STEP  -  1 

READ  DATA  FROM  NICOLET  TO  CYBER  SO  VALUES  AT  A  TIME  (ISO  CHARACTERS). 
THE  CYBER'S  SURFER  IS  ONLY  190  CHARACTERS  LONO  AND  IT  NEEDS  TIME  TO 
PROCESS  DATA  FROM  THAT  BUFFER. 

ITL  -  20 

1ST  -  I  START 

ILP  -  I TOTAL  /  ITL 

IREM  -  I TOTAL  -  ILP* ITL 

DO  39  1*1.  ILP 

WRITE!*.  29)  NWAVEN.  1ST.  ITL.  ISTEP 
9  FORMAT! *0.0.  *.  IS.  *.  *.  16.  *.  *.  16.  ’,  ’.  13.  •!■) 

REAO  *•  LERRI 

READ  *.  (1DATA(K).K-IBT*1.  IST*ITL).  LERRS 
»  I8T  -  1ST  ♦  ITL 

IF  (IREM.  EO.  0)  00  TO  49 
WRITE!*. 29)  NWAVEN. 1ST.  IREM. ISTEP 
READ  *.  LERR 

READ  *.  (IDATA(I). I-IST+I, I8T*IREH).LERR 

READ  NORMALIZATION  DATA 

I  WHITE!*.  30)  IWDAT INHAVEN.  3) 

I  FORMAT! ’N.  *,  IS.  *!’) 

READ  *.  LERR3 

READ  *.  IVAL.  I  SWEEP*  ICHAM..  IDI8P.  VNORM.  HNORM.  IVZERO.  IHZH.  IHZL. 

1  IRNVZ.  IRNHZ.  LERR 4 

WRITE!*. '(*0!*)') 

READ  *. LERR 9 

WRITE!*.  90)  LWR1.LNR2. LERRI. LERRS. LERR3. LERR4.LERR9 
FORMAT! ’ERROR  CODES  ’.713) 

WRITE  ALL  DATA  OUT  TO  FLNH 

OPEN (2.  FILE  -  FLNH) 

REWINDIS) 

WRITEIS.  *)  NWAVEN.  ISTART.  I  TOTAL.  ISTEP,  IWDAT  I  NWAVEN.  3). 

1  (1DATA!!).  I-IBTART*!,  ISTART*  I  TOTAL).  IVAL.  ISWEEP.  ICHANL.  IDISP. 

8  VNORM.  HNORM.  IVZERO.  IHZH.  IHZL.  IRNVZ.  IRNHZ 
CLOSE! 2.  STATUS  *  'KEEP  ' ) 

WRITE  1 3. 40)  FLNH 
FORMAT  (’REPLACE,  ’.A.  •) 

RETURN 

END 


noon  on** 


DISK  TO  CYBER  TRANSFER 

SUBROUTINE  DI6KC 
DIMENSION  IWDATI33.  7) 

LOGICAL  FOUND 
CHARACTER  0PER*4.  F*22 

READ  HAVE  FROM  DISK  TO  NICOLET 

OPER  •  'R»  O.  * 

CALL  DISKH(OPER.MEHSEC) 

AFTER  THE  WAVE  HAS  BEEN  READ  IN.  ITS  WAVEFORM  NUMBER  HAS  TO  BE  FOUND 
80  ITS  DAT AROINTS  CAN  BE  TRANSFERRED  TO  THE  CYBER. 

00  TO  ! 10.  20.  30.  30. 10.  20) .  MEHSEC 
10  NWAVEN  -  1 

00  TO  70 

20  NWAVEN  -  2 

00  TO  70 

30  WRITE!*.  M“:®SA>  •)'» 

WR 1TE < *._' !^C.  4. 0!  “)  #) 

Read  *.  lcsi 

WRITE!*.  '!“C.O.  1. 13!“)  '1 
READ  *.  LCS2 
WRITE!*.  M“H!“)') 

READ  *.  LWRI 

READ!*.  '! 12)')  IT0T 

ITEN  -  ITQT/10  ♦  1CHAR! '0') 

I ONE  -  MOD < I TOT.  10)  *  ICHAR('O') 

F  -  'I '//CHAR! ITEN) //CHAR HONE)// '112. I*.  212.  311).  12)' 

READ!*.  F)  ! ! IUDAT ! I.  R)»  K“l.  7). 1-1.  IT0T).LWR2 
WRITE!*.  '!“C.  0.0!“) ') 

READ  *.  LCS3 

WRITE!*.  '!“C.4.  1.32!“)') 

READ  *.  LCS4 
WRITE!*. ' ! “Q! “) ' ) 

READ  *.  LERR 

WRITE!*.  90)  LWRJ. LWR2,  LERR 
SO  FORMAT! “ERROR  CODE  “.313) 

FOUND  *  .FALSE. 

DO  33  1-1.  ITOT 

1F!IWDAT!!<  7).  EO.  MEHSEC)  THEN 
NWAVEN  >  IWDAT ! I. 1 ) 

FOUND  ■  .TRUE. 

END  IF 

33  CONTINUE 

IF  !.  NOT.  FOUND)  00  TO  60 

70  CALL  HEHC! NWAVEN) 

RETURN 

80  WRITE!*.  '! “INVALID  WAVEFORM  NUMBER")') 

RETURN 


> 


c 

C  DISPLAY  DATA  ABOUT  WAVEFORMS  IN  MEMORY 
C 

SUBROUTINE  WAVEF 
CHARACTER  F*2 2 
DIMENSION  IWDAT <35. 7) 

WRITE!*. '<":09A>  ")') 

WRITE I*. ' 1 "C.  4. 0! ") ' ) 

READ  *.  LC8I 

WRITE!*.  ' < "C.  0.  I.  I3!*)M 
READ  •«  LC82 
WRITE!*.  '!*W!") ') 

READ  *.  LWR 1 
READl*j_l<.I2>'>  ITOT 
TtEN  »  1T0T/10  ♦  I CHAR! *0*1 
I ONE  -  MOD! ITOT,  JO)  *  J CHARI '0') 

F  m  '(  '//CHAR ! ITEM) //CHAR HONE)//'! 12.  16.  212.  311).  12)' 

READ  I*.  Fill  IWDAT  !  I.  K).  K“l.  7).  1-1.  IT0TI.LWR2 
WRITE!*.  ' ! "C.  O.  0!  " )  ' ) 

READ  •<  LCS3 

WRITE!*.  '!"C.4.  1.  32!") ') 

READ  *.  LCS4 
WRITE!*.  '! "O! ")  ') 

READ  *.  LERR 

WRITE!*.  IS)  LWR1 ,  LWR2, LERR 

19  FORMAT ! "ERROR  CODE  *.313) 

WRITE!*.  20)  ITOT 

20  FORMAT! 'TOTAL  NUMBER  OF  WAVES  '.  13. /> 

WRITE!*.  30) 

30  FORMAT! 'WAVE  NO.  OF  POINTS  NORM.  SET  NORM.  STEP  CHANNEL 

1REFERENCE  TITLE') 

WRITE!*.  40)  ! ! IWDAT !I.  K). K*l«  7).  1-1.  ITOT) 

40  FORMAT! IX.  12.  AX.  IA.  9X.  12.  IOX.  12.  9X.  I1.9X.I1.BX.  ID 

RETURN 
END 

11.41.96.  UCLP.  AA.  PTRL09. 


O. 282KLNS. 


**  END  OF  LISTINO  •• 


PROOrfta  CONVERT (PNIC.  TAPE2-PNIC. PEXP. TAPE3-PEXKI 


e 

C  THIS  PROGRAM  IS  USED  WHEN  TRANSFERRING  DATA  FROM  THE 

:  NICOLET  TO  THE  CYBER  .  IE  AFTER  USING  PROGRAM  NICOLT. 

E  THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  READ  THE  OUTPUT  FILE 

E  GENERATED  BY  NICOLT  IN  FREE  FORMAT  (  HERE  IT  IS  CALLED  PNIC) 

E  AND  TO  CONVERT  IT  TO  REAL  VALUES  IN  A  DE8IRED  FORMAT. 

C  THE  OUTPUT  FILE  OF  THIS  PROGRAM  18  CALLED  HERE  PEXP 

E 

E  THIS  PROGRAM  HA8  TO  BE  WRITTEN  IN  FTN9 

E 

E  NOTATION: 

E  TIME ( I )-  REAL  TIME 

E  VOLT <  I  >«  REAL  AMPLITUDE  OF  THE  SIGNAL 

E  RTIME( I >-  RESET  TIME  (IE  WHEN  USING  THE  RESET  ON  NICOLET) 

E  RVOLTCI >-  RESET  VOLTAGE  AMPLITUDE  (IE  IDEM  ) 

E 

E  DIMENSION  I  SHOULD  BE  GREATER  THAN  NF 

C  NF^NOtN 

E  NO-  NUMBER  OF  THE  FIRST  DATA  POINT  TRANSFERRED  (SEE  NICOLT) 

E  N-TOTAL  NUMBER  OF  DATA  POINTS  TO  BE  TRANSFERRED. 

C 

DIMENSION  TIME <9000).  VOLT (3000), RTIHEiSOOO).  R VOLT (9000) 
DIMENSION  1 DATA (9000) 

E 

E  READ  OUTPUTFILE  (TAPED)  FROM  NICOLT 

E 

READ (2. •)  NUAVEN. ISTART. I TOTAL. I STEP.  NR SET, 

KIDATA(I).  I-ISTART. IBTART+ITOTAL-1 ) , IVAL.  I SWEEP. 

2ICHANL.  1DISP.  VNORM.  HNORM.  IVZERO.  IHZH.  IHZL. 

3IRNVZ.  IRNHZ 

B 

E  CONVERT  TO  REAL  VALUES 

E 

DO  1  I-ISTART.  ISTART+ITOTAL-1 
TIME  <  I  >-(  ( I-I )— (65936+IHZH+IHZL)  )+HNORM 
VOLT ( I )•( I DATA ( I )-IVZERO)+VNORM 
RTIME( I )-( ( 1-1 ) -IRNHZ )+HNORM 
RVOLTd  )  — <  IDATA(  I  -IRNVZ)+VNORM 
1  CONTINUE 

C 

C  WRITE  IN  FORMAT  100  THE  OUTPUT  FILE  (TAPE3)  (HERE  PEXP) 

C  CONTAINING  ONLY  TIME(I)  AND  VOLT(I) 

C 

DO  2  I-ISTART.  ISTART+ITOTAL-1 
MR  I  ft' <3.  1 00)  TIME  ( I).  VOLT  ( I ) 

S  CONTINUE 

100  FORMAT (9X.E12.  4.  9X.EJ2.  4) 

STOP 

END 

10.  43.  07.  UCLPt  AA.  PTRL05.  0.  066HLN8.  **  END  OF  LISTING  ** 
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APPENDIX  D 


THE  INHOMOGENEOUS  VISCOUS  WAVE  EQUATION  FOR 
LASER-INDUCED  SOUND 

The  wave  equation  describing  the  pressure  field  of  a  lossless  medium 
containing  a  heat  source  was  first  derived  by  Ingard7  Starting  from  this 
linear  inviscid  thermoacoustic  wave  equation,  Westervelt  and  Larson 
studied8  the  special  C8se  of  a  sound  field  generated  by  the  heating  of  a 
laser  beam  penetrating  into  a  liquid,  end  they  discovered  the  fescineting 
properties  of  whet  is  now  called  a  thermoecoustic  array.  In  this  appendix, 
we  rederive  from  the  hydrodynamicel  equations  of  motion  the  so-called 
Westervelt-Larson  wave  equation  for  laser-induced  sound,  including  now 
the  effect  of  viscosity. 

Let  p=p0+p  be  the  total  pressure  in  the  medium,  where  p0  is  the 
ambient  pressure  end  p  is  the  acoustic  pressure  generated  through  the 
thermal  mechanism.  Similarly,  let  p=p0+p  be  the  total  density  of  the 
medium,  T=T0+T  be  the  total  temperature,  and  s=s0+s  be  the  total  entropy 
in  the  fluid  medium.  The  subscript  o  refers  to  an  undisturbed  quantity. 

Since  the  acoustic  disturbances  generated  by  the  thermal  mechanism 
are  expected  to  be  small,  we  start  directly  from  the  linearized 
hydrodynamicel  equations  of  motion,  the  conservation  of  mass,  and  the 
conservation  of  linear  momentum. 


5 


mess: 


pt  +  p0  V.  u  =  ° 


(D1) 


momentum: 


p0  +  Vp  =  (\+2p)  V2u 


(D2) 


Here  u  denotes  the  particle  velocity,  (X*2y)  represents  the  viscosity  of  the 
medium,  end  the  subscript  t  denotes  a  time  derivative.  It  is  implicitly 
assumed  in  Eqs.  (Dl)  8nd  (D2)  that  the  medium  is  at  rest,  steady,  uniform, 
and  irrotational,  and  that  heat  conduction  can  be  neglected. 

The  thermal  input  from  the  laser  intensity  induces  a  change  in  entropy 
in  the  system,  so  that  an  adequate  equation  of  state  describing  the  process 
is 

st8te.  p  =  p  (  p  ,  s  )  (D3) 

Forming  a  Taylor  series  expansion  of  Eq.  (D3)  about  the  equilibrium 
quantities  p0,  p0,  and  s0,  we  have 

p  =  p0  +  (3p/ap)s  (p-po)  +  (ap/ds)p  (s-so)  +  ...  (D4) 

By  definition  the  smell  signal  sound  speed  is  c  =  V  (dp/dp)s  ,  so  that 
Eq.  (D4)  can  be  rewritten  as 

p  =  p/c2  +  s  (dp/ds)p  (D5) 

The  partial  derivative  in  Eq.  (05)  is  the  result  of  the  change  in  density  due 


to  the  change  in  entropy  induced  by  the  laser  heating.  It  is  therefore  related 
to  the  coefficient  of  thermal  expansion  p  of  the  medium.  By  definition45 

p  -  -  (1/p)  (dp/*T)p  (D6) 

or 

p  =  -  (1/p)  (dp/ds)p  (ds/dT)p  .  (D7) 

But  it  is  well  known  from  thermodynamics46  that  the  specific  heat  at 
constant  pressure  cp  can  be  expressed  as 

cp  =  T  (ds/dT)p  (D8) 

Combining  Eq.  (D8)  with  (D7),  and  linearizing  the  result  about  the 

undisturbed  values,  yields 

(dp/ds)p=  -p0pT0/cp  (D9) 

Therefore  Eq.  (D5)  becomes 

P =  p/c2  -  sp0pT0/cp  ,  (DIO) 

and  the  last  term  in  Eq.  (DIO)  represents  the  effect  of  the  laser  heating  on 
the  change  of  density  in  the  medium.  Since  we  seek  a  wave  equation  for  the 
acoustic  pressure  p,  we  eliminate  u  between  Eqs.  (Dl)  and  (D2).  It  is  found 


V2P  -  Ptt  +  (*+2nVPo  ^2Pt  =  0 


(Dll) 


Eliminating  9  between  Eq.(QIO)  and  Eq.  (Dll)  yields  a  wave  equation  in  terms 
of  the  acoustic  pressure  p  and  the  input  entropy  disturbance  s, 

V2p  -  (c“2)ptl-  (p0pT0/cp)su  + 


(X+2p)/p0  V2  [  (c“2)pt-  (p0pT0/cp)  st)  =  0 


(D12) 


Now  let  us  relate  the  change  in  entropy  to  the  heat  q  added  to  the 
medium  by  the  laser  beam.  The  conservation  law  which  governs  that 
relation  is  the  second  law  of  thermodynamics  for  a  reversible  process. 


p  T  ds  *  q  dt 


(D13) 


or,  in  linearized  form, 


Po  T"o  st  =  Q 


(D14) 


Combining  Eq.  (D14)  with  Eq.  (D12)  to  eliminate  the  entropy  fluctuation  s,  we 


obtain 


V2p  -  c"2pu  +  (X+2y)/p0c2  V2pt  = 


-  (p/cp)  qt  +  (p/Cp)  (X+2y)/p0  V2q 


(D15) 


This  is  the  Westervelt-L8rson  wave  equation  for  a  viscous  (but  non- 
heet-conducting)  irrotational  fluid  medium  containing  a  heat  source  of 
strength  q.  For  a  lossless  medium,  Eq.  (D15)  reduces  to  its  standard  form 

□2P  =  V2p  -  c“2  pu  =  -(p/Cp)  qt  ,  (D16) 

where  O2  denotes  the  O'Alembertian  operator. 

It  is  also  convenient  to  relate  the  energy  q  added  to  the  medium  to 
the  intensity  of  the  laser  beam  illuminating  the  medium.  Let  us  denote  by 
<  I  >  the  intensity  of  the  laser  beam  averaged  over  8  period  of  the  very 
rapid  (  ~  10~14  s  )  optical  oscillation  of  the  laser  light.  <  I  >  is  the 
Poynting  vector  of  the  process,  and  the  flux  of  intensity  V.  <  I  >  represents 
the  heat  absorbed  by  the  medium  per  unit  time  and  unit  volume.  In  other 
words 

q  =  -  V  .  <  I  >  (D1?) 

In  practice  the  laser  beam  is  shining  in,  S8y  the  +z  direction,  so  that 
qs  -  dl/dz  .  Equations  (016)  and  (017)  are  the  starting  point  of  Chapter  III. 

It  is  interesting  to  note  that  Eq.  (D16)  has  the  same  form  as  the 
Westervelt  inhomogeneous  wave  equation4^  describing  parametric  end-fire 
radiation.  In  the  wave  equation  for  the  parametric  array,  however,  the 
source  term  represents  virtual  sources  generated  by  nonlinear  interaction 
in  the  medium,  whereas,  in  Eq.  (016)  ,  the  source  term  represents  the  time 


and  spatial  derivative  of  the  laser  intensity  distribution  in  the  laser  beam. 
There  is  nevertheless  a  close  mathematical  analogy  between  the 
parametric  array  and  the  thermoacoustic  array,  and  this  analogy  has  been 
used  by  Novikov,  Rudenko,  and  Timoshenko4®  to  study  the  pressure 
waveform  radiated  by  a  transonic  thermoacoustic  array. 
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APPENDIX  E 


HALF-ORDER  DERIVATIVE  AND  IMPULSE  RESPONSE  hL(t) 

AND  CONNECTION  WITH  THE  PARABOLIC  CYLINDER  FUNCTION 

The  concept  of  h8lf-order  derivatives  appears  as  early  as  1695,  in  a 
letter  from  Leibniz  to  L'Hospital,  opening  thus  the  8venue  to  wh8t  is  now 
called  fractional  calculus.  A  complete  historical  survey  of  the  ni8in 
contributions  to  fractional  calculus  can  be  found  in  reference  49.  In  this 
appendix,  it  is  shown  that  the  half-derivative  operator  is  closely  related  to 
the  vertical  impulse  response  hj_(t)  defined  in  Chapter  III.  It  is  also  shown 
th8t  the  h8lf-order  derivative  C8n  be  expressed  in  terms  of  a  parabolic 
cylinder  function. 

1.  Half-order  derivatives 

It  has  been  shown50  that  the  acoustic  pressure  response  p(t)  of  an 
infinite  line  source  producing  a  uniform  mass  outflow  l(t)/h  per  unit 
length  h  is  given  by 

co 

p(t)  =  (2irh)“ '  lt(t-r/c)  [r2-s2l-1/2  dr  ,  (El) 

s 

where  s  =  r0  sinfl0  (see  Fig.  9),  and  where  the  subscript  1  denotes  8  time 
derivative.  If  it  is  assumed  that  most  of  the  contributions  come  from  the 
point  of  the  line  source  right  in  front  of  the  receiver  (i.e.,  from  r=s),  then 


r2-s2  s  2s(r-s)  ,  and  Eq.  (El)  can  be  rewritten  in  the  form 

t-s/c 

p(t)  =  (2irh)_1  (c/2s)1/2  J  lt(T)  l  (t-s/c)-T]"1/2  dT  (E2) 

-oo 

Using  the  identity50 

t 

(d/dt)1/2  l(t)  =  j  tt(T)  [  TT(t-T)  ]",/2  dT  ,  (E3) 

~  CO 

the  acoustic  pressure  p  can  be  expressed  in  terms  of  the  half-order 
derivative  (d/dt)1/2  of  the  source  function  I,  as  follows. 

p(t)  =  (2h)_1  (c/2irs)1/2  (d/dt)1/2  l(t-s/c)  (E4) 

The  notation  (d/dt)^2  for  the  integral  defined  by  Eq.  (E3)  comes  from  the 
f8ct  that 

(d/dt),/2  [  (d/dt),/2  (1)]=  dl/dt  (E5) 

Note  from  Eq.  (E3)  that  (d/dt)^2  !(t)  can  be  expressed  as  a  convolution 
between  lt(t)  and  a  half-order  delta  function  81/>2(t)  or  "h8lf-order 
derivative  operator". 


p(t)  =  81/2  (t-s/c)  *  lt(t) 


(E6) 


But  we  recall  from  Chapter  111  that  the  impulse  response  hL(t)  is  defined  as 

p(t)  =  hL(t)  *  lt(t)  (E7) 

Equations  (E6)  and  (E7)  show  that  there  is  a  correspondence  between  the 
impulse  response  h^t)  and  the  half-order  derivative  operator.  It  should  be 
emphasized  however,  that  the  infinite  line  source  model  used  to  introduce 
the  concept  of  a  half-order  derivative  in  Eq.  (El)  differs  in  two  ways  from 
the  physical  problem  of  8  thermoacoustic  source  as  discussed  in  this  study. 
First,  we  are  dealing  in  our  experiment  with  a  thermoacoustic  source  over 
a  pressure  release  boundary,  and  not  with  an  infinite  line  source.  Second, 
we  are  dealing  with  8n  exponentially  shaded  line  source,  and  not  with  a 
uniform  line  source.  This  means  that  the  analogy  between  the  impulse 
response  hL(t)  and  the  half-order  derivative  operator  is  strictly  applicable 
only  for  the  case  of  very  weekly  absorbing  media.5^ 

2.  Connection  between  the  half-order  derivative  and  the  parabolic  cylinder 
function 

An  equation  which  occurs  often  (see  for  instance  reference  52)  in 
describing  the  pressure  field  of  a  cylindrical  line  source  of  Gaussian  cross 
section  distribution  is 


(E8) 


F(g)  =  J  6  uc  exp  [ -(d-gr  ]  dd 
o 

It  cen  be  shown  that  in  fact  F(g)  is  related  to  both  the  half-order 
derivative  of  a  Gaussian  function,  and  to  the  parabolic  cylinder  function. 
This  provides  a  connection  between  these  two  functions. 

Combining  Eqs.  (E3)  and  (E8),  we  obtain 

F(g)  ■  (ir/2)  (d/dg)~1//2  exp(-g2)  (E9) 

Using  the  identity53 

CO 

D_1/2(z)  =  ir"1/2  exp(-z2/4)  |  expl-zx-x2/2]  x~1/2  dx  ,  (ElO) 

o 

where  D_1/2  is  the  parabolic  cylinder  function  of  order  -1/2  and  combining 
Eqs.  (E8)  and  (ElO)  yields 

F(g)  -  (2)",/4  tt^2  exp(-g2/2)  D.t^  (-g/2)  (Ell) 

Combining  Eqs.  (E9)  end  (Ell)  gives  an  expression  which  relates  the 
half-order  derivative  of  a  Gaussian  to  the  parabolic  cylinder  function  D_1/2  . 


APPENDIX  F 


DETERMINATION  OF  THE  OPTICAL  COEFFICIENT  OF  ABSORPTION 

The  optical  coefficient  of  absorption  oc  dictates  the  penetration 
depth  of  the  laser  beam  in  the  water,  according  to  8n  exponential  law  of 
attenuation 


l(z)  =  l0e"cCz  ,  (FI) 

where  I  is  the  laser  intensity  of  the  laser  beam,  8nd  z  is  the  vertical  axis 
of  penetration  of  the  laser  beam.  Therefore  1/cC  is  a  measure  of  the  1/e 
length  of  the  thermoacoustic  source.  This  means  that  it  is  necessary  to 
determine  a:  for  both  optical  wavelengths  used  in  our  experiment:  1.06  ym 
(Nd-.Glass)  and  0.6943  ym  (ruby). 

This  was  done  in  reference  14  for  the  case  of  absorption  of  light  in 
distilled  water  and  fresh  lakewater.  In  our  study,  the  same  experiment  was 
repeeted  for  the  case  of  absorption  of  light  in  the  fresh  water  of  the  tank 
in  which  the  experiments  were  made.  The  details  of  the  experimental 
procedure  used  to  measure  a  can  be  found  in  reference  14.  The  method  is 
based  on  measuring  the  attenuation  of  light  for  different  widths  of  the 
sample  of  water.  The  value  of  oc,  in  nepers  per  meter,  is  obtained  by 
measuring  the  slope  of  the  attenuation  versus  distance  on  a  logarithmic 
scale. 

Figure  41(a)  shows  that  the  absorption  of  light  at  an  optical 
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wavelength  of  1.06  ym  is  in  the  water  tank  about  13.7  Np/m.  Similarly, 
Fig.  41(b)  shows  that  c£  is  about  1.5  Np/m  for  attenuation  of  ruby  light  in 
the  same  water.  Note  that  the  numerical  values  of  a  can  differ 
significantly  from  one  study  to  another.10^4'15'54  This  variation  is 
attributed  to  the  fact  that  a  is  highly  dependent  on  the  optical  properties 
of  the  medium  being  used. 


APPENDIX  6 


NUMERICAL  PROGRAM  "MTS" 

A  computer  program  called  "MTS-  has  been  written  in  Fortran,  in  order 
to  compute  the  pressure  radiated  by  a  moving  thermoacoustic  source.  The 
numerical  code  is  based  on  the  analysis  presented  in  Section  IV.A. 

The  input  data  to  the  program  are  variables  such  as  the  geometry 
between  the  source  and  the  receiver,  the  optical  coefficient  of  absorption, 
the  source  velocity,  the  modulation  frequency  of  the  laser  intensity,  the 
laser  pulse  duration  and  the  optical  energy  going  into  the  water,  the  laser 
beam  radius  on  the  water  surface  and  the  intensity  distribution  (either 
Gaussian  or  uniform)  across  a  section  of  the  laser  beam.  The  output  of  the 
program  is  a  plot  of  the  pressure  waveform  at  the  receiver,  with  its 
corresponding  acoustic  level  in  decibels,  from  either  peak  to  peak  or  rms 
pressure. 

A  listing  of  the  program  "MTS”  is  given  below. 


PROCNAM  MTS  (OUT.  INPUT.  OUTPUT.  TAPES* OUT.  RESULT.  TAPE3^RE8ULT. 
SPLOT.  TAPE  5-PLOT) 

DIMENSION  TIME <5000).  HI  100).  P(SOOO).  PT<5000) 

DIMENSION  PULSE <5000).  T1HE2(5000).  P2<5000).  PLSDOT(SOOO) 

EXTERNAL  PLS 
COMMON  DT.  RF.  TT.  T.  XIO 
DATA  LABX/*TIM£<MS)*/ 

DATA  LAB V/ "PRESSURE*/ 

LUNP-6L OUTPUT 
C 
c 

C******************************************************************** 

C  THIS  PROGRAM  COMPUTES  THE  PRESSURE  WAVEFORM  RADIATED  BY  A 

C  MOV I NO  THERMOACOUSTIC  ARRAY. 

C  (SEE  NOTES  5/18/83)  3-D  PROBLEM 

C 

C *********** ********************************************************* 

C  REVISED  02/07/85 

C 

C 

C  PLS  DEFINES  LASER  PULSE)  STORED  IN  PULSE 

C  MAX  NUMBER  OF  POINTS- 500 

C  TIME  DERIVATIVE  OF  LASER  PULSE  IS  STORED  IN  PLSDOT 

C  P(U>— PRESSURE  RADIATED  BY  ARRAY  IN  POSITION  I  WHEN 

C  EXCITED  BY  AN  IMPULSE  PLSDOT(I) 

C  P<J>  IS  COMPUTED  INSIDE  A  LOOP  ON  THE  POSITION  I  . 

C  AND  EACH  CONTRIBUTION  IS  PROORESSIVELY  ADDED  IN  AN 

C  ARRAY  PT<J>. 

C 

C  J  REPRESENTS  THE  (HORIZONTAL)  TIME.  . 

C  PT (J)— TOTAL  PRESSURE  RADIATED  BY  THE  MOV I NO  SOURCE. 

C 

C  N- NUMBER  OF  POINTS  REPRE8ENT1NQ  PULSE 

C  TT-TOTAL  TIME  DURATION  OF  THE  PULSE 

C  DT-TT/N-TIME  RESOLUTION 

C  RF-RADIAN  FREQUENCY  (FACTOR  OF  2  BECAUSE  MODULATION  18  A  SIN**2) 

C***** ****************************************************** ********** 

C 

C 

C  INPUT  VARIABLES 

C 

PRINT*.  "****  WELCOME  AND  BONJOUR  •*«•" 

.PRINT*.  "  * 

PRINT*.  •  PRESSURE  RADIATED  BY  A  H.  T.  A.  * 

PRINT*.  *  * 

PRINT*.  "ENTER  THE  FOLLOW I NO  PARAMETERS" 

C 

PRINT*.  "INITIAL  RANQE  R  <M)  -• 

READ*.  R 

PRINT*.  "INITIAL  THETA  (DEO)-* 

READ*.  TETAD 

PRINT*.  "INITIAL  PHI  (DEO)  -" 

READ*.  PHID 
PI-3.  141592654 
TETAR«TETAD*P 1/180. 

PH1R— PHI  D*P  I  / 1 80. 

PRINT*.  "ABSORPTION  </M)-* 

READ*.  ALFA 

PRINT*.  "MACH  NUMBER-" 

READ*.  XM 

PRINT*.  "MODULATION  FREQUENCY  (KHZ)" 

READ*. F 
F-F*10O0. 

PRINT*.  "PULSE  DURATION  (MS)-" 

READ*.  TT 
TT-TT/1000. 

PRINT*.  "AMPLITUDE  AT  PHOTODETECTOR  <V>«" 

READ*.  APHID 

PRINT*.  "BEAM  RADIUS  (CM)-" 

READ*.  ACM 
A— ACM*0.  01 
C— 1486. 


ESTIMATE  THE  IMPORTANCE  OF  FINITE  BEAMWIDTH  EFFECTS. 


PRINT*.  *DO  VOU  WANT  TO  INCLUDE  FINITE  BEAMWIDTH  EFFECTS?" 
PRINT*."  YES-1  NO— O  (DEFAULT-1)" 

READ*.  IB 

IF ( IB.  EQ.  0)60  TO  1 
DT-A/(5.  *C) 

IF<IB.  Ed.  DPRINT*.  "TIME  RESOL.  DT  SHOULD  BE  LESS  THAN  ”.  DT 
IF  BEAMWIDTH  EFFECTS  IMPORTANT.  NEED  10  POINT/CROSS  SECTION 
PRINT*.  "TIME  RESOLUTION  (MICROSEC)  -" 

READ*.  DT 
DT-DT*1.  E-OA 
IF (F.  EQ.  0.  >60  TO  2 
DTLIM— 1.  /O.  *F) 

IF  THERE  18  LESS  THAN  S  POINTS  PER  MODULATION  PERIOD 

A  WARN INS  16  PRINTED 

IF. Dr:  LE.  DTLIM)60  TO  2 

PRINT*.  "***WARNIN6  -  DT  18  TOO  LAR6E" 

60  TO  1 


RF-2.  *PI*F 
V— XM*C 

EO-1.  03* APHID 
SO— P I *A*A 
PO-EO/TT 
ATRANS-O.  98 
XIO-PO/SO 
BETA-2.  07*1.  E-04 
CP— 4190. 

XKC-ATRANS*B£TA*ALFA*C/(4. *PI*CP) 

IF ( IB.  EO.  0) XKC-XKC*SO 

IFLAO-O 


COMPUTE  RN.  TETAN  (FINAL  RAN6E  AND  ANCLE  AT  END  OF  PULSE) 

DT-TIME  RESOLUTION 

N-NUMBER  OF  POINTS  IN  LASER  PULSE 

N2— NUMBER  OF  POINTS  TO  BE  STORED  FOR  ANY  P(J> 

SEE  NOTES  FOR  THE  NOTATION 


RN-SQRT ( ( R*R  > ♦ ( V*TT ) **2- ( 2.  *V*TT*R*SIN( TETAR >*COS(PHIR ) ) ) 
TETAN-ACOS ( R*COS ( TETAR ) /RN ) 

N-INT (TT/DT) 

IF( XM-1 .  )6> 6.  7 
AP-R *S I N ( TETAR ) /C 
EPS-1. 

60  TO  8 

A“- ( RN*SIN ( TETAN ) /C ) *TT 
EPS— 1. 

PRINT*.  "TIME  WINDOW  FOR  THE  COMPUTATION" 

PRINT*.  "INITIAL  TIME  (MS)-* 

READ*.  TO 

PRINT*.  "FINAL  TIME  (MS)-" 

READ*.  TF 
T0-T0*0.  001 
TF-TF*0.  001 
TF-TF+DT 

IF(TO.  LE.  AP)00  TO  11 

PRINT*.  "INIT.  TIME  SHOULD  BE  LESS  THAN  ",AP 
60  TO  8 

N3- INT  ( ( TF-TO )  /DT ) +1 


NZEHO-INT{ <AP— TO)/DT) 

NZEROP-NZERO+1 

N2-N3-NZER0 

I F(N.  GT.  9000.  OR.  N2.  CT.  9000)  I  FLA 0-1 

N  SHOULD  BE  LESS  THAN  9000  < VERTICAL  NUMBER  OF  POINTS) 

N2  SHOULD  BE  LESS  THAN  9000  (HORIZONTAL  NUMBER  OF  POINTS) 

PRINT#.  “  ■ 

PRINT*.  "RESOLUTION**.  DT 

PRINT*.  "N  (VERTICAL  TIME)  «“.  N 

PRINT*.  "N2  (HORIZONTAL  TIME)  N2 

PRINT*.  "NZERO  -“» NZERO 

ENOFILE  LUNP 

IFdFLAC.  EO.  1)00  TO  999 

PRINT*.  "IS  THAT  OK?  (YES-1.  NO-O) “ 

READ*.  I STOP 

IFdSTOP.  EO.  0)00  TO  999 


PRINT  HEADINGS 
PRINT (2.  100) 

F0RMAT(20X.  "MOVING  THERMOACOUSTIC  ARRAY") 

PRINT(2.  101 )R.  TETAD.  PHID.  ALFA.  XM 
PRINT (2. 102)F< N.  TT.  DT 

FORMAT ( ///.  9X.  "R--.E8.  2.  9X.  "THETA  (DEC)-".  EB.  2.  9X 
8"PHI— “.  EB. 2. 9X> "ABSORPTION-".  EB.  2. 9X. “MACH-".  EB.  2) 
FORMAT <9X,  "MOD.  FREQ.  -".  EB.  2.  9X.  “N-".  14.  9X. 

8“TT»".  EB.  2.  9X>  "DT-".  EB.  2) 

PRINT (2.  103) 

FORMAT (///.9X.  "TIME".  10X.  "PRESSURE") 


COMPUTE  PLSDOT 

DO  3  H-l.N 
X -FLOAT (M) 

IF(M.  OE.  2)00  TO  92 
PLSDOT ( 1 )— 0. 

00  TO  3 

PLSDOT ( M )  —  ( PL8 (X)-PLS(X— 1.  )>/DT 
PULSE  <  M ) — PLS ( X ) 

CONTINUE 


INITIALIZE  ARRAY  PT(U) 

DO  30  J-1.N3 
PT  ( J)-0. 

CONTINUE 


COMPUTE  RI.THETAI  FOR  EACH  POSITION  OF  THE  ARRAY. 
THEN  COMPUTE  Pd.J)  LOOPS 


DO  9  I-l.N 
RI-(R#R)*(V*I#DT)**2 

R 1 -SORT < R I  —  ( 2.  *V*1*DT*R#SIN(TETAR)*C0S(PHIR) ) > 
TETI-ACOS ( R*COS  <  TETAR ) /RI ) 

RCI-R1*C0S(TET1 ) 

R8 1 -R I *S I N ( TET I ) 

THE  IMPULSE  RESPONSE  IS  INDEPENDANT  OF  PHI 
EXCEPT  THROUGH  RI 


C  COMPUTE  THE  IMPULSE  RESPONSE 
C 

KFLAO-O 

T2I-(RSI/C>*(I*DT> 

T1I-AP 

T3I«(RI/CJ*<I*DT) 

T3MI-T3I-DT 
DO  10  J-l.  N2 
T»T1I*(U*DT> 

T1ME(J)-T-(I*DT> 

IF<T.  LT.  T2I  >00  TO  36 
IF((T-T2I>.  OE.  1.  E-OS >00  TO  31 
P(J>-(X  C*EXP (— ALFA*RCI I >*PL8DQT ( I I/RBI 
00  TO  10 
36  P(U>-0. 

00  TO  10 

31  R00T-80f<T(  <  (C*TIME(  J)  >**2>-(RBI*R8I  >  > 

IF  (ROOT.  0£.  1.  E— 06)00  TO  32 

PRINT*.  *****  ROOT-O  ****** 

32  IF(T.  OT.  T3HI100  TO  33 

P  (  J)*2.  *XKC*EXP  ( -ALFA*RC  I  >  *C08H  ( ALFA*R00T  )  /ROOT 
P( J>-P( J)*PL8D0T( I > 

OO  TO  10 

33  KFLAO-MFLAO+1 
IFOu-lAO.  NE.  1100  TO  34 

C  IF  KFLAO-l.T-RI/C.AND  P-P(RI/C> 

P(  J)-XHC*EXP(-2.  *ALFA*RCI )*PLSDOT(I)/RI 
OO  TO  10 

34  P(U>— 2.  *XKC*EXP(-ALFA*ROOT)*BINH(ALFA*RCI)/ROOT 
P  (  J  )  -P  (  J  )  *PLSOOT  (  I ) 

10  CONTINUE 


ADD  BUCCECBSIVELY  IN  PT(J>.  ALL  THE  ELEMENTARY 
CONTRIBUTIONS  P(JI.  AS  I  IS  INCREASED  (IE  AS  THE 
SOURCE  18  NOVI  NO). 

DO  33  U-1.N2 

PT< J>-PT(U>+(P(J)*DT> 

CONTINUE 

MOVE  TO  THE  NEXT  POSITION  OF  THE  ARRAY  (NEXT  II 


ADD  NZERO  ZEROES  AT  THE  BEOINNINO  OF  THE  RECORD 
(IE  FROM  T-TO  TO  T-AP> 


DO  600  KK-1.N3 

TIME(KK)-(TO+( (KK-1 )*DT) >*1000. 
IF (KK.  OT.  NZERO >0O  TO  601 
P2(HH>-0 
00  TO  600 

P2(KK)-PT(KK-NZER0) 

CONTINUE 
DO  606  1-1.  N3 
PT ( I )— P2( I >*EP8 
CONTINUE 


NUftt*. «  0AU6SIAN  Cl)  CM  UNIFORM  <2)* 

READ*.  I SHADE 

CONST-BORTCPI)*C*A*SINCTETAR> 
TB0«CR-CA*8INCTETAR> ))/C 
TBF-CR+CA*8INCTETAR> ) )/C 
1-0 

T-TBO-DT 

T-T+DT 

1-1*1 

IFCT.  OT.  TBFIOO  TO  41 
0-CR-CC*T> )/<A*SINCTETAR)  > 

IF  (  I  SHADE.  EO.  2)H(  1  >  — 1.  *CONST 
IFCIBHADE.  EO.  1  )H(  I  >-CONBT*EXP  <  -0*0) 

00  TO  49 
NPB-I-1 

CONVOLUTION  FOR  FINITE  BEAHMIDTH 
DO  42  I-1.N3 

sun-o. 

DO  43  J-l.NPB 

IF<  J.  OE.  I>00  TO  44 

sun-cPT<i-u)*HCJ> >*suh 

OO  TO  43 

8UH-0. 

CONTINUE 

P2CI)-SUH*DT 

CONTINUE 

DO  46  I-1.N3 
PT<  I )— P2C I ) 

CONTINUE 

PT  INCLUDES  NOW  THE  FINITE  BEAHWIDTH  EFFECTS. 


COMPENSATE  FOR  THE  PHASE  8HIFT  OF  THE  CONVOLUTION 

DO  614  I-1.N3 

IF ( I .  HE  N3I00  TO  619 

PT C  I )— O. 

00  TO  614 
PT ( 1 )— PT  C  1*1 ) 

CONTINUE 


PRINT  THE  RESULT 

DO  602  I-1.N3 
UR  I  TEC  2.  300>TIMEC I ) >  PTC  I  > 
HR  1  TEC  3.  300)  TIME  C  I  ).  PTC  I  > 
CONTINUE 

FORMAT C2C 9X>  El 2.  4>) 


COMPUTE  THE  RMS  PRESSURE  AND  THE  ACOUSTIC  LEVEL 

PRINT*.  "COMPUTATION  OF  THE  ACOUSTIC  LEVEL* 

PRINT*.*  Cl)  PEAK  TO  PEAK  OR  C2>  RMS  ?* 

READ*.  IDB 

IFC1DB.  EO.  1)00  TO  613 

PRINT*.  *D0  YOU  WANT  TO  SPECIFY  A  NEW  TIME  WINDOW  FOR  THE1 
PRINT*.  “COMPUTATION  OF  THE  RMB  PRESSURE?  C YES- 1 .  NO-O ) * 
READ*.  IRMSW 

IFCIRMSW.  EQ.  1)00  TO  611 

TWO-TO 

TWF-TF 

00  TO  612 


PRINT*.  ‘’WINDOW  FOR  THE  COMPUTATION  OF  THE  RMS  PRESSURE" 
PRINT*.  “INITIAL  TIME  <HB>  FOR  THE  WINDOW-" 

READ*.  TWO 
TWO-TWO/ 1000. 

PRINT*.  "FINAL  TIME(MS)  FOR  THE  WINDOW-" 

READ*. TWF 
TWF-TWF/ 1 OOO. 

TW-TWF-TWO 

NPO- I NT ( < TWO-TO I /DT ) *1 
NPF-INT  < <TWF-T0)/DT>*1 
RMS-O. 

DO  603  I-NPO.  NPF 
RMS-RMS+  <  PT  < I >  **2 ) 

CONTINUE 

RMS-SORT  <RMS*DT/TW) 

PRINT*.  "*********  RMS  PRESSURE  -  ".RMS 
PMAX— 100. 

PHIN-IOO. 

DO  39  I-I.N3 

IF<PT ( I ).  LE.  PMIN)PMIN-PT(I> 

IF(^r«I).  6E.  PMAX JPMAX-PT ( I > 

CONTINUE 
PAMP— PMAX—PMIN 
1F(IDB.  EQ.  lJARCMT-PAMP 
IF( IDB.  EO.  2)AR0MT-RMS 
DB-20.  *AL0010(AR6MT*1.  E+06) 

PRINT*.  "*********  ACOUSTIC  LEVEL-".  DB 
PRINT*.  "  " 

PRINT*.  "DO  YOU  WANT  A  PLOT  <YES-1» N0-0>?" 

READ*.  IPLOT 

IFdPLOT.  EO.  0)60  TO  999 


IHAJm MT  <  C  TPF-TPO  )  /DT  > 

DO  993  1*1.  2 WAX 
P2CX)-PT<X*NDL) 

TIME2C I )-TIHEC I*WDL) 

993  CONTINUE 

DO  994  1*1.  I WAX 
PTU)*f2(I) 

TIWECI)-TIHE2<I) 

994  CONTINUE 
XHIN-TP0*10OO. 

XHAX«TPF*1000. 

997  IFCXNQRM.  EQ.  1)  DO  TO  991 

PRINT*.  "YHI^* 

READ*.  YW IN 
PRINT*.  “YMAX-" 

READ*.  YMAX 
DO  TO  992 

991  YHXN— 1.  9 
YNAX-1.  9 

992  XDEL-CXHAX-XHIN>/10. 

YDEL- < YHAX-YH I N ) /B. 

AXL-9. 

AYL-3. 

XTITL-AXL/2. 

YTITL-AYL-K).  9 
DX- < XNAX-XHJN) /AXL 
DY- <  YHAX-YM I N > / AYL 
CALL  PLTLFN ( 4LPLOT ) 

CALL  PLTOROCXORO. YORO) 

CALL  PLTDIMCli.  *6.  9.  1*4.  > 

CALL  PLTAXI8C0.  .  0.  .  AXL.  O.  •  XWIN.  XWAX.  XDEL.  LABX.  -8.  4.  !•-.!> 

CALL  PLTAXIBCO.  .0.  .  AYL.  90.  ,  YHIN.  YMAX.  YDEL.  LABY.  8.  1.-.  1.-.  1) 
CALL  PLTAXX8<0.  .  AYL.  AXL.  O.  .XWIN.  XWAX.  XDEL.  LABX.  0*  O.  1.-.  1) 
CALL  PLTAXI8CAXL.  O.  .  AYL.  90.  .  YWIN.  YMAX.  YDEL.  LABY.  O.  O.  .  1.  .  1  > 
CALL  PLTDATACTJME.  PT.  IMAX.  O.  .  O.  XWIN.  DX.  YHXN.  DY.  O.  ) 

CALL  PLTL I NE ( XT XTL.  YTITL.  14) 

WR1TEC9.  lOOO) 

1000  FORMAT  I "LASER— INDUCED  THERWOACOU6TIC  PRESSURE" ) 

CALL  PLTLINEC XTITL. 8.-.  1) 

WRITE ( 9/  100DR.TETAD.PHID.  XM 

1001  FORMAT C "R  CM)-". FA.  2.  4X. "THETACD)-".  F4.  1. 4X. "PHI  CD)-". FA.  1. 
84X.  “MACH— ,F4.  1) 

F-F/1000. 

TT-TT*1000 

CALL  PLTLINEC XTITL.  -1.  . 1) 

WRXTTC9.  1002) ALFA. F,  TT 

1002  FORMAT <  "AiS.  CNP/M)-",  FA.  2.  4X.  "FREO.  C  KHZ  >— ".  FA.  2.  4X. 
S"TP(HS)-".FA.  2) 

XFCIDB.  EQ  1  )00  TO  990 
TW-JOOO  • I TW-DT > 

CALL  PLTLINEC  XTITL.  -1.  2.-.  10) 

WR1TEC9.  1003JPAMP. DB.  TW 

1003  FORMAT  C  "PAWP  CPA)  - " .  FB.  4.  3X .  "LEVEL  CDS)  -".FA.  2.  3X 
8.  "TAVE  CMS)  -".FA.  3) 

DO  TO  969 

990  CALL  PLTLINEC  XTITL.  -1.  2.-.  10) 

WRITE ( 9.  1004 )PAMP. DB 

1004  FORMAT C “PAWP  CPA)-".  FA.  3.  9X.  "LEVEL  CDB )-"»  FA.  2) 

C 

969  CALL  PLTENDCO.  .  0.  ) 

999  STOP 

END 
C 

C  DEFINE  EXTERNAL  LASER  PULSE 

C 


FUNCTION  PL8CY) 

COMMON  DT. RF. TT. T. XIO 
T-V*DT 

989  PI-3.  141 992 A 94 

A- 10.  8 
B-5.  O 

ENVLP-X I O*  C  A*T /TT ) *E  XP  C  — B*T/TT ) 

PLS-C  <8lNCRF*T/2.  ))**2)*ENVLP 

IFCRF.  EQ.  0.  1PLB-ENVLP 

RETURN 

END 


11.43.  39.  UCLP.  AA.  PTRL09. 


0.  919KLNS. 


«•  END  OF  L18T1N0  *• 


APPENDIX  H 


SIMPLIFIED  ANALYSIS  FOR  A  LASER-INDUCED  MACH  WAVE 

The  esse  of  a  thermoacoustic  source  moving  8t  transonic  velocity 
deserves  special  attention  because  it  yields  very  high  pressurr  transients. 
In  this  appendix  we  seek  an  analytical  expression  for  the  pressure 
waveform  radiated  by  an  MTS  moving  at  velocity  v=c/sin0ocos$o  towards 
the  receiver. 

In  order  to  simplify  the  analysis,  it  is  assumed  that  the  sound  field 
r8di8ted  by  the  transonic  source  is  the  same  8s  the  one  radiated  by  a 
stationary  source  excited  by  a  short  impulse  of  light.  The  reason  for  this 
similarity  lies  in  the  fact  th3t  all  the  acoustic  disturbances  emitted  by  a 
source  moving  8t  the  speed  of  sound  arrive  nearly  simultaneously  at  the 
fixed  receiver,  so  that  the  receiver  sees  in  fact  the  result  of  an 
instantaneous  input  of  heat  in  the  medium.  This  approximation,  which  is 
valid  only  in  the  ferfield  so  that  the  angle  of  observation  between  the 
source  end  the  receiver  remains  approximative^  constant,  has  been 
discussed  in  the  Soviet  literature  .WJ2,26 

It  has  been  shown  (Eo.  3.4)  that  the  acoustic  pressure  seen  at  the 
receiver, is 


pU)  =  hL(t)  *  he(t)  *  lt(t) 


(HI) 


AD-A160  900  GENERATION  OF  UNDERHATER  SOUND  BV  A  MOVING  HIQH-POHER 
LASER  SOURCE(U)  TEXAS  UNIV  AT  AUSTIN  APPLIED  RESEARCH 
LABS  V  H  BERTHELOT  81  AUG  S3  ARL-TR-83-21 
UNCLASSIFIED  N00014-82-K-0425  F/G  20/1 


where  the  asterisk  denotes  a  convolution  and  the  subscript  t,  a  time 
derivative. 

Assuming  that  the  effect  of  moving  the  source  at  transonic  velocity 
is  simply  to  shrink  the  effective  laser  pulse  as  seen  by  the  receiver  at  the 
retarded  time  r0/c,  to  a  very  short  impulse  of  light,  we  rewrite  Eq.  (HI)  as, 

p(t)  ~  hL(t)  *  h0(t)  *  8t(t-r0/c)  ,  (H2) 

where  8  is  the  Dirac  delta  function,  it  has  been  shown  in  Section  1 1 1.6  that, 
in  the  farfield,  l\(t)  tends  to  behave  like  the  time  derivative  of  a  delta 
function,  so  that  Eq.  (A2)  becomes 

p(t)  ~  d2/dt2  ( h0(t)  1  (H3) 

The  acoustic  pressure  received  at  the  receiver  is  therefore  proportional  to 
the  second  time  derivative  of  the  impulse  response  of  a  cross  section  of  the 
laser  beam.  Assuming  a  Gaussian  intensity  distribution,  we  can  express  the 
received  pressure  as 


p(t)  ~  (2t2-l)  exp(-t2)  ,  (H4) 


where  t  *  (r0-ct/a  sin$0).  Taking  the  Fourier  transform  of  Eq.  (H3)  gives  on 
estimate  of  the  spectrum  of  a  laser-induced  Mach  wave.  For  a  Gaussian 
intensity  distribution,  it  takes  the  form 


p(u)  ~  Q2  exp(  -Q2/4) 


where  Q  *  f0  /  (e  sinfl0/c)  is  a  nondimensional  frequency  characteristic  of 
diffraction  effects  imposed  by  the  finite  dimension  of  the  laser  beam  cross 
section. 

This  lost  result  is  of  practical  importance  as  far  as  long  range 
propagation  of  thermoacoustic  signals  under  water  is  concerned,  because  it 
shows  that  the  acoustic  signal  contains  some  low  frequencies  which  will 
propagate  long  distances  under  water  before  being  absorbed. 


APPENDIX  I 


DIRECTIVITY  PATTERNS  AND  SPREADING  CURVES  WITH  A  RUBY 

LASER 


In  this  appendix  we  present  some  experimental  results  that  were 
obtained  with  a  ruby  laser.  The  optical  coefficient  of  absorption  of  ruby 
light  (optical  wavelength  0.6943  pm)  in  the  water  tank  was  measured  to  be 
1.5  Np/m  (see  Appendix  F),  so  that  the  thermoacoustic  source  was  fairly 
long. 

Figure  42  shows  the  directivity  pattern  of  such  a  source,  measured 
in  its  vertical  plane,  at  a  radius  of  4  m  from  the  initial  point  of  impact  of 
the  loser  beam  on  the  surface  of  the  water.  The  modulation  frequency  was 
9.75  kHz  and  the  source  velocity  was  v=0.5  c  ,  where  c  is  the  speed  of  sound 
in  the  water  (1486  m/s).  The  experimental  data  on  Fig.  42  is  compered  with 
two  theoretical  models.  The  solid  line  in  Fig.  42  represents  the  analytical 
expression  given  by  Lyamshev  and  Sedov  in  reference  18.  In  this  model, 
which  is  valid  only  in  the  farfield,  the  motion  of  the  source  is  taken  into 
account  by  a  simple  ad  hoc  Doppler  shift  of  the  modulation  frequency.  The 
dashed  line  in  Fig.  42  was  obtained  from  the  numerical  model  described  in 
Chapter  IV.A.  It  can  be  noted  that  this  model  seems  to  predict  the  presence 
of  sidelobes  in  the  directivity  pattern,  whereas  the  prediction  by  Lyamshev 
and  Sedov  does  not  show  this  feature.  The  reason  for  this  difference  is 
that,  contrerily  to  our  numerical  approach,  their  model  is  not  valid  in  the 


nearfield.  The  forfield  prediction  for  the  tilt  angle  of  the  mein  lobe  is  ebout 
89°,  whereas  the  experimental  result  indicates  a  value  of  about  72°.  The 
numerical  model  predicts  a  value  of  82°,  and  although  it  yields  a  better 
estimate  than  the  ferfield  theory,  it  still  overestimates  the  tilt  angle  by 
10°. 

Figure  43  reports  the  result  of  a  similar  experiment,  the  major 
difference  being  that  the  source  was  moving  at  Mach  1.  The  modulation 
frequency  was  19.5  kHz.  Again  the  nearfield  directivity  pattern  shows  some 
sidelobes  structure  end  a  very  large  tilt  angle  of  the  main  lobe,  most  of  the 
acoustic  radiation  being  radiated  at  20°  from  the  interface. 

Figure  44  also  shows  a  vertical  directivity  pattern,  but  this  time  the 
source  was  moving  at  supersonic  velocity  (M=1.5).  The  modulation  frequency 
was  about  29.25  kHz.  The  distance  of  observation  was  9.4  m  so  that,  in  this 
case,  the  receiver  was  far  enough  from  the  source  to  ignore  nearfield 
effects.  Experimental  results  show  that  indeed  the  main  acoustic  lobe 
narrows  significantly  to  only  a  few  degrees  for  the  half-power  beamwidth, 
and  approaches  the  interface  between  air  and  water.  This  is  qualitatively 
confirmed  by  the  numerical  model  (program  MTS)  but  there  is  an 
unexplained  significant  quantitative  difference  between  the  predictions  and 
the  experimental  data.  A  possible  explanation  for  the  presence  of  some 
irregularities  in  the  experimental  beam  pattern  could  be  the  presence  of 
some  spurious  harmonics  in  the  signal. 

Experimental  data  for  a  horizontal  directivity  pattern  are  shown  in 
Fig.  45.  The  radius  of  observation  was  2m,  the  modulation  frequency  was 
9.75  kHz,  and  the  source  Mach  number  was  0.5.  As  discussed  in 


FIGURE  43 

VERTICAL  DIRECTIVITY  OF  A  TRANSONIC  MTS 
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Section  IV.B.2,  there  is  a  strong  reinforcement  (more  than  10  dB)  of  the 
acoustic  radiation  in  the  direction  of  motion  of  the  source. 

Finally,  the  sound  level  dependence  on  range  of  observation  was 
determined  for  an  MTS  moving  at  Mach  1.5  towards  the  receiver.  The  results 
are  given  in  Fig.  46.  The  Doppler  shifted  received  frequency  was  about  61 
kHz  at  all  ranges.  The  results  confirm  that  close  enough  to  the  source,  a 
slope  of  3  dB  characteristic  of  cylindrical  spreading  can  be  expected.  But 
when  the  receiver  is  far  enough  from  the  source,  the  classical  spherical 
spreading  lew  of  6  dB  per  doubling  of  distance  dominates. 

The  experimental  data  presented  in  this  appendix  were  obtained  with 
a  "mechanical  modulator"  in  order  to  achieve  the  modulation  of  the  laser 
beam.  The  optical  modulation  system  (Pockels  cell  and  polarizers) 
described  in  section  II.A  became  ineffective  because  one  of  the  polarizers 
was  damaged.  A  mechanical  modulator  was  built  to  replace  the  optical 
modulation  system.  The  mechanical  modulator  consists  of  a  series  of 
wooden  sticks  placed  evenly  over  the  water  surface,  along  the  laser  beam 
path,  so  that  the  laser  beam  could  penetrate  into  the  water  only 
intermittently.  The  thickness  of  the  sticks  was  equal  to  the  spacing 
between  the  sticks,  end  this  spacing  was  equal  to  the  laser  beam  diameter 
on  the  surface  of  the  water.  This  modulation  system  turned  out  to  be  very 
effective,  although  the  modulation  frequency  was  dictated  by  the  laser 
beam  scanning  velocity.  The  relationship  between  the  modulation  frequency 
8nd  the  source  velocity  W8s,  in  our  design,  given  by  f0  =  19.5  M  ,  where  M  is 
the  source  Mach  number  and  where  f0  is  expressed  in  kHz. 


FIGURE  46 

SPREADING  OF  AN  MTS 
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